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ABSTRACT 

The basic Mitrovic's method for analysis and design of linear feed- 
back control system is studied. Some further properties of the basic 
Mitrovic's equations and curves are investigated and applied to analysis 
and design of linear feedback control systems. The method developed in 
this paper gives more flexibility to the basic Mitrovic's method. 
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ts Introduction. 

The Mitrovic's method is an analysis and design technique of linear 
feedback control systems using algebraic equations and graphical methods. 
The method is based on conformal mapping of the S-plane into the real 
domain through the characteristic equation of a feedback control system, 
and the basic principles involved in the Mitrovic's method are derived 
by applying the basic theorem by Cauchy in the mapping procedure. 

The fundamental theories are adequately developed and applications 
of the method are demonstrated in references 1 and 2. The main advantages 
of the Mitrovic's method seem to lie in the fact that all operations are 
carried out in the real domain and that it has potentiality of controlling 
or adjusting multiple variable parameters. 

The analysis of a system in the Mitrovic’s method is done by analyz- 
ing the relative location of the Mitrovic's working point M with respect 
to Mitrovic's curves, and design of a system is to reshape the configura- 
tion of Mitrovic's curves and the M point associated with the system of 
interest so as to meet the design specifications. The analysis and 
design techniques are deduced from the fundamental properties of the 
BQ vs By, curves for $= 5, denoted by 64 with the M point as listed 
below. 

1) If the M point is located at sari on the Bo vs B, curve 

for S= %, them the characteristic equation associated with 
the Ist curve and M point has a pair of complex conjugate 
reots Ts-%,tn,+ |- 3? Ont 
2) If tangent lines are drawn from the M point to the Bo vs By curve 


for S= 1.0 (the leeaneeye then the negatives of the slopes of 


the lines are the real roots of the characteristic equation. 


1 





3) If the M point is located on a region enclosed by a ;, curve 
encircling in the clockwise direction, and if the M point is 
in the first quadrant on the By vs By plane, then all of the 
roots of the characteristic equation have 5% greater than “1, 

In this paper further properties of Mitrovic's equations and curves 
on the Bo vs By, plane are investigated to extend the method to problems 
which are not treated in the literature, and to give the method a better 
flexibility. 

In Section 2 the fe curves are studied to provide a procedure to 
evaluate frequency response of closed and open loop systems on the 
Bo vs B, plane so that one can obtain information for frequency response 
and time response simultaneously. Constant band width curves are intro- 
duced on the By vs By plane to design a system to meet the band width 
specification precisely with other specifications such as S% and Kv 
specifications. 

In Section 3 fourth order system curves are studied to provide 
fourth order system charts to fit any fourth order system. The analysis 
and design procedures developed in references 1 and 2 can be worked out 
with preconstructed curves on the charts and also a different approach of 
design technique is tried with fourth order charts. 

In Section 4 an analytical design technique is derived so that the 
labor required in constructing Mitrovic'’s curves can be avoided or mini- 
mized. 

Thus each section can be regarded as independent from another; 
however, the principles developed in Section 2 may go with any of the 


sections to give information on frequency response. 






Be Frequency responses of closed and open loop systems on the Mitro- 
vic's Bo vs By plane. 
Mitrovic's curves for S= 0, the lo curves, are obtained by mapping 
the imaginary axis of the S-planme in to the Mitrovic's plane. 
For a characteristic equation of a system 


EQS + aes Tome. + aye a,— 0 (2-1) 


Mitrovic's B, vs B, curve for 0& C641 is defined by 


k-1 
= ae. 
24 (2-2) 
= Pe (5 oe 


where @, and © are the natural frequency and damping ratio 
respectively of a root of the closed loop system, and 
@($) is a function of $ only. (See Appendix I.) 
which is a curve on the Mitrovic's Bg vs By plane independent of the last 
two coefficients ag and a) of the characteristic Equation (2-1). 
Mitrovic's working point M on the Bo vs B, plane is defined by the 


coefficients ay and ay such that 


M=M (4, , 40) 
thus, for a given characteristic equation, a set of Bo vs B, curves and 
an M point are uniquely determined and conversely. 
For a particular value of damping coefficient C= 0 the equation 
defining the [> curve is, (using proper values of B.(0) and replacing 
Q),, by @ the frequency on imaginary axis of the S-plame in Equation (2-2) 


Bee Ge acw- Sn = 
B= (2-3) 


Aw a,a* See eee, 


Because of the nature of the lo curve, the frequency dependence of 
Equation (2-1) can be represented by the lengths of vectors from the M 
point to the in curve. This indicates that frequency response of a 
closed loop system can be evaluated on the Mitrovic’s plane. 

In this section relations between the frequency response (magnitude 
and phase) of a closed loop system and the lo curve with the M point on 
the Mitrovic's Bp vs By plane are investigated. Band width of a closed 
loop system is evaluated and the locus of the M points for a system to 


meet the specified band width is constructed. 








2-1. Interpretation of equations (magnitude of frequency response) on the 
Bo vs By plane. 


Consider a closed loop transfer function of a feedback control 





system 
e NS) M(S) 
ee S = - i) a= pa a eee 2-4 
aC mans PC rae eer: c ee): 


where N(S) is a polynomial of S of order less than n. 

Note that F(S) = 0 is the characteristic equation from 

which the associated ie curve and M point are determined. 
The magnitude of the frequency response of the system represented by 


Equation (2-4) is 


A=|§ Gal a (2-5) 


[Ref Fie) +j ont FU)}| 





Re{ F (jw) } and Im{ FG} can be represented in terms of a, By and 


By respectively where B, and Bp are the Mitrovic's variables defining 


ay: 
the iF curve which are functions of W defined by Equation (2-3). 

Rewriting F(S) in Equation (2-4) as the sum of even F(S) and odd 
F(S), 


F(S) =(4,t 48+ a,S°+ ++ )+S(A, +4;5° + asS% +++) 
then FGUW) = Ref FGw)} +j ImfFlja)} 


= g F a - y 
=(a,-(a,W-4,W +r + jala, —( a, aeW + = as) 
therefore by using Equation (2-3) 
Ref F(jw)! = a,- (4,0 -aaw + ae =) =p -—Bo 
Ton{F Gay} = ol Aa aetea cui ae mI = w(a, ~ 2B) 


thus, Equation (2-5) can be written as 


jan) ihe (2-6) 


he |(4o-By) +) (a,-B) | 











The interpretation of Equation (2-6) in By vs By plane is as fol- 
lows. The magnitude A at frequency ® is the ratio of the length of two 
vectors N(j@ ) and (ay - Bo) +/ O(a, ~ By). 

Since a, and ap are the coordinates of the point M(a, 54) on the 
Bo vs By, plane, and By and Bg at any frequency are determined by the co- 
ordinate of a point at @ on hs » the numerical value of the denominator of 
Equation (2-6) can be determined easily. 

Suppose the le curve and M point configuration is as shown in Fig. 
2-1 and choose a point P on the [5 curve at which the frequency is 0, 
then the coordinates of P are (By (@ ), Bg(@)). Locate a point Q at 


(a), Bo(w)) so that Q is at the same ordinate as P and is on the Line 
By = a, , then the distances from M to Q and Q to P are 

MQ =|4.- B,a)| ‘ Op =la,-B,(w)| respectively. 
Locate a point R on the line segment QP such that QR=W-QP and asso- 
ciate a frequency ® to the point R, then QR= w-2P = w|a,-B,(w)| 
and MR=(MQ+4e ) = |Ca.-B(w)) +f (a,-B,(o)) | whieh is the 
value of the denominator of Equation (2-6) at frequency 4). 

As @) varies, the point P moves along the IF curve and the point Q 
moves along the By, = 4 line maintaining the same ordinate as P. The 
point R will describe a curve which is represented by the dotted curve 
and denoted by the G) (a, - By) plot. The O(a, - B,) plot is a well defined 
curve if the Ib curve is well defined and has the following common proper- 
ties for all systems: 

1) It has frequency on the curve as a parameter such that two 

points constituting a pair of points one on the i curve and the 
other on the 0)(a, -B,) plot have the same frequencies with the 


same ordinates. 








Bo 


we 
ols ~B_) milo 
oA a! 


a 


"C,, (Magni tude circle) 


M (a ia 


ue 
Oo 
ae 





als 


Fig. 2-1 Illustration of obtaining Frequency 


i] 
response for a-syStem with no’ zeros. 





2) It starts (or the originating point is} at (a, ,0) and passes 
through the intersection point of the |; agra sane the 
B, = a, line. : 

3) The points for @= 1 on the |, curve and O(a, ~ By) plot are 
identical. 

4) It lies between the ik curve and the B, = a) line for W<1 
and dually. 

If any circle centered at the point M passes through a point at 


@ = @, on the @(a, - By) plot with radius of r, then r, represent the 


1 
value of the denominator in Equation (2-6) at @=@,. The circle cen- 
tered at the M point is called the magnitude circle. The numerator of 
Equation (2-6) can not be evaluated from the [5 curve as the denominator. 
It is possible to draw separate Mitrovic's curve for N(S) and obtain 
N(j®) but this is labirious if N(S) is of high order. 

The simplest case is for NC 5 ) independent of frequency. This 
is the case for a system that has no zeros in its closed loop transfer 


function. 


Suppose a closed loop system has no zero, then 
a 
A=; —~ 
|(4,-B + jw (a,-B,)] we 
ay-B,) +j (4,-B,)| = 207A 
The magnitude A, at W= 0) is ao/Ty and conversely if a magnitude 
A= A, is given then the frequency to give A = Ay is found by drawing a 
magnitude circle with radius equal to ag/A,> and locating the point or 
points where the circle intersects with (a, = B,) plot. 


Some particular points in the frequency response of a closed loop 





system such as the resonant point, zero db point and half power point 
(band width) can be found very simply. 

Drawing a magnitude circle tangent to the Wa, - By) plot gives 
the resonance point. Magnitude circles with radii of ag and (2 ag give 
the 0 db point and band width respectively. 

For a system with N(S) in Equation (2-4) of order other than zero, 
the procedure is slightly modified but still the basic idea is the same. 

This is illustrated with Fig. 2-2. Modification must be made to 
take care of the frequency dependence of the numerator in Equation (2-6). 
Since the denominator in the equation defining the magnitude of fre- 
quency response bears the same form regardless of the order of N(S), the 
construction of the w(a, - Bi) plot is the same for all systems. 


Rewrite Equation (2-6) as 





(2,- By) +5 (4,)- Br) = JING @)) (2-7) 
and interpret as follows; 

The quantity at any on the left side of Equation (2-7) is deter- 
mined by drawing a magnitude circle, and it is represented by the length 
of radius of the magnitude circle. On the right side of Equation (2-7), 
| N (jo) | has some finite value for every @. Then @ number which 
multiplies to | NGjo) | to give the equality in Equation (2-7) is the 
reciprocal of the magnitude. 

In Fig. 2-2, the value of IN G@)| is represented on a line drawn 
from M point in an arbitrary direction. This line is denoted by In (joy 
and at each point on the line, there associate w such that if the point 


N is at &) on the line then 


MN =|N Goo)| and A=MN/T 








lo 


{ 

| 

2 S03, =B, ) plot. 
74 





Fig. 2-2 Illustration for obtaining frequency response 


for a system with zeros. 
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Conversely, Wz: a4 for given magnitude A = A, is determined by search- 


1 
ing such @ that for two magnitude circles with radii of r and Ar cross 
the @ (a, - B,) plot and the |ncjory| line respectively at the same fre- 
quency. This can be done by a few tries and errors. 

The arguments developed above are not restricted by the order of 
system equation. The only restriction is to have i” curves plotted on 
the Bg vs By, plane with the same scales of B, and By axes to provide 
simplest graphical solutions. Bo vs BL curves can always be plotted on 
such axes scales by using suitable frequency scaling technique. 

As mentioned previously, the higher the order of the numerator 
polynomial of a system, the more complicated the method is. The problem 
is to provide a convenient method to represent the frequency dependence 
of the numerator of a closed loop system equation. 


In the following paragraphs, such methods are presented for the 


cases of first, second and third order numerators. 


11 





2-2, Representation of frequency dependence of N(S) on the Bo vs By plane. 
In the previous section, it is seen that the frequency dependence 
of the numerator polynomial N(S) in a closed loop system equation must 
be represented as a distance from the M point. 
In Fig. 2-3 it is illustrated for N(S) of first order. When a 
closed loop system has one zero then N(S) is of the first order as 


N(S)=bS +a, and IN Gao)|=|a,+j bo 


The real term of N(j® ) is independent of frequency. By construct- 
ing a rectangle with a fixed side of length equal to ay and another side 
which is in right angle with the fixed side varying linearly with W , 
INGoo| can be represented by the length of the diagonal. 

In Fig. 2-3, a point T is located on the Bo = a line such that 
MT9 = b)- If the line segment MT, is calibrated linearly with ® such 
that 4) = 0 at M and @= 1.0 at To; and if the same calibration is ex- 
tended beyond the point T>, then at each point on the line Bo = ay to 
the right of M, there associate © linearly with the distance from M to 
each point. For any point T; at W= Qj, MT, = by ®,- Noting that 
TT, = a> the length of diagonal of the rectangle MT, Ta, is equal to 

Fit =|2.=j burl [un geon| 
For a system with two zeros in its closed loop system equation, N(S) is 
of second order polynomial as 
NGS)EBSthS+dy and NG ake |(ao- bak) +j 60 | 
The imaginary term can be represented in exactly the same way as in case 
of first order N(S). 
The real term is now also frequency dependent, and can be repre- 


sented on the B, = 4) line with quadratic calibration. 
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TO} a=] 


J [w(sq) = |(ag-b pts )+30 00 | 


aT 





| 
i 
I 
{ 


1 
Fig. 2-4 Representation of [1( 50) | =fag=b,6°)+Jb,co] 
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In Fig. 2-4, b;® is represented as in Fig. 2-3. To represent the 


Oo the line B, = ay is calibrated as follows. The point 


Poa, 0) is labeled G)= 0. From the point Po up éach point on the 


term (a) - b 


By = a; line is calibrated with ® such that if a point P, on the By = a) 


1 


line reads ()= a4 then the distance PoPy = byw, Neting that 
ME, s 2,» 
MP =MR-BR = 4o— bau 
By forming a rectangle MT,|TP, ; the length of the diagonal MT is 
MT =|(2.-bw}+j bw = |v Go) 
For a system with N(S) of third order or higher, Inga| can be repre- 
sented in a similar way as in case of second order N(S) with some modifi- 
cations as the case of first order N(S) modified to account the second 
order N(S); however, it requires laborious calculations. For example, 
for third order N(S) 
NAS) Spc ce eet aa (2-8) 
NC; ) = (a,-b.w ) + j(ba- bs a) 
which indicates that the Bg = 4 line must be calibrated to acccunt 
(b, W re) which requires cubic calculations, and quadratic calcula- 
tions for fourth order N(S) and so on, 

Those calculations can be avoided by applying Mitrovic's basic 
equations for the y curve to Equation (2-8) and by proceeding the same 
way as was done to represent F( jo ) by means of W(a, - By) plot. 

This is illustrated with Fig. 2-5. 

Since N(S) is of third order, the [} curve obtaimed from Equation 
(2-8) is a straight line and this is denoted by ee - Let Mitrovic's 
variables associated with Equation (2-8) be Biy and Boy. Then 
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Fig. 2-5 Representation of [*(3e)[=|(a,-R, )+ju(d, -Byy ) 
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ie is defined by 
EB, bos es 
Bov= hos - 

The point My(b1 > ay) which is analogous to the point M(ay » aq) for 
the characteristic equation F(S) = 0, is determined by the coefficients 
by and ag of Equation (2-8). 

If lin and Boy @xeS are chosen identically with B, and Bg axes, then 

has a slope of bo/b3 and My is located on a horizontal line passing 
through M on the By VS By plane. 

If Wb; - Byy) Plot is constructed from i” and the line B, = Ein = 
b;, then by the same reasoning as for @ (a, - By) plot, for any fre- 
quency @, the quantity ING 0); ) | = oe Boy) +j @y CoB) | 
where Byy and Boy are coordinates of a point at Ws Q on im Ls 
equal to the distance from My to T, where the point Ty is at We , on 
the 4) (db, - By y) plot. 

Since | nj >| and (ag - Bg) +) @ (ay - B,)| must be compared to 
get A the magnitude, and since N(j@ ) is measured from My whereas the 
latter is measured from M, it is necessary either the a)(b, = Eiw) plot 
is translated horizontally to the right or the W (a, - By) plot to the 
left so that My and M coincide. In Fig. 2-5, the G)(by - Byy) plot is 


translated to the right by amount of a, - by » then 


Ngwp|= Mz =MT 
The method representing Incjo >| for N(S) of up to third order is 
presented. For high order N(S), |ncj@ >| can be represented by an analo- 
gous method as in case of third order N(S); however, it is inevitable to 


be laborious to represent | Ncjo as the order of N(S) increases. 
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On the By vs B, plane for feedback control systems with adequate 
relative stability, the M point is usually located in che region near to 
the origin as compared with whole region covered by the ie curve, This 
is due to the fact that [; curves for higher §$ are biased towards the 
origin. Therefore in obtaining frequency response usually a first part 
of the Bs curve is required. 

It must be noted that, since the basic Mitrovic's equations are not 
restricted to the characteristic equation of a feedback control system 
but apply to any polynomial with positive coefficients, the principles 
developed above apply to any transfer function which is expressed in the 
form of a quotient of two polynomials provided the Bo vs By curve for 
~$ = O and the M point associated with the transfer function are defined. 
Therefore it is obvious that the magnitude in the frequency response of 
an open loop system can be evaluated in a similar way on the same 
Bo vs B, plane where the magnitude of the closed loop system is evalu- 
ated. 

In general the open loop frequency response can be obtained by con- 
structing the i curve and by locating the M point associated with the 
open loop transfer function and then by proceeding exactly the same way 
as for the closed loop system. In Section 2-7, some examples for evalu- 
ating the phase of the frequency response for both open and closed loop 


systems are shown, which will well illustrate the above arguments. 
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2-3, Determination of band width of a closed loop system on the 

By vs By plane. 

For any given system, band width defined at magnitude of 1//2 or 
half power point can be determined precisely by the method described in 
Section 2-2. 

For a system with no zeros, band width is found at the intersection 
point of QM (a, - B,) plot and a magnitude circle with radius equal to 
[2 ay. 

For a system with zeros, band width can be found by trying a set of 
2 concentric magnitude circles with radii of r and ¥2 r_ such that the 
two circles cross the @)(a, - B,) plot and [ncjo >| respectively at a 
same value of frequency. 

Frequently in the study of feedback control system the point of band 
width alone of the frequency response is of interest and also sometimes 
upper bound and/or lower bound of band width is of interest. 

By the property of the equation defining band width and associated 
geometry on the Bp vs By plane it is possible to set up upper and/or 
lower bound of band width for a given system. Among the various points 
in a frequency response curve of a given system, if band width only or 
approximate value of band width is required, then by setting up the lower 
and/or upper bound of band width, the effort in plotting (a, - By) can 
be economized or even without plotting any portion of W(a, - By) the 
approximate value of band width can be determined. 

In this section the three cases for a system with (1) no zeros, 

(2) with one zero, (3) with two and three zeros are considered. 


For a system represented by Equation (2-4), and using Equation (2-6), 
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the band width defined at half power point is the frequency to satisfy 


the equation 


Paris ING) nen 
IS jw\|- |(Ap- Bo) +j@ (d1-By)| 2 


or | (a,-B,) +j@ (4,-B,>| = /2|NGo)| (2-10) 


All that has to be done is to interpret Equation (2-10) on the By Vs By 
plane. The left side of Equation (2-10) can be interpreted by the 
W® (a, - By) plot as shown in Section 2-1, and for the right side is 
needed a modification to the method described in Section 2-2 to account 
for the factor fae. 

Consider a closed loop transfer function with no zero. 


Let G) be the band width, then from Equation (2-10) 
| (@p~By(w,) +j@,(4,-B(@,)|=/2 do (2-11) 


must hold, where B,( @,,) and Bo( @,) are the coordinates of a point on 
the Is curve at W= Q,-> 
Equation (2-11) implies that a magnitude circle with radius of 

Pez ag intersect with the Ga; - By) plot at @, which is described in 
Section 2-2. Since the ly curve tends to infinity in some direction as 

Q@ tends to infinity, Wa, - B,) plot also tends to infinity. Noting the 
finite radius of a magnitude circle and that the a)(a, - By) plot starts 
from a point inside the magnitude circle with radius {2 4, it is seen 
that there exists at least one intersection point between the circle and 
the W (a, - B}) plot. Depending upon the relative location of the M 


point with respect to the D curve and also depending upon the geometric 


shape of the lh curve there may be more than one such intersection point 
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implying more than one band width to exist, 

On the By vs By plane with the if curve, by drawing a magnitude 
circle of radius DS a and by inspecting the values of @ at all inter- 
section points of the circle with the Ip curve if there are any, and by 
estimating the shape of the @ (a, - B,) plot, information about the band 
width can be obtained. Also the point at W=s 1 on i and the points 
where the if curve intersects with By = ays Bo = a and Bp = (1 + (2)a 
lines provide additional information. 

Suppose a system is represented by the I curve and M point as in 
Fig. 2-6. A magnitude circle Cap with radius of V2 ag is drawn. The 
points P,, Py, P3, Py, and Ps represent the points where the Ie curve 
intersects with Cap : Bo = a » By =- @ +/2)a, or By = a, - 

Suppose at each point of those intersections, frequencies are read 
to be > Ws » @,, @, and OF as shown. 

The latter part of the he curve is neglected since it is assumed 
that for that part of re > the W (a, - B,) plot stays outside of Cap, 
In an actual situation this can be judged by looking at the quantity 
Ge (a, - B,(@,_)) and @7(a, - B (@7)) where @, and W, are the fre- 
quencies at the points where the latter part of e intersects with 
By = (C 1b Se {2)aq and Bo = ap respectively. 

The O(a, - Bi) plot resulting from the portion of the b curve 
above the line Bo = ( L+¥2 ay can never cross Ca, , therefore, then 


there are four possibilities for band width as 


(1) W3= =), Ce), = oe 


(2) 0,= Op= a, (4) dip =) 
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Fig. 2-6 Determining approximate band width 


for a-system with no zeros. © 
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Each of the four possibilities above can be evaluated by looking at the 
point at @)= 1 on the ik curve. | 

Suppose ®3>1 then the ®(a, - B,) plot for Q= ; Stays outside 
of the 5 curve in turn outside of Cy,, indicating W,< 3. 

Suppose @,<1 then Q3P3 >, (a, - B (@,3)) and the W (a, - By) 


plot must cross C,, somewhere between OW. and W,, therefore 


Ab 
CD fla Fo (6 or On Op 1a, 

Suppose @-= 1 is located between ©, and W, or W,< 1<@3 5 
then, since the i curves for Dy < W<W, are inside of Cap » the 
Ww (a, 8) plot must cross Cy, between “= 1 and Ds » therefore 
1l< w.,< 03° In this case there is still a possibility that w,< W, 
but this can be checked by looking at he quantity G)(a, = By) for some 
® less than Wb . 

Suppose a<i< Os » then the a) (a, - B,) plot stays outside of 
Cap for 1S W<W, and W? W3 The portion of the @(a, - B,) plot for 
Bo < W< W3 may be inside of Cy . If this is the case then Gy < Gy,< Qs s 
again the possibility that a, < Wy must be checked. If Wa, - B,) for 
Wo < 0 <0 has no point inside Cap then da; - By) for 21 is out- 
side Cab and G<1 . If the points Py and Py are very close, then 
Mpy=1 = W, = We 

Suppose a, > 1, then there is no possibility that w), is in the in- 
tervals (1, W 2) and W2W,. Since the @(a, - B,) plot at @W= 0 is 
at the point (a,, 0) and the O(a, - B,) plot at @z 1 is outside of 


Cc band width is definitely o< ai< due Wy and possibly 


Ab? 
O)9 < a, < 003 


If either one of the frequencies w, and 6), is equal to 1, then 


22 





Q), = 1 and other possibilities can be investigated as above. 

Thus the four possibilities listed above are saeiyeed! The geometry 
involved in Fig. 2-6 is so simple that the relative locations of 
O(a, - B) plot for OO > Wy, » O3S Wa), , Gy < GS 3 ; a, < RS a, 
and Os W<a), with respect to Ca, are determined almost by inspection 
using the point Ws: 1 as a guide, resulting in determination of upper 
and/or lower bound of band width. 

It must be noted that, depending upon the system of interest, there 
will result various configurations of the ys curve and the working point 
M (a) 289) however, the basic principle of reasoning is the same as 
above. 

For a system with one zero, the band width defining Equation (2-10) 
is to have frequency dependent terms in both sides of the equation and 
requires some modification in the method applied for a system with no 
zeros. If a system has one zero in its closed loop transfer function, 
then the equation defining band width is 

|(a5-By (9) + j Wy ( a,-B(O,n|=/2 la, +} Py | (2-12) 
where b, is the coefficient of the s' term in the numerator of 
the given closed lcop transfer function. Bo( ®,) and By (@,) 
are the same as in Equation (2-11). 

Suppose the le curve and M point configuration is as shown in Fig. 
2-7. To account for the frequency dependent term {2 b,® on the right 
side of Equation (2-12) the line Bg = ( 1+/2) aq (or Bo = ag line) 
is calibrated with @® linearly such that the points Q and Ty correspond 


to @s=: 0 and G= 1 respectively, and the length QT; is equal to 


Vz by - Then the distance from Q to any point at W= 0, on the Line 
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Fig. 2-7 Determination of band width for a system 


with one zero. 





( ie a, will represent V2 byw4 and it is obvious that 


wo 
© 
u 


ES 


JZ | a,rjb.o | which is the quantity on the right side of Equation 
(2-12) for @=) 

Suppose the O(a, - B,) plot is as shown by the dotted curve in 
Fig. 2-7. Noting the nature of the W(a; - By) plet that it always 
starts at the point (a, > 0), crosses i at = 1 and at the point 
where iy intersects with B, = a; line, it is easy to guess or sketch the 
O(a, - By) plot. 

The band width @, is some | such that Equation (2-12) holds. 
This can be interpreted in Fig. 2-7 as follows: Suppose a magnitude 
circle C,,? is drawn such that it pass through the point W at =o’ 
on the B, = (lL +V¥ 2) a,line and suppose it reads Ws Ww" 
on the oa, - B,) plot. If w':w” then Equation (2-12) holds at 
frequency equal to a » therefore ®, 200i 

To determine band width one must find a Cy such that it crosses 
both O(a, - B,) plot and By = ¢ 1+/2 Jap at the same frequency. 

For a Cary if w+ oo” then the two frequencies are either w> a 
or W<W” , By drawing another Cy.” until the inequality between 
a) and a is reversed from the previous set of ® and aw” the bounds 
of band width can be set up. 

For example, for the configuration of curves as shown in Fig. 2-7 


W, <1 is assumed. 


For the two magnitude circles Cao and Cay 


Frequency on Bg = (1 +1 2)apy line Frequency on Wa; - By) plot 
Cao d= 0 < O = Wo _ 


Cal = > w= a 








in which the inequality sign is reversed, therefore O< a,< Wr< Das || 
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It is obvious in Fig. 2-7 that there is no possible band width for 
W< Dy . There may be another band width for W> 0, This can be in- 
vestigated by trying some magnitude circles. 

C,0 in Fig. 2-7 is identical to Cab in Fig. 2-6. If by s 0 in Fig. 


2-7, then Cy, Stays on Cro for all QW, which is the case for a system 


it 
with no zeros and agrees with Fig. 2-6. It is readily seen from Fig. 

2-7 that the effect of zero and forward gaim of a system on band width 

can be analyzed. 

For two systems, one for no zero and the other for one zero, if the 
characteristic equations are the same, then band width of the first system 
is less than that of the second. This is confirmed in Fig. 2-7 as for 
the first system G * Wg and for the second system W), >Wo - 

The effect of varying forward gain and the zero on band width for a 
system with one zero is seen from Fig. 2-7 by writing the numerator of the 
closed loop transfer function as 
N(S) = b)S+ap_ = KS + KZ 


Increasing either K or Z is to expand both C,, and Cal which is to 


AO 
increase Op and @, resulting in increased band width. Since K and Z 
do not enter into any coefficient other than a, OF a of the character- 
istic equation, the ty curve remains unaffected by either K or Z. 

Suppose K is fixed and 2 is variable. Since a, amd b,; are independ- 
ent of Z, the @(a, - By) plot and the lines B; = a, , By = a; + v2 by 
remain unchanged. Varying Z moves the M point and the line Bg « 49 ver- 


tically, thus increasing Z expands Cyg and Ca) resulting in increased 


band width. 
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Suppose Z is fixed and K is variable, since aj. a and b, are 


0 
linearly varying with K, all of the lines im Fig. 2-7 ahd ene @) (ay ~ By) 
plot are moved. The M point is moved along a line with some positive 
slope and Cag and Cy, are expanded with increasing K resulting in in- 
creased bandwidth. 

In a system which has positive coefficients in both numerator and 
denominator polynomials of its open loop transfer function, b, is always 
less than ay assuming unity feedback. This means that the bandwidth of 
a system with one zero is always less than that of a system which results 
from replacing b, with aye Therefore for 4 system whose le and M point 
configuration as shown in Fig. 2-7 bandwidth must be 

Oy< ip < Wo 
where Wp is the frequency on the a) (ay - By) plot at which the Wa, - By) 
plot and magnitude circle Ca? intersect. Cy 9 is drawn such that it passes 
through Ty where To is located on the Bo = (i+ (2 jay line at the dis- 
tance of ie ay measured from Q. 

Above inequality is useful to set up the bounds of a, when b, is a 
variable and when it is required to locate M point to meet some required 
bandwidth specification. It must be noted that the bounds determined by 
the above inequality become poorer as the difference between b, and ay 
gets larger, 


It is interesting to note that if by = a, then the open loop 


transfer function is of the form, assuming unity feedback, 


aorta 
s)= > SOOO 
&¢ ) s” te BE SD eee + aoo0 


which is a type 2 system, and this system has greater bandwidth than a 
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system whose open loop transfer function is 


7 bs+4, 
GG) Sots oe ---+ LS eee 


For a system with two zeros, bandwidth can be evaluated in a simi- 
lar way as illustrated in Fig. 2-7. One modification is required in 
Fig. 2-7 to account for the frequency dependence of (2 (ag- by"), 


2 term in the numerator of the 


where by is the coefficient of the S 
closed loop transfer function. This is done by calibrating the line 

B, = 4, as in Fig. 2-4 to represent (a - by) a ) by the distance from M 
to the point at @= a, on the B, = 4) line. To account for the factor 
2 , the origin of calibration or the = 0 point is chosen at a point 


on the B, = 4 line below the B, axis so that it is apart from the M 


1 

point by a distance of {2 a). Then the line By = a; fs calibrated with 
2 

@ such that any point at W= W, on the line will measure (2 boy 

by the distance from the Ws 0 point (originating point) to the peint 


where We: Od, + 





For a system with 3 zeros, the quantity ya (ag ~ by a) 4je0(b1 - by") 
which is the right side of the bandwidth defining equation cam be repre- 
sented in a similar way as in Fig. 2-5. The only difference is to ex- 
pand the translated @(b, - By) plot radially from the M point by factor 
of 1D : 

Then for both cases of the above, bandwidth can be evaluated by 
applying the same analysis as for the system with ome zero in Fig. 2-7. 


The following numerical examples illustrate the arguments developed so 


far. 
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Example 2.1 


Obtain frequency response (magnitude) of the closed loop system 


6 
a(S) Se 


The system represented by Equation (2-13) is a third order system. 


(2-13) 


Since it has no zeros the bandwidth will be less than the frequency at 
the point where the ie curve and Bo = Gi ay line intersect. If 
the frequency response is required up to W-: Oy? the portion of the 
ie, curve needed is below the Line By = (l+/f/2 ap = &) 45), 

In Fig. 2-8 the frequency response of the system (2-13) is repre- 
sented on the Bo vs By plane. 

From the characteristic equation of the system (2-13) the i. curve 


is defined by 
B=o 
2 (2-14) 
B)=3® 
which is a straight line with slope ef 3 passing through the origin. 


The magnitude of the frequency response is defined by 
Aj@)= 


nt ees fe eles ES on 
(4,-8,)+j@(4,-8,)| |¢-B,) +jo(2-B,)| 


where BY and By are defined by Equation (2-14). 


The point M (2, 4) is located, and the lines B, = a, = 2, Bo = ag = 4 


(2-15) 


and Bg = (1 +[2 Jay = 9,65 are drawn as shown. The &) (a, - B,) plot is 
constructed as illustrated in Section 2-1. 

The @)(a, - By) plot starts from the point (2, 0) and crosses the 
[Qo curve at G)= 1 and the point (2, 6) at which the I curve and the 


line By = 2 intersect. 


For 0 <@ <1 and for >1.42 the Ma, - By) plot lies to the 
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right of the ( curve and for 1.0< W< 1.42 te the left. The point at 
G) = .5 on the Ww (a, - By) plot is located on a horizontal line passing 
through the point at @)= .5 on the i curve and the distance from the 
line B, = 2 to the point at ® = .5 on the (a, - By) plet is one-half 
of the distance from the line B, = 2 to the point at Ws .5 on the Ve 
curve so that the distance from the B, = 2 line to the point at Ws .5 
on the @) (ay - B,) plot is equal to @(2-B,(W)) Wes which is the 
value of imaginary term at @W= .5 in the denominator of the equation 
defining the magnitude of the frequency response. 

Some points on the (a, - By) plot for some pertinent values of 
frequencies are determined likewise and those points are commected by a 
smooth curve to give the W(a, - B,) plet as shown. Thus each pair of 
two points on a horizontal line, one on the i curve and the other on the 
@ (a; - By) plot has one value of frequency and the distance from the 
point M (2, 4) to any point at frequency @ on the (a, ~ B,) plot is 
the value of the denominator of Equation (2-15). 

Magnitude circles are constructed as shown. Noting that the numer- 
ator of Equation (2-15) is constant at ag = 4. it is seen that the ratio 
of ay = & to radius r of a magnitude circle which passes through a point 
at frequency of @) on the O) (ay - By) plot is the magnitude A (WwW). 

At Q@= 0, the radius of magnitude circle is 4 which gives A = l. 
As the frequency increases from zero, the radius of the magnitude circle 
decreases, in turn increasing A. At Ws 1.18 magnitude circle is 
tangent to the Wa, - B,) plot and has radius r= .75 which implies 

resonance peak My) = Geto 2.33 


at resonance frequency WD, = 1.18 
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For frequencies from 4) = 1.18 on, the radii of magnitude circles in- 
crease, in turn decreasing A. 
A magnitude circle with radius of {2 ay = 5.65 crosses the 
(O(a, - By) plot at @=: 1.76 which implies the bandwidth ts 
Wp = 1.76 
The frequency response up to (= 1.8 is thus represented on the Bo vs 


By plane. 
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Example 2.2 
Given an open loop transfer function of a unity feedback servo as 


___ 42100068 +37) 
G(S)= S(St+h1)(S +1670) 


obtain the magnitude vs frequency plot of the frequency response. 


The closed loop transfer function is 


C (gx mq IOS 4 HE DOD | 
R ue S72 41671.154% 4+ 1226718 + 448000 (2-16) 


Since Equation (2-16) deals with large numbers, it is not convenient to 
work on the Bo vs By plane associated with Equation (2-16). 
By using a transformation or frequency scaling 
S = 100s 
Equation (2-16) transforms to 


£02 (2.1/5 + 4.48 (2-17) 


eS Siac 7a eee 
Then the equation defining frequency response is 


| L.£® +; fast a, | 


sO hy \= 
Re) (496-3, | 0, (/2.27- Bz )| Ca) 


where Bit and B are Mitrovic's variables of B. vs B, curve 


Ot 0 1 
for %S = 0 associated with the system (2-17) and QO. is re- 
lated with @W, the frequency associated with the system (2-16) 
by 
W= 100 , 
In Fig. 2-9 the c curve, (ay, = Bi) plot and magnitude circles 
representing the numerator of Equation (2-18) are constructed. 


The frequency dependence of the imaginary term in the numerator of 


Equation (2-18) is represented on the Bot = 28 line to the left of 
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Bie = 4, line. (Note that the frequency dependence of the imaginary 
term in the numerator of Equation (2-18) is represented on a line dif- 
ferent from that illustrated in Section 2-1.) 

As illustrated in Section 2-2, the distance from the point 


M (12.267, 4.48) to a point at frequency @,on the Bor = 2a 8.96 


Ot 
line is equal to [4.48 +j 12.267 @, | 

The magnitude circles are drawn so that at any frequency the numer- 
ator and denominator of Equation (2-18) can be compared. 

Since the M point is located at the point where the coordinates are 
not convenient numbers, reading off the radii of magnitude circles is 
not convenient. 

This can be simplified by providing a scale which is calibrated 
with the same scale of coordinate axes so that the radius of any magni- 
tude circle is read off easily. This is shown in Fig. 2-9, for example, 


at (), = -5, the value of Inc jo, is read to be 7.5 on the 





Wt 2 45) 
scale provided on the line B,, = 12.267 by drawing a magnitude circle 


passing through the point at 0), = .5 on the line Bo; = 2a and simi- 
larly for the denominator ere - Boe) +) O_ (4,7 By) | Qe 5 6.0, thus 


the following table is constructed: 
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Example 2.3 
Suppose for the system in Example 2.2 only the bandwidth is of in- 


terest. The equation 





4Y8-Bor(My) +/O+ (12.27-B. (w+) mB Ke8 tj/2.10¢ (2-19) 
defines the bandwidth One of the system in Equation (2-17). 

By drawing a magnitude circle with radius equal to 1.414(4.48+)12.1) 
which is the quantity on the right side of the equation above at W, =e 
and noting the frequencies in the neighborhood of the point on the 1 
curve where the circle intersects with the iB curve, and also noting the 
relative location of the W. (ay, = Bi) plot with respect to the F 
curve, some information can be obtained about the bandwidth. This is 
illustrated with Fig. 2-11. 

In Fig. 2-11, the i curve is drawn for the system (2-17), the work- 
ing point M (12.27, 4.48) is located, and Bit = fit = 12.27 is laid off 
the same way as in Fig. 2-9. The line Boy = (1+1.414)ap, = 10.8 is 
drawn. 

From the intersection point of Bot = 10.8 and Bit * 12.27 lines to 
the left, a distance of 1.414 bit - 17.1 is measured off on the Bor = 
10.8 line and labeled a, = 1.0. The radius of a magnitude circle pass- 
ing through the WD, = 1.0 point on the Bo; = 10.8 line is equal to 
[1.414(4.48 +5 12.1)| - [1.414 NGj1)}. This circle denoted by Ca, crosses 
the Ie curve between the points at w, = 1.0 and Ww et Iauke 

It is easy to predict that Ores 1 is not possible. If it is sus- 
pected that One< } is impossible, it can be checked as follows; 


The check is based on the fact that, for the bandwidth, the magni- 


tude circle must cross the @,(a,, - Boz) plot and 1.414 by, 4%, scale 
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represented on the line Bot = (1+/2 aoe = 10.8 at the same frequency, 
and the fact that if any magnitude circle is drawn and if it crosses the 
O (ayy - Boe) plot and 1.414b, “scale (Bp, = 10.8 line) ae Oe = ce 

and Ww, = We respectively with either Ww <0” or a> on » then for 
existence of a bandwidth frequency, there must be another magnitude circle 
to give a new set of w® and a) with inversed inequality from the previous 


set of w and WW. ‘ 


The magnitude circle Cal gives 


wo, >| and a, =| then Ww, Z w, 
Since Cay is drawn with ne 1 and it crosses the Pe curve between 
W,_ = 1.0 and W, = 1.1 and 1.0<@,<1.1. This is reasoned without 
drawing any portion of the (ayy - By, plot. 
Draw a magnitude circle with W= 0. This circle is of radius equal 
to 2. Age tangent to the line Bot = (1 + {2 aoe = 10.8 and represents 


[2 |ngool 


PD) aor by its radius. Since Gy (a), - By,) plot ap- 
proaches to the point at w,. = l on the ie curve from the point (12.27, 0), 
the circle must cross the WO (ay, = Bit) plot at some frequency greater 
than zero to give an inequality 

@.=0 | WE >0 then OO, eae 
The inequality sign is not reversed from the previous set of Wy. and Ww, s 
and this implies there might be no frequency such that, for 0 < W, <1, 
a magnitude circle crosses the Wilaye - BL) plot and the {2 bit We 
scale at the same frequencies; in other words, w, : w, = Whe may not 
be possible for the frequency range from zero to 1.0. However, there 


may be some frequency between O and 1.0 at which the inequality is re- 


versed since the interval of Oy from 0 to 1.0 is equivalent to a rela- 
tively large interval of original system frequency z= 100 a), 
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Locate the point at a, = .5 on the ae - Bi) plot which is 
half way from the line B,, = aj, = 12.27 to the point at we = .5 on the 
ie curve. This point is located inside the magnitude circle correspond- 
ing to y= 0, therefore, no bandwidth is possible for w. coe 

For the magnitude circles corresponding to .7 < WwW. < 1, the inequal- 
ity a <a. is not changed. This can be checked by comparing the dis- 
tances from the M point to the IF curve and the {2 bat w, scale and 
noting that the (ar, - Bi.) plot lies to the right ‘of the iF curve for 
A) ce a), < POR 

For example, a magnitude circle for i, - ./ will cross the if curve 
at a point where a, is greater than .8 to give Oe 7.8 and WO. < We ; 
Thus for the bandwidth, .7< Wied 1.0 is not possible. 


The only part left to check is the interval .5 < We < -7. Knowing 


that ’ F ¥ 
for &s' .5 wo <Ww and 
for Ws yi a <@ 


it is quite sure that for the interval .5 < we < »7 the same inequality 
will be maintained. If a magnitude circle for We = .6 which is the mid 
point of the interval of the frequency from .5 to .7 is drawn, this cir- 
cle will pass the ls curve at the point where wo. = .6 implying the 
point at W, = .6 on the (aj, - By,) plot will be inside the circle 
in turn 0,<W, which is the same inequality as the cases of Qs o2 and 
W, = -7. Thus the same inequality is maintained for @, from 0 to 1.0 and 
W..<1 is not possible. 

The tests made above are carried out without plotting any portion of 
the ()¢(az_ - Byp) plot except the point at W), = .5 on the W (are Bit) 
plot is located. The tests can be carried out with a sketch of the 
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W ,(az_ - Bip) plot and the sketch can be made simply with reasonable 
accuracy by noting that the horizontal distance from the line Bit = 4y¢ 
= 12.27 to the 0, (are ~ By_) plot is a fraction of that to the me curve 
and the fraction is the value of O,- 

For an interval of frequencies ( W,, Yo); 1f the inequality 

w'<o” (er w'>w”) 
holds at both ends 0, and We then the inequality is maintained through- 
out the interval unless the ln curve or the @ (ar, - By.) plot is oddly 
shaped in the interval, or in other words, unless the rates of change of 
the distance from the M point to the W,(ajt¢ - By;_) plot is not fluctuat- 
ing in the interval. 

In this example, if the WO, (ay, - Bi.) plot is sketched, then one 
can see that the rate of change of distance from the M point to the 
W, (are - By_) plot is steadily increasing except in the region of small 
frequencies. 

Knowing that @. >1.0, the next thing to do is to find another 
magnitude circle to insure an upper bound of Ore as illustrated in Sec- 
tion 2-3. 

Noting that Ca, crosses the I curve at the point where the fre- 
quency is slightly less than 1.1 and also noting that the W (aye bo) 
plot lies to the left of the [ curve for w,.74; it is predicted that 
(Oy, will not be very far from W, = 1.1. Therefore it is natural to 
try a magnitude circle which passes through OO, = 1.1 or WO, =< 1.2 on 
the /2 bye Qt scale. The magnitude circle which passes the point at 
Q), = 1.2 on the line Bop = (1+/2 )ag, is tried and denoted by Cyo. 


The circle Cag represents the right side of Equation (2-19) at Wr = l.25 
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It passes the ie curve at a point where the frequency is less than 1.2 
which implies it will pass a point on the @ (are - By¢) plot where fre- 
quency is less than 1.2, thus Ww, = 1.2 clearly indicates an upper bound 
of bandwidth and 
1.0°< a. < ta2 

To determine Ope precisely, the ,(a,, - By,) plot for 1.0< a, < 
1.2 is needed. 

A magnitude circle Cap which passes through the point at a). = Tee 
on the /2 by, , scale crosses the W,(aj_ - Bye) plot at a), = 1.1 and 


ic 


bandwidth is determined to be 


Ode = 


hose 
This agrees with Fig. 2-10. 

Since Fig. 2-10 or 2.11 is for the system represented by Equation 
(2-17) which is related to the original system by S = 100s, the original 
system bandwidth 0, is 

0), = 100 Ope = 110 
By substituting s = Ore = 1.1 andS = js = 110 into Equation 


(2-17) and (2-16) respectively, the magnitudes calculate to 


= ./07 








ROL sa 


= tO? 





(3) 





=i td 
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2-4. Locus of M points on the By vs By plane for constant bandwidth. 

In Section 2-3, bandwidth of a given system is deceuminad For a 
given system with fixed M point which is determined by the coefficients 
ay and ag» bandwidth is determined by varying © until some frequency 
(or frequencies) (Ws WO, is found to satisfy Ejuation (2-10). 

Suppose, conversely, ® in Equation (2-10) is fixed at bandwidth 
@,, then there must be some a, and ag to satisfy Equation (2-10), and 
the pair (a, » ag) defines the M point on the Bg vs By plane for a system 
to have bandwidth equal to We If all of such pairs (a; » 49) satisfying 
Equation (2-10) are plotted on the Bo vs By plane, then a curve which 
represents the locus of M points on the Bp vs By plane for fixed band-~- 
width @, is obtained. This curve is called the constant bandwidth curve 
for bandwidth equal to Wy or merely constant w,, curve. 

If a constant Oo, curve is constructed on the Bo vs By plane, and if 
the M point is chosen at any point on the curve, then the resulting sys- 
tem will have bandwidth equal to Ws: 

Since the synthesis of a feedback control system by Mitrovic's 
method is to choose a suitable M point to guarantee the required system 
performance specifications, the constant @, curve will provide informa- 
tion in choosing the M point to meet specified bandwidth requirement. 

By fixing = @, in Equation (2-10), the constant Ww. curve is 
defined by 

|(a,-P, (@x)) +). (a, -E,(@) | = 2 | NG | (2-20) 
where By (Wy) and By (W,) are the coordinates of the 
points on ie curve at @)= W,. 


On the left side of Equation (2-20), By(@,) and Bol @,) are fixed and 


those values can be read off on the lis curve. 
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On the right side of Equation (2-20), noting that N(S) is the numer- 
ator polynomial of the closed loop transfer function, [NC jp) contains 
a which is treated as a variable on the constant Op curve. Moreover, 
if N(S) is of order other than zero, it may contain by which is the co- 
efficient of the s* term in N(S), and by is not independent of ay which 
is the other variable on the constant Ds, curve. Thus, even though 
Equation (2-20) is a perfectly general expression for defining constant 
Gy, curves, it does not clearly indicate the geometrical shape of the 
constant Oy, curve on the By vs By plane. 

In the following, the cases of zero order, first order and second 
order of N(S) are considered. 

N(S) of zero order. 
When a closed loop transfer function has no zeros, N(S) = ag and 


Equation (2-20) becomes 
| ( 2o~B,(p)) +), (4-8, (OS) |= By 


by removing the absolute sign, and manipulating 


(4+ Byw,)) - w, (a,- By(a)) = z Bm.) (2-21) 


Since Oy» Bo( @))> B,C @) are fixed, Equation (2-21) indicates a hyper- 
bola with following geometrical properties on the By vs By plane: 

the center of symmetry at (By, , -B,(@,)) 

the vertices at (Bas), C14/2)B, (We) 

major axis is B, = B,(,) 
The two branches of hyperbola (2-21) lie symmetrically with respect to 


the Bo = “BoC W,) line. The upper branch lies above the B, axis with 


minimum point at (B,(W), .414 Bo(Wp)) which is one of the vertices. 
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For a stable system, only the portion of constant Wp curve which lies 
in the stable region is of interest. 
It is interesting to note that the locus of the vertex of the branch 


Ee 


of the hyperbola which lie in the first quadrant is similar to the : 


curve, and the locus of the center of the hyperbola is the mirror image 
of the ie curve, 

Those can be seen as follows. As shown previously the vertex of the 
upper branch of the hyperbola is at (8,(@))> -414B)(@,)) and the center 
point is at (B,C); ~By( ©,))- 

Considering Wy as a continuous positive variable, then the locus 
of the vertex point is defined by 


xX 


B, (9) 


Y= .414 B,(@) 


and in like manner the locus of the center point is defined by 


X By(W) 


Y ~By(@) 


where X and Y axes are identical to B, and Bg axes respectively. 
The two sets of equations can be compared with the E curve defined by 
By = B,(@) 
By = Bo( ®) 
Consider the consistency of principles involved in constant Ww, curve 
and in determining bandwidth described in Section 2-3. 
In Fig. 2-12, suppose a constant @,, curve is plotted from Equation 
(2-21) assuming the ’ curve as given. The point P is at @=: Wy on 


(o) 


the K curve. The constant @,, curve is denoted by C & The vertical 


b° 
axis denoted by Y is on the By = By (Hy) line which passes through the 
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Fig. 2-12 Showing consistency of Section 2-3 and 


Section 2-4 about bandwidth. 
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point P and V is the vertex. The horizontal axis and the lower branch of 
the hyperbola which lie below the By axis are not shown. 
The portion of the constant @),, curve which is of interest for a 


stable system is the portion of C under the ie curve. 


Wb 


Suppose an M point is chosen arbitrarily on C p 38 shown. Then a2 


Ww 
Ww) (a; ~ Bi) plot results as described in the previous sections and a 
magnitude circle Ca with radius equal to (2 a, must pass through the 
point Q at which Ws: a, on the W (ay - By) plot to have a consistency 
for the principles involved in the constant Oy, curve and determining 
bandwidth as described in Section 2-3. This is to show that 
MQ =/2 a, 
Locating a point R at the intersection point of the By = ay line and 


V2. 
ee Yn 5 ea < 
MQ=(R2+ ee) =((B(w,0- 4.) +(e (2,-B (ony 


But (a ag) satisfies Equation (2-21) and it rearranges to 


1 ? 
Es . d Z 2 
(4,-B,(@,)) + Wp (4,-B,(,)) = 2c) 


and MQ =/z an 


N(S) of first order. 
When a closed loop system has one zero 
N(S) =b,S+ dy 
For unity feedback control system, b; enters into a, and a, and by are 
not independent of each other. 
Suppose an open loop transfer function of a unity feedback control 


system as 


; Ee ays 
GEia EetrennS ae 


48 





where bis ay and ay are variable and dys are fixed constants. 


Then R eo : as i. aa Bae ne a, 
G tds b,S +4, 
aaa =a 5 wyakl_ = a ee = 

an pole", a oF ae ent? 
where doy & di-l ee a ay = dy and ay =dj+b, 


N(S) can be written in the form of 


With this N(S), Equation (2-20) is manipulated to 
2 2 2 
(4,+B,(o,)] + Wy (a, +(B,(, -2d,)] =D (sp, d,) (2-23) 


where D(®s, di) = 2[ Be(Wp) + OF (B,(@,)-d.'] + 470, 


Equation (2-23) defines the constant om curve for any system whose closed 
loop transfer function has one zero. 
Noting that for D(W,,.d,) 70, Equation (2-23) is an ellipse centered 


at (2d, - By(@4) “By (,)) with semi-axes of lengths equal to 


DW and /D(,,d;) . (Horizontal and vertical axes respectively. 
b 


For each dq, » as Cy varies there results a family of ellipses. 

It is seen in Section 2-3 that for a system represented by Equation 
(2-22), larger by results in larger bandwidth. 

For any unity feedback system assuming all positive coefficients in 
its open loop transfer function, b, is always less than ay and its upper 
limit is a, which is the case of a type 2 system, whereas its lower limit 
is zero which is the case of a system with no zeros. 


For the upper Limit of by = a) corresponds to d, = 0, and letting 
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dy = O in Equation (2-23), the constant p curve is then 


(Bot By(Ws)) + Op (A, +B (opr) = 2 (Bp) +p Be) (2-24) 
Equation (2-24) can be applied for a type 2 system with one zero in its 
closed loop transfer function. For any system with bi< a, or dj>0, 
if the constant Wy, curve is plotted by Equation (2-24), and if the M 
point is located on the curve, then the system will have bandwidth less 
than Oy» thus Equation (2-24) defines a locus of M points on the 
Bo vs By plane for some upper bound of bandwidth for a system whose 
closed loop transfer function is represented in the form of Equation 
(2-22). 

This will be illustrated with a specific system 


35 +a 
et ae ee Se ee 2-25 
R(S)= cal gFy $25 aS +a ee 
with associated constant ©, curves as shown in Fig. 2-13. 

For the system (2-25) the i curve is defined by 


@ 


2 
By = ob) 
> A 
Bo s@M-W 
and this curve is plotted in Fig. 2-13 as shown. 

If the system had no zero, then constant G, curves would be defined 
by Equation (2-21) with By and Bo defined as above. For this case, con- 
stant ®, curves are constructed for a. = .l to @} = .5 as shown in 
Fig. 2-13. If M point is chosen at any point on C 5 which is a constant 


@,, curve for ,, = .5 then the system 


a 
2 (2-26) 


Co wee Oe eS 
NS e+ S° +S as ee 


will have bandwidth equal to .5. 
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The dotted curve which is a parabola is the locus of the vertices 
of the constant 4), curves for the system (2-26). 
A constant Ws, curve of wy, = 1.0 for a system 


(Ce Stone 
R(S)= SH 4574574 4,8 + a, (2-27) 


which is obtained by letting by = a, in Equation (2-25) is plotted by 


1 


Equation (2-24) and denoted by C in Fig. 2-13. If the axes scales 


1.0 
are made the same, the Cio would be an are of a circle, since, for 
W,, = 1, Bo(1) = 0, B, (1) : 1 then Equation (2-24) is 
a* +(Aa,+1) 22 
which is a circle centered at (-1, 0) with radius equal to /2. 

If the M point is chosen at any point on Cio then the system 
(2-27) will have bandwidth equal to 1.0. 

Suppose the point M(.407, .089) is chosen for both systems (2-26) 
and (2-27). Since the point M (.407, .089) is on both C , and Ci o> 
the systems (2-26) and (2-27) have bandwidth @, = .5 and @ = 1.0 re- 
spectively. 

For system (2-26) 

OFF ] 


ee! =.70 
Co, CE 029 Ignis ia 


on 


For system (2-27) 


E | 4078 + .089 — 
eho De SF 4 4075 Poe op ey 
=i I. 


Then the point M (.407, .089) for the system (2-25) guarantees that 
0.5 5 W, < tao 


for 0 < by < ay 
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N(S) of second order and up. 

When a closed loop system has 2 or more zeros 

NSN, 4 DSS pees 

If the coefficients bj, S for k?72 are fixed (in order for a sys- 
tem to be analyzed or synthesized by Mitrovic's method on the Bg vs By 
plane, the coefficients aS for k22 must be determined and likewise 
for byS ), then by proceeding in the same way as the case of N(S) = b,S 
+ @, the equation defining the constant as curve, which is similar to 
Equation (2-23) with due modification to account for the coefficients bo, 
b35 - - will be obtained and likewise for the equation defining constant 
@,, curve of certain upper bound for bandwidth. 

It must be noted that for a high-velocity constant system by is very 


close to ay since the velocity constant K, is related to a,, ag and b, by 


_ Saco 
Ky a 4,-b 


Therefore, for a high K, system, the constant (), curve obtained by let- 
ting by = aj will be a good approximation to the actual bandwidth. 
Also in lag compensation, aj, = by and as shown in references 1 
and 2, lag compensation problems are solved by letting a, = by >» conse- 
quently, the constant Wh curve obtained from a, = b, will be good in 


lag compensation problems. 
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2-5. Bandwidth of third order systems. 

It is shown in references 1 and 2 that the equations of third order 
systems can be normalized to give a universal Mitrovic's chart for third 
order systems, and this chart is constructed in references 1 and 2. 

In the normalized chart, i is a straight line extending from the 
origin to infinity with slope of unity on the Bg, vs Bin plane, and the 
frequency Den in the normalized system is related to the frequency ® in 
the unnormalized system by 

W= RO 
where & is a constant, the normalizing factor, and sub- 
script n denote normalized quantity. 

Assume a third ofder system with no zeros in its closed loop trans- 
fer function: 


jee ao 
g(S)= Bo alSs ace 


It normalizes to, by transformation S = a9s 


Cp. Bon i 
pO 7? ug2 + a ndt Gn Cae 


where 4,,, = 2, Tee 


3 (2-29) 
a = ayes 


Equation (2-29) defines the relation between coefficients of the 
normalized and unnormalized systems, consequently the Mitrovic's variables 


are related by Equation (2-29) between normalized and unnormalized sys- 


tems, a 
namely By, =B,/ az 

: (2-30) 
By, = B,/ 4, 
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Applying Equation (2-21) to the system represented by Equation 


(2-28), the constant OO, curve is defined by 


(2. +B,(@un)) + Oana (Wa = 2 By an (2-31) 


where Won is the bandwidth of normalized system and is re- 


lated to the unnormalized bandwidth on by @ 


b * ora: 
Bin and By, are Mitrovic's variables for the /f curve in the 
normalized third order chart, and are related with unnormalized 
variables By and By by Equation (2-30). 

Equation (2-31) defines a family of constant Wn curves on the 
Bon VS By, Plane if a set of ) is supplied. 

If this family of constant Onn curves are superposed on a third 
order chart, then it will provide additional information about the third 
order system. 

Fig. 2-14 shows the family of constant bandwidth curves for third 
order systems whose closed loop transfer function is represented by Equa- 
tion (2-28) on the normalized By vs By, plane. Note that the portions of 
constant bandwidth curves defined by Equation (2-31) which lie in the 
unstable region are omitted. 

Fig. 2-14, the constant bandwidth chart for third order systems, can 
be superposed on Fig. 2-15 which is the third order Mitrovic's chart con- 
structed in references 1 and 2 to provide additional information concern- 
ing system bandwidth in both analysis and synthesis of third order sys- 
tems, 

To determine system bandwidth is to locate the M point of given 


system on the constant bandwidth chart and consult with the constant 
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bandwidth curves. Also it is an easy matter to choose an M point such 
that it guarantees a given bandwidth specification. 

From Fig. 2-14, the effect of moving the M point (im the stable 
region) on system bandwidth can be analyzed as follows: 

Suppose a point M (apap) is located in the stable region. Fer 


fixed a increasing a 


On 


easily seen from Fig. 2-14. This agrees with common experience. 


in? results in increasing the bandwidth as is 
For fixed a ~? the effect of changing ai, 0m bandwidth is dependent 
upon the value of Aone This can be analyzed by noting the slepe of con- 
stant Onn curves. 
It is seen in Section 2-4 that the lecus of the vertex peint of the 
family of constant bandwidth curves, the hyperbolas, is 
X = Bin(Wen) 
Y = 0.414 By, (Wen) 
where X and Y axes are identical to the Bin and Bo, axes réspectively 
and Binf tn) and Bon 6 O? are the equations which define the . curve. 
Noting that the iB curve for the normalized system is a straight line 
emanating from the origin with unity slope in the first quadrant, then 
the locus of the vertex points of the constant Wn curves is a straight 
line with slope of .414 defined by 
Bon = o414B,, (2=32) 
on the Bon VS Bin plane. 
Since a vertex point of a hyperbola is the point at which the slope 
of the hyperbola changes its sign, every constant Wn curve Lying in 
the stable region is such that the portion of the curve which is above 


(to the left of) the line defined by Equation (2-32) has negative slope 


and dually. 
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Dividing the stable region into two regions by the line defined by 


Equation (2-32), it is easy to see that if a, is fixed, then by increas- 


On 
ing aw the bandwidth is increased as far as the M point is in the 
region above the line and dually for the region below the line. 

If a closed loop system has a zero in its closed loop transfer 


function, then the normalized system equation is 


dnd + Aon 


Cy ee Se eh 2 
moe Si tBrt+ And + Ayn pee 


where bin is normalized by » the coefficient of st term in the 


numerator polynomial of the unnormalized system, and the rela~ 
is 


tion between by and bia 


bin = b faz 
The arguments developed in Section 2-4 for a system whose closed loop 
transfer function has one zero applies to the third order system repre- 
sented by Equation (2-33). 
By letting b,, = 4;, in Equation (2-33), then the resulting constant 
bandwidth curves which are defined by Equation (2-24) in Section 2-4 


(with due account of normalized third order B ) are: 


( Ayit Ba (Oynl) Den (Amn By (pd) = 2[ Benn) + Opn Bn (| (2-34) 


where Bin and B a are Mitrovic's variables for the i curve 


0 
on the normalized third order charts. 

The portions of the ellipses lying in the stable region are plotted 

in Fig. 2-16 which can also be superposed on Fig. 2-15 to provide infor- 


mation about bandwidth for any third order system with one zero in its 


closed loop system transfer function. 
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As illustrated in Section 2-4, Fig. 2-16 constitutes a bandwidth 
chart for a system represented by 


ain + aon 
A+ 874+ Ayn Bt Aan 


a,S ae a, 


yea ee 
= Ceca oe, or a) 


3) 


and a bandwidth chart for the upper bound of bandwidth for a system rep- 


resented by 


pe bS + 4, 
<(s)= Se a ea, 


or equivalently for the normalized system in Equation (2-33). 

The effect of moving the M point in the stable region on bandwidth 
can be analyzed in Fig. 2-16 analogously as in Fig. 2-14. 

Note that the ellipse defined by Equation (2-34) has a center point 


moving along the line Bp, = B,, in the third quadrant as Orn varies, 


n 
and also note that the ellipse has negative slope in the first quadrant. 


Thus it is easy to see that by increasing either ay, oF a the bandwidth 


On? 
is increased. 

In Fig. 2-17, 2-18, 2-20 and 2-21, the constant bandwidth charts, 
Fig. 2-14 and 2-16 are plotted on expanded scales. The third order chart 
Fig. 2-15 is correspondingly expanded in Fig. 2-19 and 2-22 which go with 
Fig. 2-17, 2-18 and 2-20, 2-21 respectively. 

The following examples illustrate determining bandwidths of third 
order closed loop systems by using the bandwidths charts developed in 


this section, and also illustrate a design of system to meet exact band- 


width by using the third order system bandwidth chart. 
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Example 2.4 
For the system in Example 2.1 


one — _ 
gE)= Commas 


the bandwidth is found to be wy = 1.76. Now one can find the bandwidth 
on the third order bandwidth chart. 


The given system normalizes to, by using the transformation § = 3s, 


R hig > 222 oe 


and the normalized bandwidth Won is related with the unnormalized band- 
width ae by a = 30, Since the system has no zeros, one enters into 
the chart in Fig. 2-14 or the expanded chart in Fig. 2-17 or Fig. 2-20. 
The M point of the normalized system seems the chart in Fig. 2-17 is 
suitable. By locating the point M(.222, .148) on the chart in Fig. 2-17, 
Wphn is found to be .59 by interpolating GQ, = +55 and W,. = .6 curves. 
Then the bandwidth of the original system is 
Oe Spy =(4y59)—/-77 

which agrees with the result in Example 2.1. 

For the system in Example 2.2 


Cg) = 2Z21DDD SK + 44-7000 
ROE Go OCIS “+ (226775 eee 


the bandwidth is found to be Q, - 110 in Example 2.3. 

This system is third order with one zero in its closed loop transfer 
function. If the bandwidth is determined from the third order bandwidth 
In? bin then the bandwidth so determined must be not smaller 


than 110. Noting that a, = 122671 and b, = 121000 are not greatly dif- 


chart for a 


ferent, the bandwidth determined from the chart for bi = 4, will be 
very close to the actual bandwidth 110. 
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The system normalizes to, through the transformation § = 1671.1s, 


(0 Le 
Bet S7 +. 0498S +, 20096 


Entering into the bandwidth chart in Fig. 2-21 with the normalized 
M(.0438, .00096) point, the normalized bandwidth is found to be 
Opn =.07 
by interpolating between the two curves for 4, = -05 and Bate +O7>. 
Then the unnormalized bandwidth is 
ee asa =e 77. Coma 
which is very close to the actual bandwidth W,, = 110 and agrees with 


the arguments developed in Sections 2-3, 2-4 and 2-6. 
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Example 2.5 


Given a system block diagram below 









oe 
S(Seresre 2) 





adjust K,K, and K, with the following performance specifications: 


imme 53 
2) Ky = 10 
3) Wy, zs 20 


For the system without feeding back derivative signals K must be 
set to 90 in order to meet the Ky, specification, and the system is at 
stability limit which can be seen from the third order chart, 

The closed loop transfer function with the derivative signals feeded 


back is 


ee 
RE IS S?+(l0tk Ka) S7*+(9 +KKt)S + K 


and normalizes to, by using the transformation § = a38 = (10+ KK,)s, 


ee 
b7 + a + ain S + Ion 





(SC = 
E6)= 
where; ain =(94+KKe)/U 0+ K Kay 


ed Uae 
and a3 = LO+ KK, = normalizing factor. 
The velocity constant and bandwidth of the unnormalized and normal- 


ized systems are related by 
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Ry =(10+ K Ka) Ky, 

WOp=60 + iC Kg) Wen 

Ky/ @y = Kyn / Opn = IZ 

where the subscript 'n" denotes the normalized system quantity. 
Thus the ratio of the velocity constant K,, and the bandwidth Opn of 
the normalized system must be maintained at 1/2. 

In Fig. 2-22a the FE and [3 curves are duplicated from the third 
order chart (Fig. 2-15) and then some constant bandwidth curves are 
superposed from the third order system bandwidth chart (Fig. 2-14). 

In order for the system to have ‘%= .3 the M point must be located 
below the [i3 curve in the first quadrant on the Bo, vs By, plane. For 
the normalized system K,, = &9,/a,, indicates that for an M point located 
anywhere on the normalized Bon VS Bin plane if a line is drawn from the 
origin to the M point, then the slope of the line is equal to the normal- 
ized velocity constant K,,. Therefore for this problem the M point must 
be at the intersection point of a constant Wy, = Opn curve and a line 
drawn from the origin such that the slope of the line is (1/2) Dn and 
moreover that M point must be such that it insures 2.3. For any 
chosen constant Orn = Oye curve the corresponding straight line with 
slope 1/2 Wh can be drawn and then the line almost always intersects 
with the constant 4G), = Onn curve. then intersection occurs, location 
of the M point at the intersection guarantees that both Ky, bandwidth 
specifications have been satisfied, but may not guarantee the % specifi- 
cation. In order that all three specifications be satisfied, the inter- 


section must lie in the area enclosed by the .3 curve and the Bin axis 


in the first quadrant. 
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Fig. 2-22a Solution of Example 2.9 using 
the 3rd order system band width chart. 


T2 








If one chooses the constant bandwidth curve for On = 0.7 and 
draws the line with a slope of 0.35 from the origin as shown in Fig. 
2-22a, then the intersection point is at (.575, .204) and this point in- 
sures 2 .3 which is evident from the [3 curve configuration, thus 
one chooses the M point at M(aj,,, agn) = MC.575, .204) to give 


= 9 + K Ke = 
“In = Got Kha) 575 


K 


= ——_____._._ =, 244 
“On G0 + KRa)™ 


Since 6), = .7 is chosen the normalizing factor a3 = 10+KK, is de- 
termined from the relation 

O, =(0+K Ke) Opn =10 
from which 1O+tKK, = 28.6 is obtained. Substituting this value of the 
normalizing factor into the denominators of the two equations for the 
coordinate of the selected M point above, K,Ke and then K, are deter- 


mined and the solutions are 


kK - 4675 
Ry - 0.0987 
K, = 0.00388 
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2-6. Frequency response of the fourth order systems. 
The iE curve of a fourth order system can be normalized so that a 
single curve can represent the Ne curve of any fourth order system. 
Suppose a fourth order characteristic equation is; 
F(S) = S“4 a,S-+ 4,6 + 4,5 eno 
then the |; curve is defined by 


2 (2-34a) 
B= aww” 
Equation (2-34a) represents a concave down parabola whose axis of sym- 


metry is B, = aoa;/2 with vertex at (aa3/2, aga3/4), Bg = O occurs at 


W@W =0Oand Ws: ag : 
% 


Assume a normalized frequency on such that Ws a5 Ven then 

Equation (2-34a) becomes 

B,= 42 as Wen (2-35) 

Bo= aa (@en Sz Den ) 

2 ‘ 
Let B, = By/aza, and Bon = Bo/ a5 then Equation (2-35) can be writ- 
ten as 
Zz 
Bin = Wen (2-36) 


Bon a Den — cee 
Equation (2-36) is exactly the same equation for the iB curve as is asso- 


ciated with a characteristic equation whose coefficients of the s+, 3? 


and $2 


terms are unity. 
If a curve defined by Equation (2-36) is constructed, then with 
proper scaling factors of coordinates and frequencies, it will represent 


exactly the same curve represented by Equation (2-34a) and thus ie 


curves of fourth order systems are normalized. 
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+ 
The normalizing property of the iF curve for fourth order systems 
can be used to obtain the frequency response of any fourth order system. 


Suppose a fourth order system closed loop transfer function 


Sa ies a 
£()=- | BS 0 ee (2-37) 


§“+ 4,534 4,8°+4,8 + a, 
which has a second order polynomial in the numerator. For the case of a 
fourth order system with third order numerator, add b3S to the numerator 


of Equation (2-37). By using a transformation or frequency scaling 
co a~ 4 


Equation (2-37) can be written as 


C (s)= ede s apa + dy 
-{ “3. “2 
or E a) b, a, b+ B42 4 + doto (2-38) 
R + ants a+ 648 ae ee 


The systems represented by Equation (2-37) and (2-38) are equivalent 
Since the transformation § = ais is one to one. Note in Equation (2-38) 
that the coefficients of the S* and s4 terms in the denominator are made 
unity, and also, the coefficients in the numerator and denominator of 


Equation (2-38), by comparing those of same power terms, are scaled by 


the same amount. For instance, bj and a) both are multiplied by ore 


> 
and likewise for ag and bo , etc. 


Let agn = 20% 45 
ain 2, /a2dsz 
99n = 4,/a2 =! 
Q3n = 4;/d,=1 (2-39) 
bin b, / a, a3 
ban = be / a2 


% 
9or =a, / as 


1 
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then Equation (2-38) is 





big Pa eee. 
£6)= ee LM (2-38a) 


wae ee a B+ eet a ieee ase 


Z 
By substituting s =J Win? and noting that Bin = Wen? Bon 2 on 2 Bre 
which are the equations defining the normalized Ie curve for fourth order 


systems, the magnitude of the frequency response is given by 





| Cre ) ay A yy Ven | 
| (Bon Bon) ny ang, (4, -Byn) | 


Fem) (2-40) 





thus the frequency response of fourth order systems is represented in 
terms of normalized quantities. 

The methods described in Sections 2-1, 2-2, 2-3 and 2-4 can be ap- 
plied to get the frequency response of the fourth order system. The 
frequency response of the original system represented by Equation (2-37) 
can be obtained from 


(Ap—bw ) +) bw ' 
@,-8,) +j@(a,-B8) i 





Le) a 


which results from applying Equation (2-6) to the system (2-37). The 
advantage of using normalized Equation (2-40) over the original Equation 
(2-41) is that, since the normalized hs curve for any fourth order sys- 
tem is a fixed curve, it is not required to plot the i, curve for each 
system. 

Knowing how to interpret Equation (2-41) on the By vs By plane which 


is associated with the original system, Equation (2-40) can be inter- 


preted on the normalized By, vs By, plane without difficulty. 
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The quantities or coordinates in the normalized and unnormalized 
systems or Mitrovic's plane are related by Equation (2-39) plus the fre- 
quency relation 

Zz 
(ae) = a, Hen 
The M point is located on the normalized Bon VS By, Plane by the normal- 
ized coordinate (41, . @on)° 

The imaginary term in the denominator of Equation (2-40) can be rep- 

resented in the same way as for Eyuation (2-41),namely by (a, - By) 


plot; however, it must be multiplied by ay to get the a Wen (iy Bin? 


t 
plot, and likewise for the imaginary term in the numerator of Equation 
(2-41) if a system of interest has by other than zero. 

In Fig. 2-23, the normalized B curve for fourth order systems is 
plotted on the Bo, vs B,,, plane with the same scales for both axes. Any 
point M, (41 1° 4 on? located on the Bo, vs By, plane is equivalent to the 
point M(a) 2) on the Bo vs By plane -where ain? 41 and Qon? a are re- 
lated by Equation (2-39). Any frequency Wey on the normalized i" curve 
is equal to &® divided by as! where @® is the frequency on the 
unnormalLlized ig curve at the point which is equivalent to the Win point 


on the normalized ee curve. 


The following table is a comparison of the pe curves for normalized 


and unnormalized systems; 


Normalized [is curve Unnormalized ie curve 


Axis of symmetry Bin = 0.5 By = «5 8783 
Coordinates of vertex (55-222) (.5 aga3z, .25 ad ) 
Frequency at vertex Den = +707 Q= .707 fag 
Freq. at By intersect (fe = (0) Den = 1 @=0, @Me faz 
Cc dinates of B 

oe ee (0, 0), (1, 0) = (0, 0), (0, aga3) 


ae 
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The equation defining constant bandwidth curves of fourth order sys- 
tems on the normalized By vs By plane contains a parameter ane» therefore 
it is not a practical way to try to provide a universal chart of constant 
bandwidth curves of fourth order systems, because for each a9; there re- 
sults a family of constant bandwidth curves. However, it would be worth 
while to investigate the effect of varying a>, on bandwidth, since if 
the effect of Ao, on bandwidth is known, then if a constant bandwidth 
curve is drawn or bandwidth is determined for some a9, then it may be 
possible to get some information about bandwidth for another a>;. 

Assume two fourth order systems such that when they are normalized 
they result in identical system equations except they have different 
ao'Ss name ly ash = k, for one of the two systems and af, =: k2 for the 
other where kj and kz are constants, k,;Z ko. 

The two systems have a common M point on the normalized Bg vs By 
plane. Suppose this M point is located as shown in Fig. 2-23. Since 
the normalized if curve and M point are common to both systems, the 
Wen (41, > Bin) plot will also be common, but each system will have a dif- 


ferent a (a1, - By,) plot as shown in Fig. 2-23. (k; = 1 and ky = 2 


are assumed for a simple illustration.) 





Since the bandwidth of a system is obtained at the point where some 


0). 
magnitude circle intersects the = (a, 7 B,,) plet, and noting that 
ot 


the smaller Ap, Causes the kein - Bin) plot to be displaced farther 
from the Bin = ay line than the larger Boy does, any circle centered at 


the M point intersects the (ain * Bin) plot at lower frequency for 
Ze 


the smaller aie: For example in Fig. 2-23, if the system has no zeros 
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in its closed loop transfer function, then the magnitude circle Cap to 
determine bandwidth has radius equal to (2 a, = .0707 and it intersects 


with the two en By.) plots at P, and Py to give the bandwidths 
lat 


Oo = 3342 for fag, = 1 
Wyn = «35 for fae =: 2 as shown. 


It is obvious the above reasoning holds for systems whose closed 
loop transfer function has any number of zeros, thus a conclusion is ob- 
tained as follows: 


The effect of a,, on bandwidth of a normalized fourth order system 


2t 
is to increase the normalized bandwidth with increasing aot: 
It must be remembered that the normalized bandwidth Win is related with 


unnormalized bandwidth OO, by 


Op = (2, Opn 
As a simple application of the above conclusion, refer to Fig. 2-13. In 


Fig. 2-13, the hyperbolas, the constant bandwidth curves are for a system 


@ © 
Mel ce? 10>, S ae 


This system can be regarded as some system which normalizes so that the 


= 1. In other words, if for 


coefficient of the s° term is unity with ao, = 


a normalized system 





C (2) Bom 
Ro" 2 ey a oI ae, 
a a s ap fast a+ 


if constant Onn curves are plotted with any =z 1, then there will result 


exactly the same configuration of curves as in Fig. 2-13. For the 
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M(a,, »49,) point as indicated on the curve C5 5 e is known to be .5. 


b 
If there is another normalized system with the same M point but with aoe 
slightly different from 1, then Onn of this system will not be far off 
from .5 in the well-known direction of deviation from .5. Same reasoning 
can be made if the M point is located anywhere on the Bon VS Bin plane. 


The following numerical examples are to illustrate using the normal- 


ized ie curve in obtaining frequency response of fourth order systems. 
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Example 2.6 
Given a closed loop system transfer function 


1600 


SC GN ee 
ge (S) St + 22874172574 6908 +1600 


(2-42) 


obtain the frequency response (magnitude). 
The problem can be worked out without normalizing the systems; how- 
ever, it will be shown on the normalized Bo, vs By), plane. 
For the given system, the coefficients are 
ag = 1600 
a, = 680 


172 


e2 
a3 = 22 
By using relations between normalized and unnormalized quantities in 
Equation (2-39), the coefficients are normalized to 
2. 
Aon = a7 a =.0542 
and Van ena, = 57> 


The frequencies of the normalized and unnormalized system are related 


by ) =/ a2 Oey =/3. /Wen 


From Equation (2-40), the frequency response (magnitude) of the given 


system is, in terms of normalized quantities, 






05F¢2 


Ly: = 
EG | (0542 — Dyn) +) SS (.18-Bin) 





In Fig. 2-24 the point M(a, > a9) = M (.18, .0542) is located and 


the lines By, = 41, = -18, Bon = 20n = 00542, Bon = a+/2 4on = .13t 
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are drawn on the normalized Bo vs By plane. 

Since the line B, = (1 + 2)ay = .131 crosses the Re curve at the 
point where the frequency is slightly less than On =z 4, and since the 
system has no zeros, the bandwidth Ww. must be less than .4. From the 
ie curve and M point configuration on the normalized Bg vs By plane as 
shown in Fig. 2-24, it is seen that the first one-half part of the 
curve is sufficient to see the frequency response of the system. 

In Fig. 2-25, the ie curve, M point and the Bin = a Bon = ?0n 
and Bon = (1 +f2 aon lines are drawn to expanded coordinate scales. 
Frequency response is represented in the same way as in Fig. 2-8, the 


only different thing from Fig. 2-8 is the Fe (ain - Bin) plot. 
oy 


a) 
an ; 
As illustrated previously, the e (43, - Bi»? plot is con- 


structed such that the horizontal distance from the line Bin = 41, = - 18 


Q), 
to any point on the — (an - Bi plot is equal to the horizontal 


I oe 


distance from the line Bjy = ajp = -18 to the e curve multiplied by 


Hen S Men The Yen 


Va 9 595 (ae 


iy if the ordinates are the same as in the O(a, - B,) plot. 





(41, ~ B,,) plot bears the same frequency as 


Some magnitude circles are drawn to show the magnitude A =| <gean) 
Dyn 
at some frequencies. The magnitude circle tangent to the re (4, Bi? 
V doe 
P +n, 
plot which is for resonant point touches to the ipa (ayn- BL.) plot at 


Wtn = -24 with radius equal to .05 to give 


the resonance frequency W,, 024 or W, = (13.1), = 3.14 


the resonance peak Mpw = -0542/ .05 = 1.08 


The magnitude circle with radius equal to fo @g = -/65 crosses the 


hen 


——— (a), - By,,) plot at @,, = 388 to give the bandwidth 
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Wyn = +388 or =O, = (13.1) (.388) = 5.09 
Substituting S = (0 in Equation (2-42), the magnitude for @= (1) = 


3.14 and Q= @, = 5.09 calculate to 


a 
|=garl= 1.1 for resonance peak re 


and | ease .71 for bandwidth 


Magnitude circles are drawn in Fig. 2-25 for several frequencies and mag- 


nitudes are indicates, thus the frequency response is obtained. 
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Example 2.7 
Given an open loop transfer function for a unity feedback control 


system as 


_ 3/.6 (S+i)(S +5) 
GS)= SCS +.3323)(S2 yer) 


determine the bandwidth of the closed loop system. 
The closed loop system transfer function is 


- WES +1908 +158 z.. 
ROS) = “SF 43.3338 455.6S4/ 96.665 FISD (2-43) | 
which has a second order numerator. 


The quantities to be normalized to work on the normalized By vs By 


plane are, using the relations in Equation (2-39), 


aon = 158/55.6 = .0512 
a1, = 196.66/(12.333) (55.6) = .286 
bin = 190/(12.333) (55.6) = .277 
bon = 31-6/55.6 = .568 
J 49¢ av 55.6/12. 3332 = -603 
{42 - 7.44 

Q) = 7.44 Win 


The frequency response of the system in terms of normalized quantities 


is defined by, using Equation (2-40), 


(-0512--568 Men ) ty e277) 


(0512 - ise) ty O22 (.296 - Bn) 





Ao 
es Wm) 
For the bandwidth 


(os12- Bn) +] 22 4 (.286-Byn) = (2 \( (.0512- 568 Din ) +; SE (277) (2-44) 
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The left side of Equation (2-44) can be represented by the usual 


Ee (.286 - By) plot. On the right side of Equation (2-44), both real 


and imaginary terms are frequency dependent. As illustrated in Section 
2-3, the real term with the factor /2 is represented on the Bin = 4in = 
-286 line and the imaginary term with the factor (Bl is represented on 
the Bo, = agn = -0512 line. Those are shown in Fig. 2-26. 


In Fig. 2-26, the normalized iF curve with the point M(.286, .0512) 
2 


is shown. On Bj, = .286 line, Vo bo, Wen = -805 ee is calibrated as 
follows. The point corresponding to Win = QO which is located below Bin 
axis is at the distance of /2 a = .0723 from the M point. From this 
point up, the B,, = .286 line is labeled with Den such that the distance 


from the Wn = 0 point toa We S Oey point is equal to 


ery Oe = «805 Oe so that at Orn = Dei » the distance from the 


M point to the @,, = Wr, point on the By, = .286 line is equal to 
iS: z 
{2 (a= Bien wes ) =|.0723-.205 dy, 


For example at Oren Ce Ges {2 (aon - bn 2) is equal to the dis- 
tance from the M point to the Win = 9 point on the Bin -286 Line, 
which is equal to .73. 

The imaginary term in the right side of Equation (2-44) is repre- 
sented on the Bon = #0n = -0512 line as follows. On the By, = ao, = 20512 
line, each point to the right of the M point bears frequency Wen such 
that a point corresponding to @,, = @,, OF the Bo, = .0512 line is 


located at a distance of /2 bin Dei / [ay = 05 Des : For example, 


at Den = .9 on the Bon = .0512 line, the distance from the M point is 
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(.65)(.9) = .585, thus the line Bon = .0512 has linearly displaced Den 
whereas the line Bin = -286 has quadratic Den 
With those two lines, the quantity on the right side of Equation 


(2-44) can be read off at any frequency. For example, at Wen = 79 etne 
2 » Wen : 
quantity {2 |¢.512- .568W,") +) 73 (.277)| is the length of the diagonal 


of the rectangle formed by the point M, Wea = .9 points on Bin = 286 
line and Bp, = .0512 line, which is the radius of the magnituce circle 
OF 

Then the bandwidth is determined in the same way as in Example 2.3. 


In this example, as seen in Fig. 2-26, 


a) Q of 


bn 
Noting that frequencies are scaled by (z= 7.44 Orn ,» bandwidth of the 


original system is 
GQ, = 6.7 
By substituting §S = J 6.7 in Equation (2-43), the magnitude calculates 


to 


C 
eG 6.7)|= 707 


90 








2-7, Evaluation of phase angle. 

In the previous sections the magnitudes of the frequency response 
of a system are evaluated on the Bg vs By plane. Since the numerator 
and the denominator of the equation defining frequency response can be 
represented as vectors on the Bg vs By plane, it is obvious that the 
phase angles in the frequency response of a system can also be evaluated 
on the Bg vs By plane. 

In this section the method for evaluating the phase angles in the 
frequency response of a system is illustrated with some of the numerical 
examples which are considered in the previous sections of the magnitude 
evaluation. 

Again writing the general equation for the frequency response of a 


system as 
NG) 
(4.-B.) +; 0(4,-B,) 


: 


the phase of the frequency response is 


a ®) = 


ee a4) ae NG wd) Os. ty a (A, ~B,)) (2-45) 


Since the vectors N(j® ) and (a, ~ By) al - By) are measured 
from the M point and the real and imaginary parts of the vectors are 
represented on the Bo vs By plane in such orientations that they are 
parallel to the Bo and B, axes respectively, a proper reference in meas- 
uring the angles must be set up. 

Consider the angle Z((a, - By) +50 (ay - Bo) with the aid of 
Fig. 2-27. The (a, - B,) plot is constructed as illustrated in the 
previous sections with the I, curve and the M point configuration as 
shown. Select a point such as P in Fig. 2-27 at frequency of @, and 
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in the frequency response. 


a2 


__Bo = a, line 





draw the vector PM. The components of this vector define the real and 
imaginary terms (a) - By) and @,(a, - B,) in magnitude and sign. The 
horizontal component (ay - B,) is the imaginary part and positive when 
directed to the right. The vertical component (a9 - Bo) is the real part 
and positive when directed upward. Considering the signs of the compon- 
ent vectors, the angle of the vector PM = (a9 - Bo) +) Oi(ay - B,) is the 
angle 9, which is measured from the vector (a) - By) to PM, and this is 
equal to the angle measured from the line B, = a, below the M point to 
wP and positive when measured to the clockwise direction. The angle of 
the vector N({) is naturally B5 which is measured from the line B, = a, 
above the M point to the vector N(j® ) with the positive direction as in- 
dicated, thus the phase angle of the system in the frequency response at 


frequency @) is 
Le 4 ®,) = 6, - 8, 

It must be noted that the reference line of the angle of N(jo) de- 
pends upon how the vector N(j@ ) is represented and it is seen that the 
vector NC jo ) can be represented either to the right or left of the line 
By = a, and either above or below the line Bo = ag; however, for the 
angle of the vector (a = Bo) +) (ay = By), the reference line is always 
the line B, = a, below the M point and the positive angle is measured 
from the reference line to the vector iw in the clockwise direction, thus 
the 0°, 90°, 180 and 270° positions are as indicated. 

Evaluation of phase angles in the frequency response of an open loop 
system is also possible as the case of obtaining magnitude in the fre- 
quency response. The following examples illustrate evaluation of phase 
angles in the frequency response for closed loop systems as well as for 


open loop systems. 
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Example 2.8 


For the system equation in Example 2.1 which is rewritten below 


Ye 
6G) =—s (Sy (St3) 


obtain the phase characteristics in the frequency responses for both 
closed loop and open loop systems, 
Closed loop system 


The closed loop system has the transfer function 


@ = ¢- 
RC)” Score 


in which the numerator polynomial is of order zero and the angles 
vs N(j@ ) = 0 for all frequencies, therefore only the angles 
Z (Ca - By) + w(a, - BL) are frequency dependent. 

In Fig. 2-28, the Ts curve and the (a, - B,) plot are duplicated 
from Fig. 2-8 which is for obtaining the magnitude of the closed loop 
frequency response of the same system. 

The reference angles Cue 90°, 180° and 270° are indicated. The 


radial vectors directed from the M point to the G) (ay - By) plot give 


the angles of the denominator of the closed loop frequency response equa- 


tion which is the vector (4 - By) + j@(2- B,) . Noting that the angles 
measured from the 0° reference line to the vectors are of the denominator 
and also noting that the angles of the numerator are zero degrees for all 
frequencies, the phase of the closed loop system at each frequency is the 
negative of the angle of the vector so measured, thus at (z=: 0 the phase 
is zero degrees and as () increases the phase increases in the negative 
direction. For example, at @sz= 1.0 the angle of the vector measures 


45° and the phase of the system is -45°, The phase assumes -180° at 
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(Example 2.8) 
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qs 





a 1.4 and at this point the ie curve, ©(a, - B,) plot and the line 


By 


aj = 2 intersect. If a protractor is overlapped on the figure such 
that the center point of the protractor coincides with the M point, then 
one can read off the phase at each ® conveniently. 


Open loop system 


Rewriting the open loop transfer function as 


4 
G(S)= eae ares, 


one can construct the i curve associated with the denominator polynomial 
of the open loop transfer function. By comparing the equation with the 
closed loop transfer function, it is obvious that the li curve for the 
open loop system is identical with that of the closed loop system. The 
My point defined analogously as the M point is, by noting that the co- 
efficient of the s° term in the denominator polynomial of the open loop 
system is zero, at By = 2 on the By axis and 
Mp =M, (2, °) 

is located in Fig. 2-28 as shown. Since the My and M points are on the 
same line By = 2, and the IE curves for the two transfer functions (open 
and closed loop) are identical, the G@(a, - B,) plot is not altered. 

By applying the same principles to the open loop system the phase 
can be evaluated with the vectors directed from the Mo point to the 
@ (a, - By) plot. Naturally the 90°---270 line must be shifted down to 
the B, axis. Noting that as @® tends to zero the ®(a, - By) plot ap- 
proaches the Mo point from the direction of the origin the phase at 
@®=0 is -90°. As @ increases the phase increases in the negative 
direction. Some vectors representing the phase of the open loop system 


are drawn with dotted lines, for example at @ = .5 the angle of the 


96 





vector measures 130° to give a Sisto)" phase. 

At ® = 1.42 which is the point where the phase of the closed loop 
system is -180°, the open loop system also assumes -180° phase and it is 
interesting to note that this occurs for any system which has a frequency 
independent numerator polynomial. 

As @® tends to infinity, from the behavior of the in curve and the 
@ (a, - B,) plot it is seen that the phase approaches -270° which is the 
expected phenomenon. 

In the Bode diagram manipulation, the gain crossover frequency and 
the phase at that point are used for the stability check and as a figure 
of merit in a design problem. The gain crossover point which is the zero 
db or unity magnitude point can be found simply by drawing a magnitude 
circle with radius equal to ag = 4 centered at the My point, which is to 
account the magnitude of the open loop system, thus the gain crossover 
frequency is found to be 1.05 and the phase at this point is ~168° to give 
+12° phase margin. This phase margin indicates a stable but too oscilla- 
tory system which checks with the if curve and the M point configuration. 
The gain margin can be found by drawing another magnitude circle centered 


at the My point passing through the -180° phase point. 
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Example 2.9 
For the system in Example 2,2 


a [2;000(St%&7) 7 
6S) = SCS +t1)CS+1670) 


the transformed closed loop system equation is, by using the transforma- 


tion S = 100s 


ye eee 
‘ GI HIB TIS HA 2ETS HAGE 


In order to obtain the phase characteristic of the closed loop system 


the [e curve and the @ (a), - BL.) plot of Fig. 2-9 are reproduced in 


Fig. 2-29. 

The angles of the denominator vectors oy are measured in a similar 
way as in the previous example. Noting that the system has the first 
order numerator and the vector NC j%_) at each a, is represented by a 


vector directed from the M point to the line Bot = 2a. (the two vectors 


Ot 
for Ww. = .2 and 00, = .3 are shown in the figure with dotted lines) the 


a above the 


angle ENC j0,) at each 4), is measured from the line B an 


Neos 
M point to the vector N( s;W,) at each @,. For example at 
. = .2, LZ£NCj 2) = 25°, Denoting the angles of the denominator and 


the numerator vectors at the same frequency by Oy and a, respectively, 


then the phase of the closed loop system Le Cy @_) = 85 - 0, 3) for 


ca 


example at @, = .2, £=(5-2) BONete 31S > ya = 20508 


For the open loop system the transformed system equation is 


I2AB+tG.GF 
Orie DIB + 1678 


The ig curve for the open loop system is not changed from that of 


G(4) = 


the closed loop system but the Mg point is moved to the point (.167, 0) 


0 


which is seen from the coefficients of the ae and s terms in the 
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denominator of the open loop transfer function. This implies that a new 


W, C4y¢ - Bit) plot for the open loop system must be constructed with 


4? .167 in order to obtain the frequency response of the open loop 


system and the procedure to construct this plot is straightforward. 


Since the Mo 


will not give good resolutions of the angles of the denominator for 


point is so close to the origin, the W (.167 - Bit? plot 


small ,. if the plot is constructed in Fig. 2-29. This can be overcome 


by plotting the curves for small ),. on an expanded scale. 
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3. The fourth order system charts. 

In references 1 and 2, a universal third order system chart is con- 
structed on the Bo vs B, plane so that a single chart will provide solu- 
tions to any third order system problem. 

This is done by using a suitable transformation. Any third order 


system is transformed into a normalized system so that it has unity co- 


3 2 


efficients of the terms S~ and S” in the characteristic equation to give 


a set of fixed Bo vs By curves. 


In a fourth order system, assuming unity as coefficient of the st 


term in its characteristic equation, Bg vs By curves are changed by chang- 


ing either a3, the coefficient of the 33 


32 term. 


term or ay the coefficient of the 


Since for any polynomial it is possible to make any two of the coef- 
ficients unity (or the same numbers), by a suitable linear transformation, 
in a fourth order system characteristic equation any two of the coeffi- 


cients of ge > 3 and 5? 


terms which are the coefficients affecting the 
Bg vs B, curves can be made unity so that all Bo vs By curves of any 
fourth order system depend upon only one coefficient. 

For each value of the variable coefficient, there results a different 
set of the Mitrovic's curves, If the sets of the curves for some suitable 
pertinent values of the variable coefficient are constructed and if the 
behavior of the set of the curves for the variable coefficient changing 
from one pertinent value to others is understood, then the sets of the 


curves may serve as fourth order system charts which can be applied to 


any fourth order systems. 
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3-1. Transformation of the fourth order system equation. 
Suppose a fourth order system characteristic equation is; 
3 2 
F(S) = S* + 4,S°+4S +4, +4, (3-1) 


and associated By vs By curves are defined by 


3 
A, PSO + AgPlSWm + ALS) Wy, 


By = 
2 4 (3-2 
Pye, anges )ca, = ae ) 
Using a transformation or frequency scaling 
Gala 
Equation (3-1) can be written as 
a a a, ay anced, ww +- 4. =O 
F(4,4)= 4, 4 444° + AJ 4s sa 
or by dividing through by as 
2 a +“ 
$ (8) = 44+ 33 + 22/8; S°+ 8/4345 + 4,/ 4, = 
or $(6)= 5% + 674+ a, + Ato Ant =0 (3-3) 
where Bees das 
3 
dit =d,/ a3 
Ae =a,/ af 


Thus Equation (3-1) is transformed into Equation (3-3) in which the 


two coefficients of the st and 33 terms are unity. 


Equations (3-1) and (3-3) are equivalent since the transformation 


Ss a8 is linear and one-to-one. 


The Bo vs By curves associated with Equation (3-3) are defined by 


3 
Bit = Gat Fa (S Wye + Dy (S) Oye + D,C$) Wnt 
(3-4) 


% 
Bue = ALi D, (3) w,3 ot P(S) Wy 
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Equations (3-2) and (3-4) are equivalent to each other with the fol- 


lowing relations 


B.=B,/45 
Se =B,/a 
Daye = Un 123 


Comparing Equations (3-2) and (3-4) which are Mitrovic's equations 
for a system represented on different coordinate scales, it is seen that 
studying a fourth order system with the second equation is more convenient 
than with the first, because in the second, only one variable coefficient 
ao, appears whereas in the first the two coefficients ag and a3 appear. 

A fourth order system closed loop transfer function represented as 


a 
ae aS hs (3-5) 
R Soma S? + aes + 2,8 +12, 


transforms into 


bed + bits + Dot 


- (3-6) 
B+ BF + Ip b> t+ Ae d + ge 


g(4)= 


to which the Mitrovic's curves defined by Equation (3-4) associate. 
The coefficients and frequencies of the transformed system are related 


with those of the original system by 


be= b./a3 
be= b, yee 
Oe = Wy [4s 


plus the relations defined as in E,uation (3-3) for the coefficients of 
the denominators. 

The velocity constants, K, and K,, of the original and transformed 
systems respectively are related by 
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kei — Ky Jaz 


ap 


since a= oa 


aye 
ke? chy = Bastar = (Gab) /4s = hes 

Thus a fourth order system is transformed into a system from which 
resulting By vs B, curves are dependent upon only one coefficient ao,. 
Comparing Equations (3-6) and (2-38a), it is noticed that the characteris- 
tic equations are of different forms which indicates that the normalized 
ie curve in Section 2-6 is different from the [3 curve resulting from 
Equation (3-4) and any Ie curve defined by Equation (3-4) cannot be 
plotted on the normalized Bg vs By plane unless the system Equation (3-6) 
is transformed into the form of Equation (2-38a) and Equation (3-4) is 
forced into a form to give the normalized if curve. If one wishes to work 
on the normalized By vs By, plane with Ie curves, then the system Equa- 
tion (3-5) or (3-6) is transformed into the form of Equation (2-38a) to 
give the equation defining Mitrovic's curves as 

Bin = (Fe BUS) One + Py-S) One He GlOne 


(3-7) 


2 3 & 
Bon = “ne - BSI, - BCS Ane 


Equation (3-7) gives the normalized i curve, and the ie curves 
can be plotted on the normalized Bo vs By, plane; however, ao¢ appears 
in too many terms and more complicated behavior of the curves is expected. 
Thus the normalized Bg vs By plane is not suitable to study fourth order 
systems. However the normalized curve can do an excellent job in determin- 


ing absolute stability of a fourth order system. 
The coefficient a7_ appearing in Equation (3-4) indicates that Bg vs 
B, curves for a fourth order system depend upon the value of ap, , 
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therefore it seems worth while to investigate the possible range of Arte 
The coefficients a and as of Equation (3-1), the characteristic 
equation of the original system, can be written as 
a3 20% tot tet ry 


(3-8) 
ay ee ar rey, +e, + ee 


i 3 4 
where -1r,,-9, -t3 and -r, are the roots of Equation (3-1). 
The roots r; , rg, ¥%3 and ry, fit into one of the following three cases: 
1) All 4 roots are real. 
2) Two real and one pair of complex conjugate roots. 


3) Two pairs of complex conjugate roots. 


Noting that ao, = ay/as , by using a9 and a3 expressed in terms of 


the roots, it is possible to determine the range of ay, to cover most of 


the fourth order systems for Case 1 and Case 2. 


Case l 
a - fia tls + ila + lols & lee 
2t ” CEP its plein = 
file +t hi thle bls Hla ie tee 
Pret +h +2(ne neti ter Gn+G Ty) 
noting that iis r,* + ls + me 70 


CE 6 5 
It can also be seen that for a fourth order system, if its open loop 
transfer function has a numerator polynomial of order less than two, and 


if it has all real poles, then 
4226-5 


In such a system, the coefficients apg and a3 of its characteristic 


equation can be expressed in terms of the poles of the open loop transfer 


105 








function, Those poles of the open loop transfer function replace the 
roots of the characteristic equation in Equation (3-8) to get the same 
result. 

Case 2 

Let the two real roots be Ty) and Yyg , and the complex conjugate 


roots v3 and Ty, C 


TZ SWy +) fr s2 Yn 
T%] = $a), —j [3% % 

and aS = (Kit th + 2SWn) (+n y4 430, (rin) + 4s? 
a2 = le +2SWn(M +h) +a" 


to compare ag and ay 


A A= Pht ls +230,(r th) 14Se1)o, (3-9) 


Assuming a stable system, Equation (3-9) is examined for several 


situations. 
EE lryeee > Se | then 
2 2 2. 
a- 4, >(i1 $°-1) My implies that 
if 5 >. 302 then 4,- 4, >0 and << 1.0 


It must be noted that the inequality above may have considerable tol- 
erance which implies that even for § some value less than .302 a9, may 
be less than l. 

Let rpg > Yr, and by neglecting ry in Equation (3-9) 

2 
ay ~ a, > fy +2 Sn t, + (43%!) ay (3-10) 


For a well behaving system, rp is usually greater than $4), the 


real part of the complex conjugate roots. 
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Let rg = k§wW, where k is a positive real number greater than l. 


Substituting this ry into inequality relation (3-10) 
Z 
Ag - 4, >((e7+ 2k + 4) S71) Wy (3-11) 


Inequality relation (3-11) indicates that for a reasonable number k, if 
S$ is not too small, then aa - a. is greater than 0 and ao, <1.0. 

For example if k=2.5 then af - a,>0 for $7.25 and if k=3 
then aa > ay > 0 for <$>.229 

Thus for a fourth order system, if it has no more than one pair of 
complex conjugate roots of the characteristic equation and if it is ade- 
quately damped, then the coefficient aoe is less than l. 

For Case 3, which is the case of two pairs of complex conjugate roots 
present in the characteristic equation, it is hard to set up an upper 
bound of a5, as in Case 1 and 2. This can be seen as follows: 


For Case 3 


Let h=h , b=r, 
and hr ore 3,On) = fi-$? Wont s) Ig = 3. Wa -j firs Wy 2. 


2, 2 
then a; = (2 $, Wn =e So Wna.) 


a 
ad, = (Es ORE iyo, + Una 


he as dy? + $3.0) ho + Wha 
aol — SL a St aL ee 
" az 437th, + PS, $2 Wn, Wya + 4S Wye 
Let Das < Dn 2 and R= Wo / Das 
+ + R* 
a; “ZS, +83,S.R + “f# 
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a 
Thus the ratio */ a2 depends upon AT ; So and k. If both 5, 
and 55 are small, then move? assumes a large number. However, if a well 
3 


damped system is assumed, then Case 3 can be approximated to Case 2 and 
ao, < 1 is a reasonable upper bound, For example if $, = .3, $y = .7 


and @ , =Wy2 then = = 8 = afl 
3 e 


In this section, in studying fourth order system curves, aot <1is 
assumed. [It must be noted that aoe <1 is not the sufficient condition 
for a fourth order system to be well damped. Even though ao, is very 
small, a fourth order system may be unstable depending upon the location 
of the M point. [t turms out that one can find a necessary condition on 
ant for a system to have a relative stability greater than 31 a given 


system damping coefficient. This will be discussed later. 
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3-2. Three dimensional representation of fourth order curves. 
The By vs By curves for fourth order systems defined by Equation (3-4) 
depend upon a5, which is the coefficient of the se term in the character- 


istic equation of the transformed system equation. a4,, can assume any 


2t 
value within its allowable range and independent of frequency on the Bo vs 
By plane. 

If a new variable Boy = ao. is defined and if Bo, varies on Boy axis 
which is chosen to be perpendicular to both Bit and Boe axes to form a 


three-dimensional cartesian coordinate system, the equations 


Bzt = 2at 

3 
Bye = Ar B38) Wn + P($) On + Da ($) Wr; (3-12) 
Bot = an Qh - (5) ww, - BS) Wy 


define the Mitrovic's surfaces on the Bor ~ Bie - Bot coordinate system. 

For fourth order systems, Equation (3-12) defines fixed surfaces 
which can be applied for any fourth order system as the case of the third 
order chart for third order systems. 

Since the variable B,, is continuous and Bot VS By, Curves are con- 
tinuous, the Mitrovic's surfaces are also continuous. Noting that any 
plane Bor = ant = constant which cuts the surfaces gives the set of Bot vs 
By, curves for various values of §€, by applying the continuous property 
of the Mitrovic's surfaces, all of the basic property of the Bo vs By 
curves and the analysis and design principles of feedback control systems 
on the By vs By plane can be extended to the three-dimensional Mitrovic's 
surfaces. 

For example, the stable region for the M point on the By vs By, plane 
which is the region in the first quadrant enclosed by the ie curve can be 
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extended to the three-dimensional Mitrovic's surface as the region (volume) 
in the first octant enclosed by the ig surface which is the surface made 
up by the lS curves for all values of 45, 7 0. Locating the poiat 


M( aot) on By, vs Ba, plane is to locate the point M(a5,5 aie: age) in 


oes 
the Mitrovic's volume. 

Fig. 3-1 is a schematic sketch of the Is surface. B5, axis is chosen 
so that Bo,, By, and Bor axes form a right-handed coordinate system. The 
equations defining the [2 surface are 

Bor = Aer 

yen any cs (3-13) 

Bor = A Wnt a Be 
which are the equations for a parabolic surface. The surface forms a 
cave with closed end at the origin on the Bor = 0 plane. Applying the 
stability criterion to each Bot vs Byz_ curve which is formed by cutting 
the surface with each Bot = 805 plane, the stable region is the whole vol- 
ume inside the cave. 


If the surface is cut by a plane Bot = Age = Constant, then there re- 


sults a By, vs Bo, (or B, vs BL curve for €s= 0 on the plane and like- 


Pas 
wise for the Bot VS By, curve. 

It is interesting to note that, Boy vs Bg; curves for $= 0 on Bor 
vs Bo, planes for fourth order systems are straight lines whose slopes and 
axis intersections are determined by a,, only and an extremely simple 
method to test the absolute stability of a fourth order system can be 
derived. 


By eliminating Oat and as, in Equation (3-13), the if surface of a 


fourth order system is 
2 
Bor= be Bee 7 err (3-14) 
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Fig. 3-1 The e surfece for 4th order systems. 
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By cutting the surface represented by Equation (3-14) with a 


Bit = 4, Plane, the Bor vs Bay curve for § = 0 is obtained as 


2 


Boe = Zit Bee it (3-15) 
which is an equation defining a family of straight lines om the Bo; vs 
Bo, Plane (or By, = 4,4 = constant plane). 

For each value of aie? there results a different line. By reason- 
ing with the stable region in the three-dimensional volume, the region 
for the point M (a,;, age) on Bot vS Bo; plane to guarantee a stable 
system is the area in the first quadrant below the line defined by Equa- 
tion (3-15) or mathematically 

O< a < Bit » Gnr —Qie 
which agrees with the Routh criterion. 


The lo curves for some values of aj, are plotted in Fig. 3-2. As 


t 


the Equation (3-15) indicates, a i line for a,, = k cuts the Bo, axis 


lt 
at Bo, = k and it has a slope of k. The frequency on these i lines is 


constant at all points on a line which is a peculiar nature of the ie 
curve on a Bo; vs Bg; plane. Thus the iy line itself is a constant Wrt 
curve on the Boy vs Bo; plane. 


For instance, along the i line for a), = ays Ont = 446 which is 


equal to Jar as the second of Equation (3-13) implies. Thus a point 


t 
lt > Slee 


represents a closed loop system which has a pair of pure imaginary roots 


M (ao, » ae) located anywhere on a single is line, say for a 


r= + jo, = J ay't which is readily seen from the basic nature of the 


yj, curve. 
As an example, suppose a characteristic equation of a fourth order system 


is 
$6) = 5° +74 Arb + 24 + agp =o 
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Since ae? ~2 the ee line is the line labeled as a,, = .2, 


It 
G@, = -446 in Pig. 3-2. The Mitrovic's working point M (age 2 49¢) lo- 


t 
cated anywhere above the ts line represents an unstable system, and 
dually for below the line. If the M point is followed along the line 
then resulting every characteristic equation has a pair of imaginary 
roots rs + y 446. For instance, if M (.7, .1) is chosen, then the 
characteristic equation is 

F (A= 44° + S74+.98 + 28+. =0 

and Fi &j.446)=0 


The point M (.7, .1) is also on the I line for ay 5, therefore 


t =- * 
another characteristic equation 
f(o=4%+ poe esi So 54 +.{=0 
has the roots fr = tvs and it calculates to 
fa Gi/F)= 0 
From the Mitrovic's basic equation and the way of defining a three- 
dimensional plot, the concept of a Mitrovic's surface can be generalized 
to any order system and to a set of any three coefficients of the charac- 
teristic equation. Then fixing one of the three axes to a constant, 
there results a well-defined curve which is the one defined by the gener- 
alized Mitrovic's equation in Reference 3, though the generalized Mitro- 
vic's equation turns out to be undefined and uninterpretable at some 
particular values of ¢ for some particular set of coefficients. (For 
example, the IS curve for the Boz vs Bo- plot.) 
The most significant feature in the concept of the Mitrovic's sur- 
face is the possibility of controlling three coefficients simultaneously 


since the surface is independent of the three coefficients which are to 


be controlled. 
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3-3. Study of fourth order curves. 

As pointed out in Section 3-1, a fourth order system can be trans- 
formed into an equivalent system so that there remains only one param- 
eter (coefficient) in the characteristic equation on which the fourth 
order system curves depend. This feature greatly simplifies studying 
the behavior of the Mitrovic's curves for fourth order systems. 

In Section 3-2, it was shown that if a fourth order system charac- 
teristic equation is transformed such that the coefficients of the s4 
and 53 terms are unity, then the coefficient of the $2 term a9_ which is 
the only one affecting Bo; vs By, curves is less than 1.0 for most prac- 
tical fourth order systems. 

Since the Bo, vs By; curves are continucus for all values of rue 
6, and the coefficients of characteristic equation, if the behavior of 
the curves is well investigated and if the curves for some pertinent 
values of parameters are constructed, then those preconstructed curves 
may serve as a guide or charts for both analysis and design of fourth 
order systems 

In this section, the effect of aj, on Bo, vs By, curves is studied 
to provide the fourth order charts. 

In Equation (3-4) which is the equation defining Bot vs B), curves 
for a fourth order system, there are two parameters a9, and 6. For one 
value of aj,, each SG gives a different curve fe to form a family of 
curves corresponding to the particular value of Bape This family is de- 
noted by [het are » or simply [sep . Then, for each ay, there re- 


sults a different family. 
Fig. 3-3 to Fig. 3-22 show such families. On each sheet the curves 


for five values of § are shown. For one value of a9;, there are 
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associated two sheets of graphs, one for 0 < § ¢ .5 and the other for 
5 £$< 1.0 with interval of § .1. The values of a9, are from .1 to 
1.0, thus for every two sheets of graphs, aor increases by interval of 
-1. (Note that for 4 < ao, = 1.0, the first graph for each value of ao, 
contains the family re for 0 £§= 1.0 and the second graph con- 
tains the family re for -5 £5 £ 1.0.) The constant wt curves 
are also shown for some pertinent values of ott 

It is observed that as aoe increases, the curves are in general ex- 
panded in both abscissa and orxinate. This is well expected from the 


equations defining Bo, vs Byte curves. The equations are rewritten below 


3 
Be = Fat PalS) Une + BCS) Woe + ACSIA 


(3-4) 
3 $ 
Boe = Gat Wnt — BlS) Wye — PS) One 


The terms with which aj, contributes to By; and Bg; are non-negative 
for 0 <6 £ 1.0, thus increasing aot increases B;, and Bg;- 


For small values of ay each member of the family re is so 


t 
shaped that the whole region under the curve in the first quadrant is 
the region for the M point where the relative stability is guaranteed to 
be greater than the value of % of the curve. Such a region is denoted 
by Re - For example in Fig. 3-3, the whole area enclosed by the i 
curve and B,, axis in the first quadrant is the region R,|5 on which if 


the point M (a ) is located, then the system has S greater than .5. 


It? “oe 
AS a5, increases, the manner in which the Re region is formed is 

changed in such a way that the curves for high 5 do not form Re : 
The behavior of the curves depends upon the values of a5, and & . 

Examination of Equation (3-4) reveals that for a particular value of 


a there is a certain critical value of 6 such that if 6 is increased 


20 
136 





iat 
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beyond the critical value, then the le curve exhibits counter-clockwise 
encirclement on the Bot VS Byte plane so that there is formed a small Re 
near to the origin or no Re at all. Equivalently for a particular 
value of § there exists a critical value of apt. 

Investigation of the critical value of 6 or aot is important since 
in synthesis of a system one must choose a point M (aj;, 4g¢) so that the 
specified system performances are met and particularly for the relative 
stability of the system, the M point must be located on a Re where 6 
is some suitable value to meet the specified relative stability. For 
example if the system relative stability is specified by 4 = .5 then 
one should not choose the M point anywhere on the Me curve for 
ao, = 1.0. (See Fig. 3-21 and 3-22.) 

In synthesis of a system, usually 6 is chosen at some pertinent 
value according to the required relative stability. Therefore it is con- 
venient to investigate the critical value of ap, for given values of 
6, thus a study of the behavior of the curves with changing aj, for a 
fixed 5 is needed. This requires that the family of curves must be 
represented in a different fashion from those of Fig. 3-3 through Fig. 
3-22 in which the families of Boy vs By, curves with changing § for 
fixed aye (the family [sst,45-) are represented. That is, the family 
of curves must be represented with changing ao, for fixed S . Those 
families are denoted by lefaoet A 

At this point, in connection with the study of the families 

le tang » a constant € plane is defined as follows: 
A constant ¢€ plane is a Mitrovic's plane on which Mitrovic's curves 


for fixed 6 are plotted, and denoted by the $= k plane where k is a 


37 





constant andO << k<¢1.0. With this definition, for th: eh 
systems, the family of Bor VS By, curves with various valves of a 
G fixed to zero consists of all of the Bot VS By, curves for < : 
the “% = 0 plane, and likewise for any other values of 0 = € € 1.0. 


The constant @, lines on the constant <4 plane 


In general, if the Mitrovics curves are plotted on a constan® ¢ 
plane, and if any one of the coefficients which appears in the eguatic 
defining the Mitrovic's curves is varied, then all of the points at which 


the frequencies are the same are on a single straight line whose slope 47:3 
the coordinate axis intersections are determined by 5 and WwW. this 
line is called the constant W, line. (Similarly the constant oe ine 
for a transformed system.) 

For Bg vs By curves defined in general form as in Equation (2- 
by fixing 6 = €, », a constant, the equation of Bo vs By curves on the 
Si plane are: 


1“ 


= k-| 
B= 2 Fa Pe(S) Op 
Kka2 
(3~ 16} 
woe a 
0 BAe Fe (5) an 
If the coefficients a,8 » k#i are fixed and if a; is varied, 
ieaoat 2 
then on the $4 plane, a family of Bo vs By curves |e, fai} is formed 
Rewriting Equation (3-16) to separate the term containing the vari 


able coefficient ay from the terms of the fixed coefficients, 
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4- n &- 
B= OS) Wn em Ap Pel B® 
RL 


(3-17) 


oe. id 
By=-4, 4. (Dn ~ 2 Ae Pei.) Om 


=, 
. 


RFA 


Suppose any two members of the family ls tan , one for a, = ary 


and the other for a, = ajo + Let the increment in a, be 4a,= ayo 7451 
then the increments in the coordinates B, and Bo at any frequency asso- 


ciated with Aa, are 
ie 
AB =44, 616) Wn 


A B= -44; G1 05.) Wn (3-18) 


AB, 


ee — Dkt0 On SPs (51) 





ww 


Equation (3-18) implies that for any change 4 a, in the coefficient 
ay the ratio of the changes in the coordinates 4 Bo/ 4 B, is independent 
of a; at all frequencies, and that the distance between any two points 


at the same frequencies one on the curve for ay and the other on the 


curve for ay, is linearly varying with Ma; . 

For the family ex on the O41 plane, consider a point at 
Wry = Oni on each member of the family. If a, assumes a continuous set 
of values, then the locus of the point of a, = Day on the o1 plane is 
obtained by eliminating ay in Equation (3-17) and replacing O), by Ont ; 
namely 


= dj. (3: WD» n £-1 n & 
p,=- rae Dy Pe (S:) Oyj )- 2 Sheng (3-19) 


. 


RL REL 
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Des(8.) Ons B+ agus Ble )= B(5.) Hou) (3-19a) 


or a= 
ae $14, 


2 P: (51) 


Since 6 is fixed to 34 on the 4, plane, and @, is fixed at W, , 


Equation (3-19) describes a line with a slope of — Pel $.)QDny pass- 
P($:) 
é <1 ee t 
ing through the point (z A, P.(5:) Wn » 25 Bg P_(41) wnt 


Refi 
on the ey plane. Thus the equation defining the constant 4), line for 
(n= a1 on the eal plane with carying ay is obtained. 


The axis intersections of the constant a line are found by letting 


B, = 0 and Bo = 0 in Equation (3-19) or (3-19a) 


a8 £-| AS) a : R-1 
B, intersect, - Z, de P(E Oni 2 BE) eae. Q.($:) Dans 


1 h ; , k 
Rei ae eri 
a RA DiS.) Pol. - BIN Hs 
=e Ay, Des) EI BENE) (3-20) 
$a om ( FC) 
4 
Bo intersect. = = A Ons (ee) (3-21) 
4 4 


pei 


In particular if i= 2 then noting that -hQo.1-Fp-2 = Py 


and Py = -l, Equations (3-19), (3-20) and (3-21) reduce to 
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(3-285 


For a fourth order system whose Bo vs B, curves are defined by Zqua- 


tion (3-4) which is for the transformed system, the constant Dat lines 


for Ont = Wer on the oy plane with aoe varying is obtained from Equa- 


tion (3-19) or (3-22) by proper substituticns, thus 


Wet fo, eee .. 
Bot = Bye a, (i Oe wWa.)} 


with axis intersections given by 


cs eee S 
Byy intersection = i, (3 Co ea | 


Bot intersection =- ao (I= £6201) 


Equation (3-25) defines the constant wW lines 


nt 


Oe on the S41 plane for fourth order systems. The 
stant Ont lines are proportional to the frequencies 


tional to S96 Gr, = 


all constant OD ot lines on the 


1412 


(3-255 


(S226) 


(3-27) 


for all values of 


slopes of the coa- 


and inversely propor- 


26). For & = 0, the slopes are infinity 


5 = 0 plane are vertical Lines. 








In view of Equations (3-26) and (3-27), if @,;< 2%, he 


Bit intersection > Q 


Bor intersection < 0 and dually. 


For the maximum values of the B intersection and einizus value 
Lt 


Bo, intersection, differentiating Equations (3-26) and (3-27) give 


Max, By, intersect, = 1/278* at OL le 
Min, Bop intersect. = -27/20485° at Wei= 3/85, 
Thus the constant WW, line of the fourth order system moves to the 


right with increasing slope as @,, is increased. Then at uw), = 1/35 
the maximum By, intersect, is reached. As the frequency a, is increased 
beyond 1/35 then the constant G),e line moves back to the laft with = 
steadily increasing slope. 

Fig. 3-23, 3-24 and 3-25 are families of Bo; vs By, curves 
eo plotted on the Ge= Q, .1 and .2 planes respectively. The 
coefficient a9, is varied from .1 to 1.0 with Maz, = .1, The constane 
Ant lines are shown some pertinent values of frequencies, and those can 
be compared with the constant ©); curves in Fig. 3-3 te Tig. 3-22 which 
are the families ror . 

The property of the constant OD - line can be used in constructing 
& Bop VS By, curve on a constant plane for any value cf Bn when there 
are on the plane any two curves each for a known value of age. This is 
illustrated in Fig. 3-26. 

Suppose two ie curves, one for a, = .30 and the other for 
a9, = -40, are preconstructed on the S= .5 plane. The lire seguents 
joining the pairs of two points at the same frequency on the twe curves 
are the constant W,, lines. for example, the points Py and Fo are af 
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a Boe vs Bay curve by using 
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the same frequency ae = -25 on the two curves, one for 5, = .30 and 
the other for a>, = .40 respectively, end the line Peareae joining the 
points Py and Pg denoted by L 9, is the constant @,, line forW,. = .25. 
All members of the family PS en will have the WJ. = .25 points ca 
or extension of the line segmeat Li95° 


7 
Suppose the ls curve for ay, = .32 is to be constructed. Since 


t 

the increments By, and “Bo, associated with “a9, are proportional t« 
AS ao, at every frequency, the point Q which is located at one-fifth way 
from the point Py to the point Po on the constant {,, line L 95 is the 


ex 


5 curve for ao; = .32. Any other points ot 


point at Wat = «25 on the 
the curve I's for a9, = .32 are determined analogously, then connecting 
those points with a smooth curve, the desired curve is obtained as shown 
by the dotted curve. 

With the aid of the constant We line, one can visualize the changes 
of Bo, vs By, curves with changing a>,. For example a point P, at 
Ont = .25 moves along the line L.g5 in Fig. 3-26 as a-, increases frem 
-3. If a point is located on a constant Git line for Ot = a, on the 
Ss S$) plane, then the value of ao, for which the le, curve passes 
through the point can be determined from the linearity property of the 
constant @), line. 


baer 


Now consider an M point located at ore . Oey on the Sy curve, 
and one asks if the M point insures &% 7 for the system. For the 
problem to be specific, suppose the M point is chosen at ‘, = .5 on 
= 
ls curve for a9; = 1.0. The family 51a e is shown in Fig. 3-27. 


By applying the basic rule to determine system stability to the point 


M (.375, .125) with the | 5 curve for ay, = 1.0, it is obvious that the 
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system has a pair of complex conjugate roots with & less than .5. 


It is seen in Fig. 3-27 that, for ao, less than or equal to. 7, see 


2 


entire areas enclosed by the Ve curves in the first quadrant are the 
R,5 region. For a9, = .8 and .9 the areas enclosed by the curves are 


divided into two regions, one the R.5 region and the other the nen=R 


‘) 


region. For aot = 1,0 the curve gives almost no R «5 region. 


t 
The curves in Fig. 3-27 agree with the curves in Fiz. 3-13 to 2-27, 
For G= .5, the critical value of aj, (mentioned previously) seems t- 


lie between ant = ./7 and aot = -8. If this critical value of ay, is de- 


et 


termined, then the M point chosen at any point on the rs curves for ay, 
less than the critical a5, will guarantee the system 5 greater than or 
equal to .5 whenever the M point is in the first quadrant. In the fel- 
lowing paragraphs the critical value of ao, denoted by aot, is investti- 
gated for each ©. 
The critical value of ao;. 

Inspection of the figures shows that in order for a I. curve to 


form a region R the curve must encircle a point in the region RL in 


oo Ss 
the clockwise direction. This implies as one possible situation that 
the points of the extrema of the coordinates Bit and Bo, must occur in a 
definite order; the maximum point of Bo; followed by the maximum point 
of By, (and then minimum Bot followed by minimum Bits 

Let Oa = and Opes be the frequencies at which By, and By, are 
maximum for a le curve. Then in order for a : curve to form a regton 


R the condition 


S D 
Womax < Wy max 


must be fulfilled. (Note: even if the condition Op, < Gu, wax is not 
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. 


fulfilled, there may exist R, and this is considered latar.} Tt 


S 
illustrated in Fig. 3-28. In Fig. 3-28a, Wonax < oo ae and the region 


R,. is formed, whereas in Fig. 3-28b, Womax > ymax and there is no 


§ 


region R 
ns 


and the region R, is formed. 
2 


- In Fig. 3-28c the critical case Wonax = Dimax is shown 


Since Go, and Dy ase are functions of aj, and § in the fourth 
order curves, the conditions Ce 2 Q) nae can be expressed in terms 
of ay, and & thus for a given oy the critical value of ay, car be 
determined. 

The values of Wimax and Womax are found by differentiating the 
equations defining the Bot vs Blt curves Equation (3-4) which is rewrit- 


ten below. 


' aS, 
B= doe BS War + PSOne + DS One 
(3-4) 


Bot = Ae (SD. — BB) Ome 
Noting that @(s)= 0 at §=0, Qs)= Oat S- .5, DKS) 
= Oat S= 0, and .707, in determining the conditions Uh,,,. = Ole 
it is convenient to divide the interval 0 ¢ §$ © 1.0 inte 3 sub- 
intervals 04¢< .5, .5< S 4.707 and .707 <§< 1.0. The epe- 
cial cases $= 0, .5, .707 and 1.0 are considered separately. 
Differentiating B,, and Po, with respect to Ont in Equation (3-4) 


and setting to zero 


OB i. 2 
Ta, OC er Gat Pp, + ZPOne tIAWye =O 
(3~28) 
Bh 2 sf ’ _ 7 aa = 
qe Se tre, oo Bone 4B UW,» 0 
ant 
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Fig. 3-28 An illustration of behavior of 


Bot VS Fi, curves.(The region R, and extreme points) 


$ 
vos 





Lf QR#zO and #0 , then by solving each of Equation (3-28) 


for Wnt 





(3-28 


& 
Sexe? 








where oO, is the frequency at the extremum of 3,, and like- 
wise for Wp,. 
For the interval of §, 0<6§<.5, P,>0, Py? 0 and P< 0. 


Therefore 
Pp; pee dar QD 
@ Been | 73 
Imax “39, | Ph 3B 








(3-30) 


_ 38 GR 24x 
Wo max "PD Tleeor* “dh, ; 





Note that the negatives of the double signs in Equation (3-29} give 
Oimin 29d Dgmin and those are negative numbers, and also in Equation 
(3-30) Wimax and Yoo, are positive numbers. 

Next one must find the conditions on a9; which are equivalent to the 
conditions Wonax é Oy eae. This is done conveniently with the aid of 
Equation (3-23) instead of imposing the conditions Op max SW max On 
Equation (3-30). 

Since Dace: and QU yas are solutions of Equation (3-28), the sub- 


stitutions 6) Wimax and Wit = Womax into the first and the second 


ie 

of Equation (3-28) respectively must satisfy the Equation (3-28), 
a 

thus ae Pn + 22,0, max + 2PalWsmazr) =O 


2. 
2492-3 PW max -4¥B(Wemax) =O 


hey 





solving for aj, and completing the squares 





_ 3% be, 


ae 
yor 


us 


la] 
ws 


a 
ry 152+ 2250, 
a,1=2R ( Wom ax + UP ee ae 
Equation (3-31) gives the relation between a,, DW) max 29d Dona. Not- 
ing the signs of Po, b3 and Pus and -3f4 > 2% for 0 << 3 
oo 
Equation (3-31) which are parabolas can be sketched as in Fig. 3-29. 


The two equations of Equation (3-31) are plotted on the same co- 


Both parabolas pass through the origin and are concave upward since 
-3 fi! Po moO, 2 Ps, >0 for 0<6<.5. The vertices of the parabolas 
2 
are below the frequency axis for both since Py < Q and 
ie 3BPu 


-2.25 D7 /5 P, < 0 and they occur at @),,y = - B/34,>0 and 


Womax = -/5Q Jak <0. Thus the conditions 


Womax = 1 max 


vill 


can be expressed in terms of ag, and ao,, as ae Coe 
where a9;_, is the value of ay, at which Wo max = Wiiaye sey tettiss 


Dee = Dimax = Womax» Equation (3-31) forms a set of two equations 


with two unknown to give 


a = ie 
eB E WSs) 
63-32) 


jae 
ai Sees 
By proceeding as above, it is easy to see that Equation (3-32) holis 
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(The 2nd of Eq. 3-31) 


OMSK = a Mare WwW. 
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Fig. 3-29 The relation between 85407 O ay and oe 
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for the intervals .5< 6 < .707 and .707<§ < 1.0. This is because 
Equation (3-31) holds for any of the 3 intervals of 5 provided W,_,, 
and Womax defined by Equation (3-30) exist. Therefore for the two in- 
tervals of $,.5 <6 < .707 and .707 <$§< 1.0 it is necessary to 
check if the value of ate determined by Equation (3-32) gives positive 
real Gy max 2nd Wonax: 

For .5<6< .707, A< 0, D, < 0, therefore from Equation (3-30), 
Wimax is real positive for all ag_. In order for Wy,, to be positive 


real ao, is restricted by, noting that b< 0, 


eo G on i er) 

Pi 32Q, § (457-1) / 

The expressions (3-32) and (3-33) are compared as follows: 
@ a 
“ ‘ 

a enone < atl i a where .5 < § < .707 

dtc 8<*(1-$*) #(4S7—1) 
Since [= 6" So ese eee 


g 5 
(esas \( cee ee Gen (ae ) 


? ¥ é 
-$6°r 66 eG = ee ome 
2 
66 36%436"-1 = 0 
2 
*-1) £0 
2 Ae: 
Since eee oe (<*-1) <0 
95% 
and Aree < se =) 


thus ao,, defined by Equation (3-32) gives positive real Or max and 
Qomax implying that Equation (3-32) is valid for the interval 
P2707, 
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The validity of Equation (3-32) for the interval .797 < ¢ < 1.0 
can be proved Likewise. 

For § = 0, .5, .707 and 1.0, Equation (3-32) is also valid. If 
6 = 0 the iF curves for fourth order systems form the region Rg for ail 
values of aj, implying that a9. ——> o which agrees with Equation (3-32). 


If G = .5 then Equation (3-28) becomes 


Sze73Oye =O 5 S4n4-3BUne =O 


to give @, Lin) alee faz, / 3 


Oymax = DE fe 

Let Wimax = Womax>» then agee = .75 which is equal to the value of 
arte evaluated from Equation (3-32) with S = .5 

If § = .707, the solution of Equation (3-28) obtained likewise 
gives a5,, = .9 which is equal to the value of ao;, evaluated from Equa~ 
tion (3-32) with S = .707. 

If $= 1.0 then Equation (3-32) is undefined, but lim(ay,,) = 
.375 which agrees with the value of ag+, obtained by solving Equation 
(3-28) with S= 1.0 

Thus a general expression of ajtc as a function of 5 is obtained. 
Table 3.1 below which well agrees with the curves in Fig. 3-3 to 3-22 
shows 29t, and Gen for some pertinent values of S - The Boe and Bie 


are the coordinates of the point at 4.,; = Wee on the le curve for 


“ot * “2te* 
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.05 
.10 
ole 
. 20 
925 
- 30 


739 


.70 
75 
-80 
85 


.90 


1.00 


Table 3.1 


S2te Wee 
» 50250E+02 -49881E+01 
- 12750E+02 »24775E+01 
- 58056E+01 - L6347E+01 
-33750E+01 1209 7E+01 
. 22500E+01 -95238E+00 
. 16389E+01 . 77938E+00 
~12704E+01 -65486E+00 
- 10313E+01 - 56090E+00 
-86728E+00 - 48 748E+00 
- 75000E+00 -42857E+00 
-66322E+00 . 38032E+00 
-59722E+00 - 34014E+00 
-54586E+00 . 30621E+00 
-50510E+00 »27723E+00 
-47222E+00 »25225E+00 
-44531E+00 . 23053E+00 
.42301E+00 -21151E+00 
-40432E+00 -19475E+00 
- 38850E+00 . 17990E+00 
-37500E+00 .16667E+00 © 
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Boe 
-62499E+03 
. 39050E+G62 
7705 LE+C1 
~243202+01 
.99187E+00 
-47532E+00 
. 25443E+00 
. 14764E+00 
.91110E-01 
-99038E-01 
 39813E-01 
-27762E-01 
« 19923E-01 
. 14661E-01 
. L1033E-01 
-84696E-02 
-66209E-02 
-52617E-02 
~-42450E-02 


-34722E-02 


Bic 
- 29000Er 02 
-62490E+91 
o 277 58E¢O1 
~ L5594E+01 
99 584E+00 
-68930E+G9 
~ 90423E+00 
- 38401 E+00 
- 30L59E+60 
~ 24271 E+ 06 
. 19928E+00 
- 16642E+00 
- L4102E+00 
~ L2LO7E+CO 
~ LO516E+00 
-92325E-01 
-81865E-01 
. 73264E-01 
-66135E-01 


-60185E-01 





The necessary condition on the coefficient a?t for existence of the 
region Re on the Bot vs Bit plane and the design charts. 

As seen in the derivation of ao,, in terms of 5 , the le curve of 
a fourth order system for a given § on the Bot VS Bry plane has the 
points of extrema whose locations depend upon the value of acy with re- 
spect to a9, and the encirclements of the curve, which are closely re- 
lated to the existence of the region Re » depends upon the locations vf 
the extreme points. 

Another feature related to the existence of the region Re is the 
behavior of the curve at large values of Wate This depends upon the 
value of $ and will be described with the basic equation of the curve, 


Equation (3-4) which is again written below. 


2 & 
Bit = Fae BSW + PSOne + PylS) Ont 
(3-45 
\ 


Bee ala - 0 (2) Oe 
or 0's s < 25), Q(S)> 0, ZB(s)> O and G(€)<0 indicate 
that the curve tends to infinity in the third quadrant. From Equaticn 
(3-29) it is seen that there are no G); .4,, and Womin: Therefore if 
a9  ao¢- then the curve vents down in the first quadrant and enters 
into the third quadrant after passing through the fourth quacrant, thus 


the entire area enclosed by the curve in the first quadrant is the region 


(th 


Re + If ay, is increased beyond a>-,. then the curve encircles in th 


$ 


counter-clockwise direction, and there exists the region Re provided 


i) 


curve does not enter into the second quadrant. Referring to Fig. 3-28b 


which is the case of a9; > a1, for 0<§€ .5, if ay, is not to» 
much larger than ay-, then the curve enters into the third qua?ratt 


ms 
cas 


through the fourth quadrant after forming a region R, in the 
a 
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quadrant. (Note that in this case the curve forms a closed locp in the 
first quadrant and the area enclosed by the loop is non-Re -) Ef a 
greater than the value of ap,, such that the curve enters into the third 

quadrant by passing through the origin, then there exists mo region R . 

This value of a5, denoted by a9+, can be found by setting By, = Boy = 

in Equation (3-4) and solving for aj; , thus one obtains ay in terms of 


to 
S$ as 
22, 


aoto= 1/45 
Wey = 1/25 Be 


where Deg is the frequency at Bit = Boe = © when ang = 


=2te 
Thus for 0<6<¢.5 the condition on ao, for the existence of the 
region Re is O< aon < aoto° As ao, approaches Bato the region Re 
gets smaller and it is in the region near to the origin. 
Noting that Equation (3-34) holds for 0< 6 < 1.0, and by comparing 


Equations (3-32) and (3-34), it is seen that 


Bote < axt0 for OPS 2 aeg 
oer 4.40 at 6 =. 707 (3-35) 


The relation between Re » a, and a5, for the intervals of Sea, 
5 <6 <¢ .707 and .707 <6 < 1.0 can be analyzed in a similar way as 
for the case of 0<6< .5 

Fig. 3-30a and 3-30b are sketches of the curves for .5 < 6 <.707 
and .707< 6 $ 1.0 respectively with various values of aj,. The curves 
are labeled by numbers from 1 to 5 and ao, increases with the labeling 


number. The arrow on each curve indicates the direction cf increasing 


Wnt» The dotted portions of the curves merely indicate the behavior of 
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the curves after leaving the first quadrant. 

For .5<6£.707, AO)>0, BsH<o ad Ale< 
indicate the curves tend to infinity in the second quadrant. Frem Eque- 
tion (3-29), it is seen that there is no positive real frequency of 
64min but there may exist Womin for .5<5£.707. For apy < ay 
which is the case of curve 1 or 2 in Pig. 3-30a the curve bends dow: in 
the first quadrant, enters into the fourth and then inte the tify 
drant, thus the entire area enclosed by the curve in the first quadrant 
is the region Re a Abe ate < Fat < A2to then as shown ty curve 3 ix 
Fig. 3-30a, the curve encircles in the counter-clockwise direction in thse 
first quadrant and enters into the fourth and the third then into tos 


second, thus a region Re is formed below the earlier part of the curve. 


If aoe 2 20t9 (curve 3 or 4) then the curve does not enter inte the 
fourth quadrant and no Re is formed. Thus for .5 <4 €.707, the rela- 
i i / < € Bale 

tion between Re > arr, and Ar. is the same as for O< S$ 


For .707<$<¢ 1.0, the curves are sketched in Fig. 3°-30b. Noting 
that a9t9< 22te >» it is seen that if at< 22t6 (curve 1 and 2) then the 
entire area enclosed by the curve in the first quadrant is the region 
Re a ake Beg < 82t < A7t¢ (curve 3 and 4) then the curve forms a close4 
loop in the first or/and the fourth quadrant encircling in the clockwise 
direction to give the region Re inside the loop in the first quadrant. 
The curve 5 does not form the Re region, thus for .707<464 1.0 the con- 
dition on a5, for the existence of the region Re is Q< a5, ¢ Bee 

In particular for § = 1.0 ag;¢s .375 from Table 3.1 and the region 
Rj 9 is the region for an over-damped system. Thus the necessary comdi~ 


tion on az¢ for a fourth order system to have an over-Camped region 
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Ry.9 is O< aoe < -375 which agrees with the condition derived from 
the statement in References 1 and 2, that is, 
The necessary condition of existence of an over-damped region 
for the M point is to have n-2 distinct positive real roots of 
the polynomial F"(-S) = 0 where F(S) = 0 is the characteristic 
equation and the primes denote differentiations with respect to S. 
This can be shown as follows; 
For a fourth order system, the characteristic equation is 
f(4)= B94 67 + aed? + aed + A, =O 
and £4) = (eee ee Ra, =0 


The roots of the second equation above are 


Ai =. 2B (ass) ane 


In order for the roots to be two distinct positive real 


a - 
4, <(-25)(6) =-375 


The analysis of the relation between R and aj, made above can be sum- 


S 
marized as in Table 3.2. 

Fig. 3-31 to 3-44 show the families leaa on the o= «35 04, 
«5, .6, .7, .8 and 1.0 planes for aj, up to about ay... (For $= .3 
and .4 plane a+ is up to 1.0.) Since in a design problem, the value of 
S is chosen deliberately at some pertinent values, and since the curves 
for any intermediate values of a9, shown in the figures can be constructed 
easily by means of the constant @)., lines (as shown previously in Fig. 
3-26) and since a5, = ag¢_ is sufficient for a curve to form the region 
Re » the figures may be used as design charts for fourth order systems. 
An inspection of the curves reveals that the arguments developed for 


Aztec » Fat, and the region Re agree with the curves. 
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Table 3.2 


(agtc is given in Table 4.1 or by Equation (3-32) 
agto is given by Equation (3-34)) 


Interval offs 


Ole 6 < .707 


82tce < 2940 


mio7~S = 1.0 


a2te> a to 


PIE 


O< age € 42tc 


Forces FatS 22t0 


49¢ Z 7240 


O< aor < 450 


9940 < 22t< %2t¢ 


99t2> 22tc 


Re 
Whole area in the first quadrant en- 


closed by the le curve. 


Area nearer to the origin in the 
first quadrant enclosed by the [z 
(The area in the closed loop is 
non-Re ) 


Whole area in the first quadrant en- 


closed by the (e < 


Area in the first quadrant enclosed 
by the closed loop encircling in the 
clockwise direction. 
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3-4. Application of fourth order charts. 


laa} 
( 


Since the curves constructed in the previous sect! 


fourth order system and the nature of the curves analyzes 2222 ah. 
versal to the fourth order systems, they can be used Im both analyeic 
and design of fourth order systems. It must be moted that fo2 cucyes ar 


plotted for some pertinent values of the parameters @o,y atc So om whi 
the curves depend and there must be employed some finterpolation tesla 
to apply to a specific system. 

In this section, the fourth order charts are applied ts: *°me a mer 
cal examples. The procedures and methods of analyzing or designing « 
system are the same as developed in References 1 and 2. 

Example 3.1 
Given a fourth order system with open loop transfer furc*ion 


Ce?) 4 


where K, Z and P are variable parameters. 
The problem is to analyze the effects of varying K, Z and P on the 
system stability. 


The closed loop transfer function is - 


1c is ) = i, ee ee eon 
R S* + (@+P)g3+(2+2pP)S*+@P+K)S + KZ 


and transforms into, by using a transformation S = (2+ P)s, 


fK Kz 
(v= @rpys? + Grp Rees. a 
A +6) + 2 Kz ’ 


aS 
R ._ 
(2+PPp (2+p)3 (2p y4 





thus the transformed coefficients are 
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Kz 
oe = FP)™ 


2pPeri< 
2g 2 


a 2a lid 
ae 
K 
bit 7 (2 & oe 
The transformed characteristic equation has unity coefficients of 


4 3 
the terms s ands and the charts can be applied. 


Suppose K and Z are fixed. By varying P, A, 9 aye and Age ate 


varied indicating that both the curves and the M poimt move. First com- 
sider the change of the curves. Noting that the shape of the curves 
(Family oe ) depends only upon ap,, the change of the curves 


can be analyzed by investigating the relation between ay, and P. 


2t 
Assuming P is any value greater than zero, from Equation (3-38) a», 
is a function of P. To see the dependence of ay, on F, ay, is diffsren- 
tiated with respect to P 


dao, _ _2P 
Cae 


dP 


implies that ay, isa monotonically decreasing function of P. ay, maxi- 





<0 for Pro 


mum is obtained at P = 0 thus 
dae max =°? 
In view of Table 3.1, a5,, = -5 is between G= .7 and Gs .75 
(a5,, = -5 corresponds to S = .707). Noting that az, > a),, provides 
no region of Ry for G=.707, for very small values of P, there is no 


region for the M point to guarantee <4 => .707. 


As P increases from zero, aot decreases from .5 and approaches zero 


when P tends to infinity. Therefore if © 2 .707 is required then P must 
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be considerably larger than zero. For example if an over-damped tran- 

sient response is required, then Arte = »375 indicates that P> 1.995. 
The locus of the M point with respect to the le curves must be 

analyzed in order to see the system relative stability. The locus of 


the M points is defined by, from Equation (3-38), 


2PtK 
Pes (eee 


We 
Sea SE 
Suppose P is fixed, then a5, is fixed to give a fixed set of Is . 
The locus of the M point can be analyzed on the Bo, vs By, plane where 


the family en is plotted. 


For the problem to be specific, suppose P is fixed to a number very 
close to zero. Then Bop is very close to .5 and the locus of the M point 
for varying K and/or Z is analyzed on the chart for ay, = -5. For this 


case, the locus of the M point is defined by 


ae 


K2 
B= a7 


(3-40) 


which is obtained from Equation (3-39) by neglecting P. 

Equation (3-40) is analyzed as follows. 

If K is treated as a parameter then each value of K determines the 
abscissa of the M point according to the first of Equation (3-49) and the 
ordinate of the M point varies linearly with Z according to the second 
of Equation (3-40), thus the M point moves along a vertical line. Thus 
for a set of values of K, there results a family of vertical lines along 


which the M point moves. 


180 







—~ * a  e 
+ ee Cee @ 
-— <a 

(| — = ¢ ee 





If Z is treated as a parameter then by eliminating K from Equation 
(3-40), the locus of the M point becomes 
By. =-52 Bit (3°41) 
which is a family of lines with slopes of .5Z passing through the origin. 
Thus for each value of Z, the M point moves along a iime defined by Equa- 
tion (3-41) as K varies. In Fig. 3-45 which is a duplicate of the chart 


| 


Bh aor for ao, = .5 the loci of the M points (the family of the lines 


mentioned above) are shown. If K = 1.0 then the M point moves along the 
vertical line for K = 1.0 as Z varies and if Z = .10 then the M point 
moves along the slant line for Z = .10 as K varies, thus for each walue 
of K (or Z), there is a limiting value of Z (or K) for a specified system 
stability. For example if K = 1.0 then for the system to have 5 7 .5 
the value of Z must be 21 < Z2< Zo where 2) is the lower bound of Z for 
the M point to stay inside the region Ris and is found at the lower inter- 
section point of the vertical line for K = 1.0 and the T's curve and 
analogously for Z5- In like manner for a fixed value of Z, say Z = .10, 
K must be less than 1.5 in order for the system to have a stability 
greater than = .3 since the intersection point of the lines for 
K = 1.0 and Z = .10 is about on the curve. Also it is seen from Fig. 
3-45 that if K is greater than 1.5 then it would be very hard to make 
the system have $2.3 and if K > 1.52 then there is no value of Z to 
make the system have 6 > .3 because the maximum abscissa of the Ge 
curve is about .19 (K = (8) (.19) = 1.52) 

Thus for P = 0 the system stability for varying K and/or Z is ana- 
lyzed. For any value of P the analysis can proceed in the same way. 


Now suppose the system given by Eyuation (3-36) is obtained from a 
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if 











third order system by cascading a filter network, that is, given @ third 
order system open loop transfer function 


IK 


/s Peal ar pve = Se ie ae a (%uoAey 
OS) rea i 
Cascade a compensator represented by 
+ 14° 
G.(6) = K(S+z) (3-43) 


Se 
If K = 4 then the system (3-42) is at the stability limit which can be 
readily seen from the third order chart and has the velocity constant 
Ky = 2, 

Consider the problem of designing a lag compensator which is repre- 
sented by Equation (3-43) so that the system stability is greater than in 
terms of ¢ , $= -9, with minimum allowable velocity comstamt Ky = 4. 

The problem is to determine K, Z and P in Equation (3-36) g0 that 
the specifications are met. By following the analysis made above and 
from Fig. 3-45, K=1.0, Z = .l and P= 0 will give the solution. The 
value of P can be determined from the Ky restriction or im a lazy compen- 
Ssator P = .1Z (which is a reasonable pole-zero ratio as illustrated in 
Reference 1) gives P = .01 and this P provides K, = 5 which is acceptable 
in this problem, thus the compensated system is 

; : en 
GE) ~ epee oe 
Example 3.2 
Consider a unity feedback servo system with its transfer function 


[000 
eo So (caae ae 


and it is desired to design a lag-lead compensator in order for the sys- 


tem to meet the following performance requirement: 


183 





1) The velocity constant K, is to be maintained at 100. 
2) The settling time is to be less than .5 sec. 
3) The system is to have a © greater than .7. 
By cascading a lag-lead compensator to the original system the com- 


pensated system becomes of fourth order and the transfer function is 


: K (S+2,)(St22) 
G(S)= S(S+l0)(S+P (57h) 


where Z1> Py and Zo, Po are for the lag and lead sections 
respectively of the compensator. 
There are five parameters K, 2, Py> Zo and Py to be determined. 
K can be determined from the Ky specification by choosing 4 suitable at- 
tenuation factor of the compensator and P, the pole of the lag section 
can be neglected in Mitrovic's method without danger, thus there reamin 
three parameters Z,, Zp and Py to be determined. (This kind of problem 
is solved in Reference 1 and in this example it is to show how the fourth 
order charts can be applied.) 
If the compensator is chosen such that the zero to pole ratio of the 
lag section is to be the reciprocal of the zero to pole ratio of the lead 
section, then K must be set to 1000, thus the transfer function which is 


to be considered in the Mitrovic's method is 


_1000(St2 )(St22) 
S*($+10)($ +R) 


to give the characteristic equation 


G-6)= 


B= O° +(04B)S?-/000+/0R ieee 1000(2,+2,)S + |0002,2,=0 


with the restriction on zeros and poles of the compensator as 


thus there remain three parameters to be determined. 
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There seems to be two ways of applying the fourth order charts to 
solve the problem: (1) by choosing a suitable Po, the coefficients a, , 
a3 and a, of the characteristic equation are determined to give a fixed 
set of By vs B, curves (for various values of 5 ) and then one consults 
with the charts for the families ee (Fig. 3-3 to Fig. 3-22)3 and 
(2) without choosing Py first then the coefficients ay amd ay of the 
characteristic equation are variable and one consults with the charts 
for the families levcana (Fig. 3-31 to Fig. 3-44). in this example 
the problem is solved by the first method and in the next example the 
second method is shown. 

In order to apply the fourth order chart, the system closed loop 
transfer function is transformed into, by using the transformation 
S = a,8 = (10+ Py)s , 


bed + bet + Aor 
eC)" ae ee 


ss JOoO+ OR _ 1000 
where Bor ios 26 ee = aie 
(3-44) 
(000(2,\tZ2) 


a= oF Ee 


_ (00022, 
ae Oe 


From the design specification (3), the M point must be inside a Koy 
region. If one chooses $=: .8 (Rg is inside R 7 for a fixed a», if 
R.g exists) then from the design specification (2) @, the frequency at 
the M point on the I’: curve must be greater than 10 assumimg the set= 
tling time is approximately 4/6 @),,, and this frequency is, in the 
transformed coordinates system, 

Wnt = p/a, =10/(10+ FB) 
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In view of Table 3.1 for the existences of the R, regions the con- 


Ss 


ditions on ao, are 


R.q 4a4< .505 
R.s Aor 44S 
R.9 Ae <. YOU 
Rio Agee. 372 


For the existence of Rg, the equation for ay in Equation (3-44) and 


t 
the above condition on ay, requires Py > 49, 
If one chooses ay, = .3 which is a sufficiently small number so that 
it provides a good flexibility to locate the M point to guarantee 2.7 
then Py = 63.85, thus one obtains the following data 
P=63.e5 
4,= 73.85 
4e=-3 
Oe?-137 ona [7g curve 
At this point it is checked if the value of Py determined on the 
basis of the region Re is sufficiently large so that the larger of the 
newly introduced roots of the system is much greater than the real parts 
of the dominant conjugate roots, 
Next one consults with the family Em for a>, = «3 which 
is in Fig. 3-8. Fig. 3-46 is a duplicate of Fig. 3-8. The curves for 
G = .7, .8, .9 and 1.0 with frequencies on each curve are available. 
It is seen on the chart that the line By; = 4 X 1072 is very close to the 
[78 curve for ,150 <Wie < .200 which is an allowable range of W,, at 


the M point of this problem. Thus the abscissa of the M point is set to 
SS 
a,= 4X10 
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This a;,_ and Py substituted into Equation (3-44) gives 


Let Z5 = -1P, then Z5 = 6.385 and Z, = 9.615 to give 


=61.4 
From Equation (3-44) aoe = 20.9 X 1074 thus the M point is located at 
(a1, > 4,) = (4 X 10-2, 20.9 x 1074) and this point is at ®,, = .150 
on the I's curve which is a satisfactory M point for the problem. 
Since Zo = me Po is chosen, then P,; = .1 Z, must be chosen in view 
of the chosen compensator pole-zero ratios, thus Py = .9615 and the com- 
pensated system transfer function is 


_ 1000 (S +9.615)(S+ 6.385) 
GS) = geno) (Se.96IS (CEES) 


If one chose aj; = 4.25 X 1072 then using the proportionality rela- 
tion between a), and (Z,+ 29) as seen from Equation (3-44), (24+ Z9) 


associated with Els 4.25 x 1072 is 


Z,+2,=(16)(23"0 ay yas 


2517-2 217- 6.395 =10.62 
Z,2,= 68 


again using the proportionality relation between ane and 21,25 


-&% 
4, =(20.9x10 (SF) =23,!1x10 
and so determined M point denoted by My in Fig. 3-46 is also satisfactory, 


then the compensated system transfer function is 


(So a 100O Geach ages He S32) 
~"S(S+10)(S + 1.062)(8'4+63,85) 


Thus if a first try is not satisfactory then a second try can be made 
without too much calculation and the direction whether to increase or de- 
crease in aj, in the second try can be determined by inspecting the re-~- 


sults of the first try. 
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Example 3.3 

Consider the same problem of Example 3.2 with the ¢ specification 
as C = 1.0 with the settling time specification omitted. The problem 
can be solved on the charts of the families le tad ° 

By choosing some suitable pole-zero ratios of the compensator, Z 
and Za can be expressed in terms of Py which is the only parameter to be 
determined. 


By choosing a compensator such that 


Zo= 1B 
iA =.00| Ie 
then from Equation (3-44), a » a,, and a which must be determined 


2t lt Oc 


are expressed in terms of Py only, namely 


_ (000 +tloR a 
a4 = (Meet ae 


HIOR, 
ae (+B tad 


Ql 
oe (10+ BY 


The latter two of Equation (3-45) define the locus of the M points 
on the Boe vs Bit plane if a set of continuous values of Po are supplied, 
and this curve represents all locations of the M points (among which a 
suitable M point must be selected) associated with the compensator whose 
pole-zero ratios are set as above. 

As Po varies, the location of the M point as well as the \e curve 
on the 6 = constant plane varies, Since on a 6 = constant plane the 


variation of the le curves can be analyzed by the constant QW); lines 
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and since the location of the M point moves along the curve defined by 

the latter two of Equation (3-45), one can find a value of. Po = Po such 

that at Ba the M point is on the le curve for a5, determined by P5 
To plot the locus of the M points defined by Equation (3-45), the 


following table is constructed by remembering that for 4 = 1.0, 
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The locus of the M point is plotted in Fig. 3-47 from the table 
above. Fig. 3-47 is the chart of the family bas for ao; = .2 
(It is a duplicate of Fig. 3-44). The ie curves for a9; = .2, .18, 
-16, . . . are as shown. The frequencies on the ee curves for 
age = +2, .18 and .16 are labeled whereas the frequencies on the ly 
curves for a5; € .14 are not shown. However the frequencies on any curve 
can be determined by using the constant ©, lines; for example, if one 
constructs a constant @,, line for oer .06 by drawing a line which 
passes through the two points each at G),; = .06 on the ans curves for 


age = -16 and ag_ = .18 then the line crosses the on curve for 


a2?t = .14 at the point marked by X sign and this point is of Wy; = .06 
on the re. curve for aj¢ = 14. 
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On the chart with the locus of the M points plotted as shown, the 
relation between the family |g Saye} and the locus of the M points is 
analyzed. On each point on the locus of the M points there is associated 
a value of Py which determines a9,- For example if Py = 80 then 
a9, = -222. Noting that the Lats curve for aj; = .222 will be shifted 
to the right of the curve for ao, = .2, if Po = 80 is chosen cneneete M 
point will be inside the iw curve for a9; = .222 implying that Po = 80 
makes the system over-damped. Consider the M point for Py = 120. For 
Po = 120 ay = .13 and the M point will be outside of the [ae curve for 
@9, = .13 since in Fig. 3-47, the Fe curve for a>, = .14 is very close 
to the M point for Po = 120 and it can be estimated that the Fe. curve 
for a5, = .13 will be between the ee curves for apy = .12 and a2¢ = .14, 
thus P, = 120 makes the system under-damped with a fairly high &. Con- 
sider the M point for Py = 100 which gives ay, = .165. Let the point QQ 
and Q% be the points at which the locus of the M points intersects with 
the ae curves for a5, = 16 and Gop F -18 respectively. a5, of the 
Lg curves varies from .16 to .18 between the points Q; and Qj). Noting 
that the portion of the locus of the M points between Q; and Q) is nearly 
a straight line, and that the M point for Py = 100 is at about 1/4 way 
from Q; to Qo, by using a linear interpolation, one can predict that the 


I1.0 curve for aj; = .165 will pass through the point M (a; , ao) which 


implies that if Pp = 100 is chosen then the M point will be on the lo 
curve, thus Py = 100 gives the solution. 

The frequency at the point M (a), ; aoe) can be determined with the 
constant 4); line for Ont = .09 as drawn in the figure which shows 


that Woe (or c ) = .09. From the relation between transformed and 


192 





untransformed system frequencies one can obtain 

=A, S=(10+ Pe) <4 =(I1D)09) = 9.9 
which is the value of the double real roots of the compensated system, 
and the compensated transfer function is 


S _jisoo (satatsehiesenlo ju 
G(S)= esos +1 CS F100) 


which happens to be a cancellation compensation. 

It must be noted that one can construct the ear curve for ay, = 
-165 (which is for Py = 100) by using the constant Wat lines to notice 
that the point M (a,; » 4p,) in Fig. 3-47 is slightly inside of the eg 


curve for a, 165. The effect of increasing Py on the system % 


t=: 
(against G = 1.0) can be seen from the figure, that is, by increasing 
Py from 70, the system 6 is decreased from $71 and it assumes 

$= 1 at P, = 100, The same analysis can be made for Py < 70 if the 
locus of the M points is plotted on the o = 1.0 plane for that range of 


Po - 
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fy Design. 

In general, design or compensation of a given system is to reshape 
the system equation by using some compensators so that the system roots 
are placed at suitable locations on the S-plane with due considerations 
of the specified performance requirements. 

In the usual design, the steady state characteristics are specified 
by designating the minimum allowable steady state error which can be 
interpreted in terms of the velocity error constant Ky (or K, » Kp etc. 
depending upon the system type number) and the transient response char- 
acteristics can be interpreted in terms of the system dominant roots 
locations. 

In this section, assuming that the design specifications can be 
interpreted in terms of the steady state error constants and the loca- 
tions of the dominant roots on the S-plane, an analytical as well as 
graphical method of designing cascaded compensators is developed by using 
fundamental properties of the basic Mitrovic's equations aaa curves in 
order to minimize the amount of trial and error which is usually encoun~ 
tered in design problems covered in the literature. Also it turns out 
that the labor can be considerably economized and the difficulty of 
choosing some suitable values of poles of the compensators (which is 


essential in the current Mitrovic's method) is somewhat removed. 
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4-1, Design of single section compensator. 


Consider an open loop transfer function 


Ges nae (4-1) 


where Ky is the minimum required forward gain and D(S) is a 
polynomial of S with all positive coefficients, 
and the design specifications require that the system dominant roots be 
located at 
S=-$,On1= [\-8? Wnt 
without hurting steady state performance which is specified by the velo- 
city constant Ky (or K,, Kp etc. according to the system type number). 
Since the gain is adjusted to K, the Mitrovic's working point 
M(a, P ag) associated with Equation (4-1) insures the specified velocity 
constant K, but in general it does not give the roots at the specified 
locations. In order for the system to meet both K, specification and 
the desired roots locations the point M(a, , ay) which guarantees the 
specified velocity constant must be at Wy, = @,) on the e 1 curve, 
and moreover the point at a = @y, on the 51 curve must be such a2 
point that it guarantees the system 6 not less than Sy 
Let the By vs By plane be the one on which the system of interest 
is analyzed, and let the Mitrovic's variables for the original system in 
Equation (4-1) be BY and Bo with the superscript '"o" to denote the 
original system. Then the Mitrovic's curve for $= o of the original 


system is 
B= Bi ($4, Wn ) 


ee. (4-2) 
B =B ($1, ®n) 
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by choosing the coordinate axes such that By axis == B, axis and BO 
axis == By axis, the curve in Equation (4-2) can be constructed on the 
Bg vs By plane and this is denoted by lz, in Fig. 4-1 with the super- 
script "o" to indicate for the original system. 


ih 


Let the coefficient of the S~ term in the polynomial D(S) be d 


1 
assuming a type one system in Equation (4-1). Then the pair (d, : Ko) 
determines the location of the Mitrovic's working point on the Bo vs By 
plane and it is located at 


M= M Cds > Ke) 


as shown in Fig. 4-1. The radius vector OM has a fixed slope 


Ko 
dt 
and this vector is denoted by Ly, : 


= Ky 


The point at GW, = @_, on the SF curve is located at the point 
denoted by MG which in general does not coincide with M° requiring a com- 
pensation. Thus the point Mg is the required location of the Mitrovic's 


Oo 
working point on the [< 


, curve in order to meet the roots specification 


and the coordinates of Ma on the By vs By plane are, by letting 


a. = Wn in Equation (4-2), 


Ma=Mi ( Bs, » Wnt), B5($.,@ns)) 


The radius vector Omg has a slope BG (5) ‘ Wa! Br ($4>@y4q)> 
° 
The he curve, M° and Mg points configuration depends entirely upon 
the original system and the design specifications. 
The ultimate goal is somehow to make the points M° and Ma coincide 
ce] 


since if mM? = Ma 


Since there is no way to move M° and M3 in Equation (4-1), a filter 


then both K, and the roots specifications are met. 
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L (The locus of Mg) 





> (Rane = ant as 





BS 





: \ 
The locus of the point MG ata, =n) 


on the Bo vs, By curve for ¥= 4) 


of the compensated system. 


197 





network single section compensator is cascaded to the original system to 


give a compensated system equation 


= Sise2) 
GS) = Bess?) a0 


where Z and P are zero and pole of the compensator and are 
variable parameters. K is the new gain which has to be 
adjusted to meet the steady state performance requirement. 

Let & be a number such that Z =P, then from the Ky requirement 


K must be adjusted such that 
Ne a or k= Krai ° 


Substitution of this K into Equation (4-3) gives the compensated system 


equation as 


Gs)= KoA (S+PA) 


DiS)(S+P) Set? 


where Z is expressed in terms of P andd . 

The compensated system in Equation (4-4) has its own Mitrovic's 
curves and Mitrovic's working point. By cascading the filter network, 
the significant points and curves for the original system on the Bg vs 
B, plane are changed into those of the compensated system. Those changes 
can be analyzed mathematically as well as graphically. 

Let Mitrovic's curves and the points (corresponding to M°? and MG of 
the original system) be superscripted by "c'' for the compensated system. 


The characteristic equation of the compensated system is 
D(S)(S+P)+ Kor (S+P/L) =O 


or SD6) + PDS) + MAS + KP =O (4-5) 
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It is seen from Equation (4-5) that any Bp vs By, curve for the com- 
pensated system is composed of two component curves, one from the poly- 
nomial SD(S) of which all of the coefficients are fixed, thus providing 
a fixed curve, and the other from the polynomial PD(S) which is the 
original system polynomial, from which the original system curve is ob- 
tained, multiplied by P. Let By and Bo be the Mitrovic's variables for 
the polynomial SD(S) and let the coordinate axes be such that By axis = 


By axis = Bi axis and Bo axis = B5 axis = B, axis, then the equation 


0 


defining Bg vs B, curve for ‘5 O14 of the compensated system are 


B= B(8,, dn) +PBV(S,, Bn) 


? ° 4-6) 
Po =B. (> On) +PB,(5)> @,) ‘ ) 


The interpretation of Equation (4-6) on the By vs B, plane is as fol- 
lows. First of all the Bp vs By curve for C= O1 of the compensated 


A vs By curve for 4% = 54 on the Bo vs By plane, 


system, namely the B 
depends only upon P. Assume the Bo vs By curve for 1 denoted by TZ, 
(which is obtained from the polynomial SD(S)) be constructed as shown in 
Fig. 4-1. On the two curves ie and ie, » the points at any fre- 
quency W, = we can be located. If the coordinates of the point at 


O - a! 


a 
* n on the Ie, curve are multiplied by P, and added to the co- 


ee 


ordinates of the point at WJ, = @,, on the S$, curve componentwise , 


c 
then there locates the point at @),, = an on the le curve. By run- 
¢ 
ning aes from 0 to Om , the ce curve for the frequency range from 0 
to On can be constructed which is not shown in Fig. 4-1 since it is not 


c 
necessary to construct the le curve at this point; instead one focuses 


c 
attention to the particular point Ms which is at oO, = Wi on the ie 
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4 3 “ c 
curve. Mg is a moving point as P varies since Mg is at G, = aJ,4 on the 


n 


(5 c 
le, curve and the le curve moves as P varies whereas the point Mg 
(tex 
is a fixed point at = Mp4 on the 3, curve. Thus the point Mg 
is not indicated at a particular location on the Bg vs By plane in Fig. 
4-1 but the movement of MG affected by P is investigated. 
By fixing the frequency at , = @,4 im Equation (4-6) the coordin- 


n 


ate of MG on the By vs B, plane is determined by 


By =F) (5, On) +P Be (S, Ont) 
° (4-7) 
Bo = Bil Sia Wns) +P BS, @n) 


which is the equation defining the locus of the MG points on the By vs 
2 t ' 
By plane as P varies. Noting that B) (Sy: Wn)? BO (S$); G4)» 
o o 3 
By ¢ $1: Wa4? and Bo C5 On) are fixed numbers, Equation (4-7) defines 
a line on the Bg vs By plane such that it passes through the point 
' i] 
(By (5,,@.4); By ¢ 61> 0,4) denoted by Q which is the point at 
a ; ° ° 
Ges td) ,4 on the le, curve with a ge: Bo ( Si? 0) 4/3, ( oir a4) 
which is the slope of the radius vector OM - This can be visualized by 
eliminating P in Equation (4-7) to obtain 
’ ; : 
Bo - By ($9 @n2) = 2-2 —™ 1B, - Be, Wa) 
; : es, > tt) = a 


thus the locus of M, points defined by Equation (4-7) is a line parallel 


Be nnn (B 


to OMG and passes through the point Q (By CS a4)» Bo ( oe w 4) 
on the By vs B, plane. This is denoted by L in Fig. 4-1, At each point 
on the line L, there is associated a value of P according to Equation 
(4-7) and the direction of increasing P is in the direction of the radius 


—_— 


Oo 
vector OM, . At P = 0, the point at (ys W nd OP the compensated Bo 
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c ui ; pete 
vs B, curve for = Sy (Mg on i ) is at Q which is conformed fro 
«| 


Equation (4-7). As P increases from zero, My moves along the line L ii 
the direction determined by om’ . Thus the desired location cf the Mic- 
rovic's working point (to insure the roots specification) which was orig: 
inally at Ma jumps to the point Q upon cascading a single section filter 
network with the pole at the origin and moves along the line L as the 


pole of the filter network increases such that the distunce from © to it 
Le) 


is proportional to Pp, If a filter network with finite ,»ole is cascadd, 


oO 


then it jumps from the point Mj 


to a point on the line L. 

It is readily seen that the extension of the lin: L i+ the ogrocite 
direction of the vector on beyond the point G corresporoes to the néega- 
tive P which means the pole of the compensator on the ripht half of the 
S-plane. 

The orientation of the line L depends upon the locations of the potat 
Mg and Q which both may be in any quadrant on the Bg vs B, plane depend- 
ing upon the original system and the design specification. 

So far the effects of cascading a single section filter network on 
the Ma point are analyzed. Remembering that in general the Mitrovic's 
working point determined by the pair (a, 5 ay) the coefficients of the 
sl and 3° terms respectively of the characteristic equation is not 
identical to the desired M point located at 3 Beat the | , 
thus it is necessary to adjust the variable parameters t» bring the two 
points together to a suitable location. To do this the behavior of the 
point determined by the coefficients ay and ag must be investigated, 

The point vee is at (d, ; Ky) for the original system. By cascading 


the filter, this point is moved to some other location on the Bg vs By 
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plane since the coefficients of the characteristic equation are changed. 
Let the point M be the point for the compensated system analogous to 
the point M° of the original system. Then from Equation (4-5) the char- 
acteristic equation of the compensated system, the M© point is defined by 
M=aM (dP + hod , Ko P) 

which implies that the location of the M© point on the Bo vs B, plane de- 
pends upon both & and P. 

The locus of the M points on the Bo vs B, plane for dX and P vari- 
ables is defined by 


8, = d,P + Kok 


(4-8) 
ey = Kop 

By eliminating P in Equation (4-8) and noting that = K 94 » one obtains 
Boe Ky (Ban Kod) (4-9) 


which is an equation of lines with the slopes fixed at K, and the B,-axis 
intersections given by K,A which is a linear function of A , thus the 
loci of the M° rotnts are a family of lines parallel to the line Dee : 
In Fig. 4-2, some members of this family are shown with the configuration 
of the line L, the vectors om? , ong and Q point assumed to be the same 
as in Fig. 4-1. The line Ley (ome vector) itself is a member of the 
family defined by Equation (4-8) or Equation (4-9) corresponding to 

cA : 0. Any member of the family of the loci of the M° points can be ob- 
tained by translating the line Ley horizontally since the slope of any 
member is fixed at K,, and for any two members one for d= ay and the 
other for od: Ay » any two points one on each member read the same 


value of P if the two points are at the same ordinate, which is readily 


seen from the second of Equation (4-8). 
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and the locus of the M§ points. 
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Let the line L! be a member of the family. At each point on the 
line L} there associates a value of P linearly according to Equation 
(4-8). The point for P = 0 is on the By axis. If L, is divided inte 
two parts by the B, axis above and below, then the upper part is for 
P>O and the lower part is for P<0O. Now let the line LL move horizon- 


tally by varying <. Then for “>0_ L!, is on the right of 1 and for 


! 
“Kv 
A~ QO it is on the left of L' , thus the line L! is well defined for 
~Or~ c= ce 

From the definition of the point Mi» and the locus of MG points L 
if the Mitrovic's working point M is chosen at any point on the Line £ 
then the compensated system has roots at § = - aya as) W) at 
and from the equations defining the line iy , it is obvious that if the 
Mitrovic's working point M is chosen at any point on Li then the com-~ 
pensated system satisfies the K, requirement. Therefore in order for the 
compensated system to meet both the Ky and roots specifications, the Mit- 
rovic's working point M must be chosen such that the point is on both 
lines L and Ly and this implies that the Mitrovic's working point M 
must be at the intersection point of the two lines. If one is free from 
the parameter P appearing at each point on both line Li and L, and also 
free from any restriction on & the intersection point will be Located 


anywhere on the line L if the two lines are not parallel, since A can 


be varied to any number, and each intersection point determines 4 re- 
sulting infinitely many numbers of ds. However the lines Te and L both 


bear P at each point on the lines according to the linear functions of P 
defined in Equations (4-7) and (4-8). The proper point must be ches:n 
among the infinitely many intersection points (as mentioned above} sc 


204 





that at the intersection point the values of P on Ty and L must be the 
same values. This is mathematically speaking the two equations, Equa- 
tions (4-7) and (4-8), must be solved simultaneously. From the property 
of the linearity of Equations (4-7) and (4-8), there always exists one 
and only one solution point. Even in the case that the lines Te and L 
are parallel to each other, there exists a solution. Let the solution 
point be Ms, then the point M, is determined by solving Equations (4-7) 


and (4-8) simultaneously, namely 


Be (S$, War) 
hE oor 


(4-10) 
ol, = Al Bi 5, » Opn) ale ( Bis, ; Wns)—d)| 


where P, and are the value of P and & at the solution 


point M, 
By substituting P, and o, back into any of Equation (4-7) or (4-8) the 


location of the point M, on the Bg vs By plane is 


= 5 Bo Br CS, Wn) Ke B, ( Fi > Wr) 9 
M. M. (B ( sO te Ky = BS, y a Bo(%,, War) (4-11) 


=M.(di2 + Koos , Koi) 
where B, B°(S, Wp =Be(S,,n)-B, CS. On) 


So far, any restrictions to be imposed on P, <A and the location of 
M, in view of practical situations have not been considered, and P_ and 
A, in Equation (4-10) can be any number from negative infinity to posi- 
tive infinity and the location of M, defined by Equation (4-11) can be 


anywhere on the Bo vs B, plane. 
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Now consider the restrictions on P. and A; . In view of the fact 
that the filter network to be realized by a passive RC network, the 
first restrictions on P_ and «, are 

me 

x 70 
By considering a practical RC network either too small or too large « 
must be avoided. Thus one may restrict ~ such that 

NC d S10 
which is a reasonable restriction. 

In order for the compensated system to be a stable one, the Mitro- 
vic's working point must be in the first quadrant, thus the M. point 
must be in the first quadrant. If the solutions given by Equations 
(4-10) and (4-11) are within the restrictions made above, then the solu- 
tions may be the proper ones. If the solutions violate any of the re- 
strictions above, then this is an indication that the original system 
cannot be compensated to meet the design specifications by a practical 
single section filter network. Thus one must try to cascade another 
section of filter network or use some other scheme for compensation. 

The restrictions on the solutions can be interpreted on the Bo vs 
By plane. The restriction that requires the M. point to be in the first 


quadrant is fulfilled if 


because for that P, and A. , any point on Lh is in the first quadrant. 
The restriction .1<%< 10 indicates that the point M, must be on the 


area bounded by the Lines L', and Lio.g Which are the members of the 
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family of the loci of the M© points for A= .1 and «= 10 respectively, 
and those lines are shown in Fig. 4-2. Let Q' and Q" be the points on 
L at which are crossed by Ly and Lio.o tespectively. Then noting that 
the M. point must be on the line L, the allowed locations of the M. 
points are between Q’ and Q" on L. The line L', (or 110.0 °F both) may 
cross LL at a point below the By axis. If this is the case, then the 
allowed locations of the M, point must be governed by the restriction 
pe 0. 

The existence of a solution point M, between Q' and Q'' on the line 
L can be determined as follows. Let the point Q' read P = P and P = EH 
on the lines L and Ly respectively and let the point Q' read P = P,, and 
Ps oe on the lines L and Lino respectively, then there exists the point 
M between Q' and Q'' if and only if 

R>e and = P <P’ 

or noe and PB. >P 
This is an immediate consequence of the linear properties of the equations 
defining Ds and L, thus one es, check whether there is a point M, in 
the allowable region graphically; however it is much easier to check it 
by solving Equations (4-7) and (4-8) simultaneously to get the solutions 
in Equations (4-10) and (4-11), since the equations to be solved are 
simple linear equations. 

When the solutions turn out not to be acceptable, then by choosing 


some suitable o& and P, say Az A, » P=: Py , to get an intermediately 


compensated system 


G,(s)= Ky (sz: ) 


4-12 
D(S\(S+P) “ia 


where Ky is the new gain adjusted to give the specified K, 
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one can follow the same procedures as developed above for the system in 
Equation (4-12), that is the system in Equation (4-12) is treated as the 
system in Equation (4-1) and the polynomial D(S)(S+ Pi) in Equation (4-12) 
plays the role of the polynomial D(S) in Equation (4-1). By cascading 
another section of filter network to the system in Equation(4-12), the 


system of interest is then 


o: iste (S42) 
(S\= — oS YS 4-13 
GUIS)= Bey (StR Sep) aa 

The system in Eyuation (4-13) is essentially in the same form as 
the system in Equation (4-3) except that K which must be determined 


enters into the coefficient of s* 


term (affecting Bg vs By curves) of 
the characteristic equation for the system in Equation (4-13). The de- 
tails for this case including how to choose Zy and Py in the first stage 
will be discussed later. 

If the solutions are acceptable then there still remains one and 
the final check, "Are the roots at S = -$,Wy +/1-~%2 Wp; dominant?" 
and this is done by sketching the 5 curve for the compensated system 
(final) to see if the point M is on the region Res 


The examples shown below illustrate the arguments developed above 


for the single section cascaded compensator. 
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Example 4.1 
Consider the transfer function 
ae 
CC) S(S+5) 
which has a pair of complex conjugate roots at § = -2.5 9.2... Le is 
to be compensated, without reducing the velocity constant Ky = 20, so 


that the system has a pair of complex conjugate roots at 


= - 3, Vn t/i- 37 Oni =-10.5 eee 
which is equivalent to Sy = ./ and On = 15 


D(S) = S@ 


+ 58 and dy = 5, Ky = 100. The M° point is at (5, 100) and 
the line joining the origin and the M° point has a slope equal to Ky ceed 
which is the slope of the family of the lines L’ (the loci of the Mo 
point). 

The compensated system equation is, noting that d Ky = 100, by using 


Equation (4-4) 


100L( St PAK) 
(S7+5S)(S+P) 


G(s)= 
The Bg vs B, curves for the original system are defined by 
Br=B(S)W, 
Bo =On 
to give the location of the Mj point at 
=M, (B Be ae pe LS eh 
a Ja| ($=.7,O%=15), ($=. 4 =>) = My il 225) 
thus the MG point is located in the first quadrant and the radius vector 


OMG gives the direction of the Line L which is determined next. Note 
that the points M° and MG are not identical and a compensation is needed. 


The characteristic equation of the compensated system is 


(s3+ 587) + p(s?+ 5s) + 100%S+ LOOP = 0 
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in which SD(S) = sot 582 and PD(S) = P(S*t 5S). The Q point is deter- 
mined from the Bg vs By curve for the polynomial SD(S) evaluated at 
m =../ and Gi) = Isp) theateis 

Bes (8) adn + De C07 


1 tegen : 3 
ee =5 i), = G2($)W,, 


to give B'($=.7, Wy 215) = -1lt 
CS v Q),, = ils — 2600 


therefore the Q point is located at (-111, -3600) which is in the third 
quadrant. The locus of the Mg points, the line L is, by using Equation 
(4-7), 


B=-lHi+pl2n 


(4-14) 
By = -3600 +P(225) 


which define the line passing through the point Q (-111, -3600) with slope 


of 225/21 = 10.7. Note the slope 10.7 is the slope of the radius vector 


ee 


OMe . 
The locus of the point mM is, from the characteristic equation of 

the compensated system 

B, =5 P+ |00K 

By =100P (4-15) 
At this point Equations (4-14) and (4-15) may be solved simultaneously 
to get the solutions of P and < without constructing any curves or 
lines; however for the illustrative purpose some significant lines are 
drawn in Fig. 4-3. Noting that the points vee 3 Mg and Q are needed only 
for determining the slopes and the direction etc. of the lines, those 
points are not shown in Fig. 4-3 but the necessary portions of the lines 


are drawn. 
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The Line L is constructed as follows, From Equation (4-14) the By 
axis intersection is, by letting Bp = 0 to obtain nS 16 and by substi- 
tuting P, into the first of Equation (4-14), 

Bi intersect.= 225 at Bae 
thus the point y is at (225, 0). From Qy the line L is directed into 


the first quadrant with slope of 10.7. 


Dae 
QQ with slope equal ti K, = 20 since from the first of Equation (4-15) 


The locus of M° point for & = 2.25 denoted by L 5 Passes through 


by letting B, = 225 and P = 0, is found to be 2.25. The line Lj 

the locus of M° for o&= 8 is drawn by noting in Equation (4-15) that 
the B, axis intersection of the line for A = 8 is 800 and Ls 25 and Le 
are parallel. Note that on the two lines Ly 95 and Eee for each pair 
of points at the same ordinates, one on faa and the other on Ly » the 


values of P are the same. At point Q% the lines Lg and L intersect and 


We ! 2 
read Py = 131.5 on Lg and P, = 


74.5 on L. At the point Q) Pi = O and 
Pi s 16, thus 

Ear 

Fe > Ke 


implies that a solution point M, exists between Q; and Q, on the line L. 
Thus the compensation can be accomplished by a single section lead filter 
network with 2.25<A< 8 which is an acceptable one. 

By solving Equations (4-14) and (4-15) simultaneously or by using 


Equation (4-10), one obtains 
3 = 28.8 


ky = 3.5 


to give the locations of the My and M points at M, (494 , 2880) 
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and Z.= & /os = 3.23 
ee wige k= 257 


then the compensated system is 


a oa Se tat a) 
GS) = S(St+5)(S +229) 


Finally the compensated system is checked on the third order chart 
or by constructing By vs B, curves for the system to see that it is ade- 
quately compensated. Note that the characteristic equation of the com- 
pensated system factors as 

57433.0S7+ GINS +270 =(S +215 +225)(S+12.9) 
to give the roots at 
= gi) 11 


S=-|2.8 


Zs 





Example 4.2 


Given an open lcop transfer function 


é 
GS)" e(e24 2.05412) 


which is analyzed on the third order charts to be an unstable system with 
a pair of complex conjugate roots in the right half plane, design a com- 
pensator to place a pair of complex conjugate roots at 5: -3 and 

QJ, = -6 with the velocity constant maintained at Ky, =: 5. 


For the original system the significant points are 


M’=M'(1.2, 6) 
Ma = Mal 8°(.3,.6), Bo C3,-6)) 
=M; (.95, 59) (4-16) 


The radius vector OMg has the slope = .59/.95 = .63 


By cascading a filter network the compensated system is 
6A (S+ FH) 
S(Ga1 2S tla) Csre) 


to give the characteristic equation 


Elis 


(eee |.287) + P(S3+ 2.874 1.28) 4+ 6kS + 6P=O 


from which the Bo vs By curves for = .3 and the Q point are obtained 


as 


B= .72Wny +128 Oe - 9P¢ a)? 


2 
ceeewe steero 6a 


BL (S$=.3,@=.6)=.65/ 
(4-17) 


B,($=-3, Oe-6)=.09 


thus the point Q is at (.681, .09) and the locus of Mg point L is 
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defined by, using the coordinates of Mg and Q in Equations (4-16) and 
(4-17), 
B,=-69! dS) 


Ge = .0oF oy PEAY ) (4-18) 


From the characteristic equation of the compensated system the coeffi- 


0 


cients of the S! and §° terms are seen to be 1.2P 4 6A and 6; respec- 


tively, thus the locus of the yO points is 

Bel 2Pproenr 

ene (4-19) 
By solving Equations (4-18) and (4-19) simultaneously or by using Equa- 
tions (4-10) and (4-11) with the proper substitutions of the coordinates 


° 
of the points My and Q from Equations (4-16) and (4-17) respectively, one 


obtains 


ee 


Fe =.0166 
t=. 112 Zt fly =-14PF fee hore ers 


and the point M, at (.705, .0996) which are the acceptable solutions. 


The compensated system is then 


WEBS te £8) 
S (67425 +1-2\(S+.0766) 


Gis)= 
and can be checked on the fourth order chart for the dominance of the 


roots at = .3 and @), = .6 


The closed loop transfer function is 


67 3(5 4 18) 


Cave Meee Reed See 
R (s)= 675 2.0/a6 S24 |.2892 5° 4-795 Ss Faye 


R223 ps. 
R 


Eye = 
a; 
= 


from which 25-0 Noting that ao, is less than ay 


tc 
for Gz .3 (see Table 3.1 in Section 3-3) and the point M, is in the 


first quadrant, the roots at “= .3 and @, = .6 are the dominant pair. 
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4-2. Design of double section compensators. 

In a design problem, one may find that a given system cannot be 
compensated to meet the specified performance requirements by a single 
section compensator, and this is the case that the set of solutions ob- 
tained by cascading a single section filter network as in Section 4-1 
turns out to be unacceptable and a multi-section compensator is needed. 

Let the set of solutions obtained by cascading a single section 
compensator be 

aie 

dk =o 4 and. 2= & /o, 
and suppose the solutions are not acceptable. This implies that the orig- 
inal system must be cascaded with a compensator with more than one sec- 
tion, thus the designer may try to cascade a double section compensator 
and proceed with exactly the same approach as in the case of a single 
section compensation. 

The unacceptable solutions from the single section compensation are 
first analyzed before cascading a double section compensator. The unac- 


ceptable solutions fit into one of the following cases: 


fy ita | >] 
la) 70, ol, > 0 
Se, u& <0 
1c) B <1), A, >O 
id) Fe <0, ee 
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2) with lay] <2 


2a) ry 20 O<dec.l 


2c ier 10 A, >0 
2d) iF <0 2 <O 


The cases b, c, d are not acceptable because of the inrealizability of 
the compensator with a passive network. The case la is for the noise 
problem and the restriction %< 10 is not a theoretical one. In case 
2a which is a lag network, if & is too small then the sizes of the com- 
ponents of the filter network may become impractical, yet some & con- 


siderably smaller than .1 may be acceptable. 


> so 
Let G_(s)= Ke (S¥ % be a fictitious compensator 
vy ars 
a5 
where K, is the attenuation factor and 2, - Py Lok : Py and ky fit 


into one of the above cases which are the solutions obtained by cascading 
a single section compensator and Gig (8) is not acceptable. 

Te Gog (S) is cascaded to the original system G,(S) in Equation 
(4-1), then the transfer function 6, (SG. (S) satisfies the well-known 


roots locus criteria thus 


ox, (S+2,) = 


b ee (SP = =~ 4-20 
Ge(S) Gep (S) DEV(S+P) 1 (4-20) 


ac S=-5,On; + JS? Ont = 15) 


Noting that Kae 0, ty can be either positive or negative the 


angle criterion in Equation (4-20) for the designated root r, is 
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£.65(S) 





+ Ye eh 2 = 20) (4-21) 
5 


set, SHE her 


where n is an integer. 


Let f be the angle that must be added to the angle LG) ere by 
a 


some compensator in order for the compensated system to have roots at 


the designated Locations, then 


Lo)... + fee 7-!) i (4-22) 


From Equations (4-21) and (4-22) 


A (St%,) 





s=F 


Sez 
or Peek 2 if od, >o 


Se: Sar i 
(4-23) 
S +z, : 
Y= aap if Sone 


S=F; 

Since Oo fe 180° indicates that Gy(S) needs a lead angle in order for 
the compensated system to have the roots at the designated locations, 
and similarly since -180°< ¥ < 0° (or 180° < ‘f < 360°) indicates a lag 
angle, the unacceptable solutions for P and A tell the designer the 
type of compensator needed. 

With the aid of Fig. 4-4 one can obtain the following table which 
shows the proper type of compensator needed in view of the angle cri- 
terion for the possible combinations of the unacceptable solutions Ae 
and Pe. 
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jw : 
2a. P.> Oo 


la. Py>o 1 
ay >10 O<ol, <.1 


c 









ee. hia 0 
> 0 









ides “Py <0 
a) <0 





Fig. 4-4. The possible pole zerd configurations ' 
of the unacceptable single section 
compensator, (S+P;/e )/(S+P)). 
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5 i |compensator 








Positive 0°< ¥< 180° ; lead 
ee Mie 
Positive 180° £<360° , lag 




















A) >] 
| i, Cen ee ag 
| | Negative -180°<¥ < 0° lag 
tee eee: wee ces ee i ee nee 
PEPE < 0, eae Z-70 Negative 0°< < 180° : lead 
jf 
| Pe>0, Oelxe1 | Z_70 | Negative pea << OT sates | 
Ee ee en OE 
% ik Pe> 0, Aye 0 | SF | rotsive Bigs “f<360, || Tag | 
| {| Su SS Se ee a — ae 1 = aie an ia a et —_— : 
| 
y Fe 0. Ae > 0 Ze< 90 | Positive | 05 = ec 180° ! lead | 
j ! : / 
= | Pee ae. ee UU 
; ; ; t { 
| PecO, % 20! ZrO | Negative’ 0 < ¥< 180°; lead 
| eeStO aaa eee | | 


re er ee ee NT RY rN wr a 





Fig. 4-4 shows the pole zero configurations on the S-plane of the 
unacceptable compensator for all possible combinations of the unaccept- 
able solutions P¢ and A, obtained from the single section compensation. 
Thus the results of cascading a single section compensator to the original 
system gives the proper type of compensator to be cascaded. Upon realiz- 
ing that the system cannot be compensated by a single section compensator, 
one may try as a next step to cascade a double section (identical) com- 
pensator in which it is not necessary to know the proper type of compen- 
sator needed, since if the system can be compensated by such a compen- 
sator, then the solutions will now be acceptable and the type of compen- 
sator will be indicated by the solutions. 

Consider the original system in Equation (4-1) cannot be compensated 
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by a single section compensator and an identical double section compen- 


sator is cascaded to give the compensated transfer function 


(s +R 
Cea. Ss eee 4-24 
S)= Biss (s+) baal 


where D(S), P and A are defined as in Equation (4-1) with 
the same design specifications as in Section 4-1. From the steady state 


performance specification the gain K now must be adjusted to give 


ae ene 

Kis daz 

; > a 
or iS = d, Ky A = Kok 


The characteristic equation of the system in Equation (4-24) with K re- 
placed by Ree can be written as 
‘ Pg os ds 7k, 
E(Ss)= DGS) (S4 PY + Rol (S4+PK) =0 
i(s\)=S°D6) +2PS DG) 1 4S u sea Ps4 ce=O0 
or FiG) S}+2PS DS) +P DE6)tNAS 4 2RAPS + KP = (4-25) 


An inspection of Equation (4-25) shows that the Bo vs By curve for 


oF $1 associated with the polynomial F(S) consist of the curves for 


the component polynomials s*D(s), SD(S), D(S) and s?, The curves asso- 


ciated with the polynomial SD(S) and D(S) are computed already in the 


4 1 
single section compensation as the By vs By and BO vs By curves respec- 


2 


tively. The curve associated with the polynomial S“ can be computed 


simply by considering al to be a second order polynomial of the form 
2 

F(S) = S + a,Staq, 

Bg = a 2 - If the curve associated with the polynomial s*D(s) is de- 


for which Mitrovic's equations are By = P2(S)W,: 


noted by the Bo vs oF curve then the Bg vs By curve for $= O41 asso- 
ciated with the compensated system characteristic equation, Equation 
(4-25) is 


22) 






BiH. (9,,@n)+2PB(% art PRE, a4 Ro BSG), 
P : - a (4-26) 

Be=B (5, Qn) +2PB, OF, On) + PB (Sn) 1 Rak Oy, 
which is a curve that depends upon P and 4. (Note that in the single 
section compensation, the compensated system curve depends only on P.) 
The point My which is at a. = G1 on the compensated system curve in 
Equation (4-26) moves if either P or & is varied, and the equation of 
the movement or the locus of the Ma points is obtained by fixing the 
frequency at wD, = Dal in Equation (4-26), that is, the locus of the 


Ma point on the B, vs By plane is 


0 


7 7 ; c fh.6), 3 

B=B (5,0) +2 P BiG. ,Mit PBS ni) + Rok BS. )Ons 

” ) ae Peney (4-27) 
B=BG, jon 2F Gs, Our B, (6, Out) + Kod ae, 

Note that in setting up Equation (4-27) only the first and last terms 


need new calculations. 


1 


The coefficients a, and ae of the terms S~ and 5° in the compen- 


1 
sated system characteristic equation F(S) in Equation (4-25) are 
d)P*+ 2K.X P (The d,P* is obtained from the expansion of the term 
P“D(s)) and KoP? respectively, and the locus of the M° point for the 
compensated system in Eyuation (4-24) is defined by 


B,=diP?+ 2Koa P 
ee (4-28) 
ea sels 
Equations (4-27) and (4-28) can be solved simultaneously for 
and P which are the solutions for the double section compensator. Since 
Equations (4-27) and (4-28) are of second order in P and& , the simul- 


taneous solution of the two equations requires the sclution of a fourth 
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order polynomial which is in general not simple to solve; however there 
are some convenient methods to solve the equation such as graphical 
methods, approximation technique etc. 

If the solutions from double section compensation are not accept- 
able, then a triple section compensation may be tried with the same 
procedure as for the double section compensation; however the results 
require solution of a sixth order polynomial which may be solved by a 
digital computer. 


The following example illustrates the double section compensation. 


Example 4.3 


Given an open loop transfer function 


420 


GoG)= - 





Compensate so that roots are at 
So Uj, = 15 
without reducing the velocity constant K, = 28 
The given system is unstable with the complex conjugate roots in 
the right half of the S-plane which can be seen from the third order 
chart. If a single section compensator is cascaded, then the system 
transfer function becomes 


ROK (STP) 
SCs 71)(st (5 (sp) 


from which the characteristic equation is written as 


G(S)=- 


(64+ 16S 4/5 S*)4P(S24I1GS*4+ IES) +4208 +420P =O 


Soe 1582 gives the Ba vs By curve for 


The polynomial SD(S) = s4 + 16s 
==. 2eas 
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B’ (Sas Bw = 8S sy, 


) 
ae 2 
Bo ($25, Un)= 15 Wy, 16 On 

from which ges =o ca, 15) = -2150 


Bs, =15)=-50625 (4-29) 


are obtzined. 


The polynomial D(S) = $2 + 1687+ 15S gives 
o : 
6, CS=.5, Om ) = 16 Dy, 


Ov an Re ha S 
Bi(S=-5, On )='64,-ay 


to calculate Bris=.5, Op=I5 )= 2¥O 


33($=5, O-15)=225 (4-30) 


Using Equation (4-10) with the values of Bi» By BU» and Bo in 
Equations (4-29) and (4-30) and Kd) = Ky = 420, one obtains the set 
of solutions for the single section compensation, namely 

ie i oh? eo 

Dee a208 
which are not acceptable, thus the original system cannot be compensated 


by a single section compensator. 


The original system is then cascaded with a double section (identi- 


cal) compensator as 


ee 20K" (SHPAk) 
S'S (S445) CS 4p ¥ 





The characteristic equation is then 


= Ss 7 es eee ge Doras S IN one 
BS = IGS 7415S? 4 2p(Sees Ht S)t OPH bSH15E ) -420 0S 


' (4-31) 
43404 PS +420" = O 
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EBs 5) 


in which $“D(S) = S°+ 168 
The Bo vs Bz curve for 5 = .5 associated with the polynomial s2p(s) is, 


noting that the coefficient of the s? texm in s*p(s) is zero, 
B= oa 6 es 5 wnt 


B=-i5a?e +a, 


from which 
B)(é=.5, Wy Hl) =~ 3387S 


(4-32) 
E. ($=5, Wx! )= WOP7TSO 


are obtained, 
Substituting the values in Eyuations (4-32), (4-30) and (4-29) into Equa- 
tion (4-27) the locus of the Ma point for the compensated system with the 


double section compensator is defined by 


B,=-337542P(-3/50) +P*(24e) +4¥20K0°CI5) 


4-33 
B= 7Of750+2 P(-50625)+ e (225) + 420f¢ (225) ( ) 


Using K, = 420 and d,; = 15 in Equation (4-18), the Locus of the M 
points is defined by 
Bj HIS P +P40dk P 


Bo=420 p* 


By solving Equations (4-33) and (4-34) one obtains the set of solutions 


(4-34) 


R= WS 
Aes YOF 22 5.32 
which are the acceptable set. The compensated system transfer function 
is then 2 
J 


ee 700K ae ae 
G(s)= S(S+\(S415)(S 42175) 
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The dominance of the roots at § = .5 and W, = 15 can be checked 
on the third order chart by factoring the characteristic: equation of 
the compensated system into two polynomials (s* + 158+ 225) and F,(S) 
where F,(S) is a third order polynomial, or by sketching the [Ys 
curve of the compensated system as shown in Fig. 4-5. Note that the 
roots at S= .5 and ®, = 15 are not the dominant pair but all other 


roots have & greater than .5. 
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4-3. Successive application of single section compensator design 
procedures. 
When the original system Gg(S) requires a compensator with more 
than one section, the system can be compensated by a non-identical multi- 
section compensator, by simply repeating the single section compensator 
design technique (Section 4-1) with judicious choice of an intermediate 
root location. 
Returning to the equations giving the solutions of single section 


compensation 
By ($1, Qn) 
Ke By G, Ont ) 


R= 
| (4-10) 

ae = | Bi (S,, Onn) + Pe(BPCS 24m) ~ as)] 

where the Mitrovic s variables and the constants are defined 

as in Section 4-1. P¢ and Ag are unacceptable solutions 


associated with the designated root T= -$,Wnt +/i-$? On 


It is obvious in Equation (4-10) that if 5, and/or a) are changed to 


1 
' 
some other values oy » and/or («) i then Equation (4-10) may give some 
n 

acceptable solutions Py and ay and this is equivalent to compensating 
the original system to locate the root at r] ssa + fi-$* Dp 

which is different from the designated root location Ur =-5,Wp + /-3? On 
thus the system is partially compensated by so determined single section 


compensator of P, and a, - To this partially compensated system, the 


1 
procedure of designing a single section compensator with the designated 
root ry is applied and the system may be compensated to the desired 


performance requirements; if not, the procedure is repeated. 


The choice of the intermediate root = or the choice of Sy and 
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Wy, is based on Equation (4-10). Since the only restriction on P is 
P>0Q, the choice of Sy and anes must be in such a direction that the 
numerator and denominator of first of Equation (4-10) have the same 
signs, and yet that P with the other numbers in the second of Equation 
(4-10) must give an oO such that it is positive and not too large and 
not too small. Thus one must consider too many numbers simultaneously 
which is hard to define by a general rule. However, in view of the 
fact that for a smaller ¢ specification the compensation is in general 
0 
easier, by choosing a S smaller than 4%, and analyzing the Ue 
b) 
d 
if 


ni ’ 
curves one can choose a set of suitable 4, and me 


and re 
Suppose S and a, are chosen and a suitable set of Py and oc 


is determined, then the partially compensated system transfer function 


can be written as 


Koo, (St PAY ) 


- 4=35 
DC) S+p) ae 


G(s) = 


where K,/D(S) = GG) is the original system defined as in 
Section 4-1, 


which has the velocity constant at the specified value K, and has a root 


T° To compensate the system in Equation (4-35) so that it has a root 
= —S,Wn, + fie? Oy 4 which is the specified root, a single 


section compensator is cascaded and the transfer function becomes 


Kod d(S+ RA)(St+PAY (4-36) 


6CI=~56)(S4P)(St PD 





The characteristic equation of the system in Equation (4-36) is then 


FOES LS) KP*PISPE)+RPDG) + heal AS Kook (Rb F&)S+ RA P~9 


or FG)=S*DG)+ PSDE)+P(SDE)+ PDE) + KdiAS*+ (KRRAtKAPIS (4-37) 
+E P =O 
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With the Mitrovic's variables for the component polynomials in Equation 
(4-37) defined as in Sections 4-1 and 4-2, the Mitrovic's curve fer 

GS = S1 of the compensated system in Equation (4-36) can be written in 
terms of the component curves, namely 


Bo= BS, +PBNS 9) +P BS, Wn) +B BUS ony) HK oA AES, On 
(4-38) 


4y . Set . fi - 2 
By = = G ? Li, aie BG, On) 4P Bis, 9@,)+ lz BOCs, 4} +a Cn 
By fixing bee : Wnt in Equation (4-38), the locus of the My points 
¢ 
(which is at Oy = Gy On the le curve defined by Equation (4-38)) 
on the By vs By plane is defined by 


B= B)+R 8, + PCB) FRB?) +d XP, ) Wry 
Bo= BRB, 5 Pees th By) + Kot x od Cage) 
where the Mitrovic's variables on the right sides of Equation 
(4-39) are the values of corresponding Mitrovic's equations 
evaluated at 6: S1 and Ww = Ooi 6 
Note that in Equation (4-39) the last terms are from the component 
polynomial Kd. s? in Equation (4-37) as in the case of the identical 
double section compensation in Section 4-2, and if ay =o and Py, = P 
then Equation (4-39) becomes identical to Equation (4-27). Also note 
that in Equation (4-39), the values By> Bo; BY. and Bo are computed at 
the first stag: to reveal P. and ae are unacceptable, By and By have to 
be calculated. 
By collecting the st terms from the component polynomials in Equa- 
tion (4-37) and by doing likewise for the 3° terms the locus of the M° 


point on the By vs B, plane for the compensated system in Equation (4-36) 
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is defined by 


Bedi Pi thot) )P +A X 


4-40 
SSS le . : ) 


Equations (4-39) and (4-40) are linear first order equations which can 

be solved simultaneously without difficulty to give a set of solutions 

in P and 4 . If the solutions are acceptable then the compensation is 
accomplished, and if not then the procedures can be repeated. The follow- 
ing example shows the non-identical double section compensation which is 


seen in Section 4-2 with identical double section compensation. 


Example 4.4 
Consider the problem in Example 4.3 with the same design specifica- 
tions. The original transfer function is 
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Go(S) = SCSt+IDCStHS) 


which cannot be compensated by a single section compensator, The poly- 
nomials to be considered in the design problem are 
3 2 
D(S)= S°+16S +158 
SD6)= S*+ 16 S74 1557 
Jdb)= Ss + 16 Sas 


to give the associated Mitrovic's curves as follows; 


Be( SW») = /b BUS) + BO, 

2 
BE (S,Wn) =/6 By - B Wy 
3 
B, (4,2) =/5B(S)w, +b DS) OS + AE) On 
B, ($, On) = PSG =< 760(2)0, Lio 
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: a a a 
BGM) HIT Ons )Wa 4 16 P.( SI, + PK(E) aye 


} 
i} 


, 3 b abet ate 
B’(4,, a) =7EA(S)Ay 16D. (6) W, - REO 
By cascading a single section compensator to place the root at oe) 


and 4): 15 resulted in an unacceptable set of solutions 


Pi 


which indicates that the system requires a lead compensation. At the 


P. = -259.8 & = -208 
SF 


stage of getting the unacceptable solutions the following are computed; 


} ; a 
B (62.5, WyrI5) = 2UO 


By ($=-5; Mis) Fane ah 


B, ($=-5, a), =15) = = 3/50 
B, (o=-5, a, 2/5) = ~50625 


Since the system cannot be compensated to have the root at = .5 
and 4. = 15 by using a single section compensator, an intermediate root 
is chosen so that the system is partially compensated by an acceptable 
single section compensator. 

Let the intermediate root be at G = .3 and W, = 8 then for this 
root Bi ap By and Bo must be calculated. By using the corresponding 
equations, the following are calculated: 

BY ($5.3, W,2?) = 1/476 
BO(6=.3, W=8) = 7/68 
62.3, @,=")=223.6 


Wy 9°) =-6576.6 
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By substituting those values into Equation (4-10) one obtains 
Pe22.16 
= Gee, 273 

which are acceptable solutions, thus the partially compensated system 


249] S+3.73) 
SCStI)(S +15) ($422.16) 


G(s)= 
has the root at §= .3 and 4) = 8. 
To set the root at $= .5 and 4), = 15 which is required from 
the design specification, another single section compensator is cascaded 


with the partially compensated system to give the system equation 


AS 2 eee ae 
71 SCS 41) (S415) (S422.16) (54+ P) 


The new values to be computed are 
B’(S$=.5, Opals y= so 7) 


Bi ($=.5 , Wn =15) = 708750 


oUt : 
By substituting He Bo, By; BGs By and By for $= .5 and we z= 15 and 
the values of P, and Ay into Equation (4-39), one obtains the locus of 
c P 

the Mg points 

B =~ 23375 + (22,16;(-31F0) + P(-31ISO 422/240) + 249T IEA 

B= 708 15+ (22.16 )(-$062)+P (-$062.5 4.22, 16(2-5) + 249 Tazz )a 
or By=- 73/73 +2168 P+ 37¢55A 


(4-41) 
(3, = 4130/9 —£5539 P+ 561F25 


The locus of the M° points of the doubly compensated system is, by 
using Equation (4-40), 
3, = 2629 P +93 O6BK 
(4-42) 
Bo=9306 Pp 
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Equations (4-41) and (4-42) are solved to give the solutions 
p=2509 
ABM Z2=7.76 

which is an acceptable set. 


The compensated system is 


Gls) = ep Z PERSE IBS 756) 
S$ (S41) (S415 (6422.18 (512504) 


and that the original system is compensated to the design specifications 
is confirmed by the compensated system curve ls and the Mitrovic's 


working point as shown in Fig. 4-6. 
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Fig. 4—6 
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5. Conclusion. 

1) Frequency response. 

It is seen that the frequency response of both closed and open loop 
systemscan be obtained on the Mitrovic's By vs B, plane for any order 
system and particularly for a system with no zeros, the method is simple. 

By using the concept of the constant bandwidth curves on the Bq vs 
By plane a system can be designed so that it meets the bandwidth specifi- 
cation closely when the bandwidth specification is a critical one. Fer 
the third order system, a universal bandwidth chart is constructed which 
can be used in both analysis and design of a third order system. 

2) The fourth order charts. 

Universal fourth order charts on the Bo vS By plane are constructed 
which can be applied in both analysis and design of any fourth order 
system. With the preconstructed fourth order charts, the labor of con- 
structing the necessary curves which is essential in the Mitrovic's 
method is removed or minimized. 

3) Design. ‘ 

An analytic method has been developed designing the required cas- 
caded compensator by using the basic properties of the Mitrovic’s equa- 
tions and the Bo vs B, curves. The method gives accurate solutions of 
the required compensator in order for a system to have its 
roots at arbitrarily designated locations on the S-plane, with arbi- 
trarily designated gain. 

If the given system and the design specifications are such that 
the system can be compensated by a single section compensator, the prob- 


lem is to solve two linear simultaneous equations which provide solutions 
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quickly. For identical multi-section compensation the problem is to 
solve a high order polynomial in which the order of the polynomial de- 
pends upon the sections needed; however, the problem can be guided into 
the ones of solving linear (first order) equations by applying single 


section compensator design technique successively. 
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APPENDIX I 


FUNCTIONS (5) (GENERAL FORMULA AND TABLE) 


PCS) = -C28Pu_ (5) + DO 4(5)) 
Pol5) = 0 
D5) = 


(For the proof of the formula see [1] and[2]) 


Table of 2, (4) functions 


Ke 0 O28 0.2 

















1.0 
ff oa: ea 

a | o ff on ! 0.4 : ce 

i, | io | 0.96 0.84 | 0.64 | 0.360 |-0.44 Tae -3.0 
dy, 9 — -0.392 -0. 736; -0.984) -1.088) -1. 000: -0. 672 -0.056, +4000 


me |. & Se mec td St ei ie = a 
Ps | -1. 000. -0.882 -0. 546 -0. 050, +0. 510. 7“ a base: 038 -5.000 


eS eee ee Oe ae 


Ye | 0 | +0. 568 +0.954) +1. 014 +0. s80) 0 -0.824 -1.398, +6.000 


t | j 
= ee cage ata | ia as 
QD, +1. -000 +0, 768 +0.164) -0. 559 ~1.054; -1. 000; -0. 258) +0. 9181-7, 000 
—_—, - et a Pag OIE +$——}— oa ep eee 7 eve ma er enh wre - et Pewee em a nee 


2g 7 0 | -0.722, -1.020 -0.679: +0. 164) +1. 000) +1. 133] +0. 112 +8..000 











ee oe ee a 2 ae ee aioe he . 
—@g - 1.000, -0.624 +0.244. +0.966 +0.923, 0 | -1,102 -1. 075 -9, 000 
/ { a ee ee ge rr -~—i--— hb. = sn a 


Groh 9 +0. 847 +0.922 +0.099 -0.902 -1.000. +0.189, +1. 393 +10, 000 


rg nn : i i 
ear a ene cmon PO Er ne eA ANA A A th me Sl SNS tes OE oe 


| 
a AROSE A et ET NEPEAN 
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APPENDIX IT 


FORTRAN PROGRAM FOR COMPUTATION OF Bg VS By CURVES 


(Programmed by Dr. J. R. Ward) 
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